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Fundamental differential equations are derived under the unrestricted ap- 
proximation of electrical neutrality that admits trapping. Extension is 
made for applied magnetic field. The transport equations derived hold 
without explicit reference to detailed trapping and recombination statistics. 
Modified ambipolar diffusivity, drift velocity and lifetime function apply 
in the steady state. The same diffusion length is shown to hold for both car- 
riers, and a general " diffusion-length lifetime" is defined. Mass-action 
statistics are considered for cases of (one or) two energy levels. Certain 
" effective" — rather than physically proper — electron and hole capture and 
release frequencies or limes that apply to concentration increments are de- 
fined. Criteria are given for minority-carrier trapping, recombination and 
majority -carrier trapping, and for " shallow" and "deep" traps. Applica- 
tions of the formulation include: the diffusion-length lifetime for the Shock- 
ley-Read electron and hole lifetimes; linear and nonlinear steady-state and 
transient photoconductivity; negative photoconductivity; the photoconductive 
decay observed by Hornbeck and Haynes in p-type silicon; the photomagneto- 
electric effect; and drift of an injected pulse. Photomagnetoelectric current is 
found to be decreased by minority -carrier trapping, through an increase in 
diffusion length. A simple general criterion is given for the local direction 
of drift of a concentration disturbance. With trapping, there may be "re- 
verse drift," whose direction is normally that for the opposite conductivity 
type, and also local regions of carrier depletion that may extend in practice 
over appreciable distances. 
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photoconductivity also is analyzed critically and in some detail, a pro- 
cedure facilitated by results from the present formulation of comparative 
formal simplicity. This analysis involves a formalism that recurs in 
theory of time-dependent transport. A detailed treatment of the drift 
with trapping of an injected pulse is given. These applications of the 
formulation constitute an illustrative selection. 

In Section 1.2 is assembled descriptive material intended to be read 
for further preliminary orientation as to the contents of the paper, and 
also to be read piecemeal with corresponding portions of the main sec- 
tions. 

1.2 Outline of Procedures and Results 

The formulation is accomplished in two stages. By treating concentra- 
tions of added electrons and holes formally as unrelated variables, differ- 
ential equations for the transport are derived in Section 2.1 along the 
lines of previous treatments. 10 ,u t Extension for applied magnetic field 
is included. 11 These equations involve no specific reference to the detailed 
trapping and recombination statistics. Specialized to the steady state, 
the ambipolar continuity equation is formally the no-trapping equation, 
but with the sum of fixed and mobile positive (or negative) charges as 
dependent variable, and with suitably modified ambipolar diffusivity, 
drift velocity, and lifetime function, which depend in general on two 
(concentration -dependent) phenomenological differential "trapping 
ratios." The same diffusion length is shown to apply for both electrons 
and holes, and a general "diffusion-length lifetime," r , based on the 
unmodified ambipolar diffusivity, is defined. The formulation is com- 
pleted in Section 2.2 with equations for the time rates of change of con- 
centrations of carriers trapped in centers of each type. 

These rates are written in accordance with mass action, which provides 
a simple J and general § basis for trapping and recombination. || Two 



t In Ref. 11 small Hall angles are assumed, in part because appreciable mag- 
netoresistance is otherwise involved. As indicated in this reference, arbitrary 
Hall angles (and injection levels) could suitably be taken into account by theory 
involving the phenomenological magnetoresistance without added carriers. 

X See Hoffmann. 2 The mass-action approach, now widely used in semiconductor 
theory, is essentially that used in early theory of metal -semiconductor junctions: 
see Schottky and Spenke. 12 

§ Boltzmann statistics, assumed for the transport equations, imply mass-action 
relationships at equilibrium: see Spenke. 13 But, with definitions of equilibrium 
parameters suitably extended, mass-action equations apply also for degenerate 
semiconductors: see Rose. 14 

|| A treatment based on Fermi statistics that allows for degeneracy and includes 
dependence of occupation probabilities on applied magnetic field has been given 
by Landsberg. 15 
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energy levels (as well as a single one) are considered; equations are 
written in Section 2.2.1 for two types of trapping centers and, through 
simple formal modification, in Section 2.2.2 for two levels from centers of 
a single type.f Partly by way of notational convention, the levels are 
taken as acceptor and donor levels, which give negative and positive 
fixed charges. This case is the simplest for which both steady-state trap- 
ping ratios occur, these ratios being the respective changes in concentra- 
tion of all fixed negative charges and all fixed positive charges divided 
by the change in concentration of all negative or positive charges. With 
suitable interpretation of the notation, the equations apply to one- or 
two-level cases in general; results written for centers of the acceptor 
type, for example, are not restricted to this type. Moreover, it will ap- 
pear that, in the analysis of transient (or steady-state) photoconductivity 
for a given multilevel model, the trapping at a given time need usually 
be considered in detail in no more than two successive levels. Levels 
appreciably lower and higher than these may contribute to recombina- 
tion, but will not contribute to trapping, in the sense that the lower levels 
may be assumed to remain completely full (or else saturated) and the 
higher levels completely empty. J 

To facilitate analysis and interpretations, in Section 2.2.1.2 "effec- 
tive" capture and release frequencies and times that apply to concen- 
tration increments are defined a ■priori from the mass-action equations. 
The four effective frequencies or times for each energy level differ from 
the physically proper ones, which depend on the trapped concentrations 
and thus on the detailed solution of the particular problem. They satisfy 
a fundamental restriction, used extensively in the theory, which is de- 
rived from thermal-equilibrium relationships involving detailed balance. 
With this restriction, quantitative criteria are established in Section 
2.2.1.3 for ranges of minority-carrier trapping, recombination and major- 
ity-carrier trapping. These ranges may be specified in terms of the 
location of the equality level § relative to the Fermi level S F and the 
"reflected Fermi level" &/, the reflection of & F about 8, its location for 
intrinsic material. If spins are taken into account, quantities of the mass- 
action theory serve to locate the trapping level relative to 8. It is shown 

t Theory for multilevel centers is given in Landsberg, 16 ' 16 Champness,' 7 
Okada, 18 Shockley and Last, 19 Mercouroff, 20 Khartsiev, 21 Sah and Shockley, 22 
Bernard, 23 Kalashnikov and Tissen, 24 and Kalashnikov. 26 

t The influence of trapping at a given level on recombination at another has 
been calculated for the near-equilibrium steady state by Kalashnikov. 26 See also 
Mashovets. 27 

§ This is the Fermi level for which the (equilibrium) rates of electron and hole 
capture and release are all equal. 2228 The equality level is similar in purport to 
the demarcation level of Rose, 2930 which is the trapping level for which the rates 
are equal. 
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that a proper criterion for "shallow" or "deep" minority-carrier trapping 
levels is that 8/ separates these levels. Thus, levels in extrinsic material 
considerably shallower than the midgap may still be "deep" levels. 

Detailed theory is given in Section III through various applications 
of the general formulation, and consequences of the mass-action statistics 
are examined. In Section 3.1.1 diffusion length and diffusion-length life- 
time, as well as the trapping ratios, are evaluated from equations written 
for the limiting linear small-signal steady state. A "capture concentra- 
tion" is introduced, use of which is found to simplify formally much of 
the detailed theory, including that for time-dependent cases. This con- 
centration is the concentration of (single-level) centers multiplied by the 
respective equilibrium fractions of centers occupied and unoccupied. 
Values of it that are small or large result, respectively, in negligible cap- 
ture frequencies or in large capture frequencies with negligible release 
frequencies. For the case of a single energy level, the general (equilib- 
rium) Shockley-Read electron and hole lifetimes 31 are obtained in forms 
involving the capture concentration. These lifetimes are shown to cor- 
respond to a diffusion-length lifetime to whose general expression is 
formally the same as that for the common (equilibrium) lifetime 31,32 in 
the limit of small concentration of centers. This common lifetime other- 
wise applies as such only under a condition restricting the capture con- 
centration, which is frequently severe: In the minority-carrier trapping 
range, it is that this concentration be small compared with the equilib- 
rium minority-carrier concentration. From conditions for the neglect of 
quadratic terms in the mass-action equations, the linear approximation 
is shown to imply a restriction of injection level that may be much more 
severe than the familiar small-signal condition 10 based on the conductiv- 
ity change. 

The general single-level trapping ratios and lifetime functions for the 
nonlinear steady state are obtained in Section 3.1.2. These and the 
mobile-carrier concentrations, as well as the volume recombination rate, 
can be expressed in terms of trapped-carrier concentration as single 
concentration variable. The lifetime functions reduce to the Shockley- 
Read lifetimes in the linear small-signal limit and to a single limiting 
large-signal value. The familiar common lifetime function 31 for small 
concentration of centers usually does not apply in the small-signal range 
unless it is substantially constant in this range. The differing general 
lifetime functions otherwise usually apply, and small-signal minority- 
carrier trap saturation obtains. The apparent diffusion-length lifetime 
then increases to a small-signal saturation -range value equal to the 
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majority-carrier release time.f Further increase occurs in the approach 
to a large-signal lifetime, which, in this case, is also the (small-signal) 
lifetime in the limit of strongly extrinsic material of the opposite con- 
ductivity type. Such increases of lifetime can account for certain cases 
of superlinearity, or the more-rapid-than-linear increase of photoconduc- 
tivity with injection level, on the basis of a single trapping level. J 

Transient decay of photoconductivity is analyzed in Section 3.2. In 
the linear small-signal case, the decay is given by a sum of exponential 
modes with (real and positive) decay constants whose number exceeds 
by one the number of types of centers present . 2 ' 28>37 ~ 41 § For nonrecom- 
binative trapping by centers of two types, the decay constants and 
equilibrium concentrations after injection are evaluated in Section 3.2.1 
for electron and hole traps present together and for electron (or hole) 
traps only. With the latter, carriers released from one type may be 
captured in the other. The general linear case for centers of one type, 
including recombination, is analyzed in detail. The two time constants 
are given in forms involving the capture concentration. If one is large 
compared with the other, then the larger may be identified as the lifetime, 
while the smaller represents a trapping transient during which approach 
to the steady-state trapping ratio takes place. This transient has small 
amplitude for small concentration of centers, for which capture rates in 
the ratio of capture frequencies and release rates in the ratio of release 
frequencies decay with the concentration in the lifetime mode. It does 
not occur if the steady-state trapping ratio obtains initially, or if "criti- 
cal recombination" obtains, with which, because of equal capture fre- 
quencies, trapped concentration does not change from initial value zero. 

Sufficiently small capture concentration gives, with the comparatively 
short trapping transient, a lifetime substantially equal to 38 - 39,40 the 
common steady-state electron, hole and diffusion-length lifetime. The 
required condition is frequently severe: In the minority-carrier trapping 
range, it is the same as the common-lifetime condition. Capture con- 
centration large results in decay times equal to 38 - 40 the steady-state 
electron and hole lifetimes and given by the electron and hole capture 
times. If one of these is large compared with the other, then the smaller 
represents the transient for practically complete trapping of the carriers 
of one kind, and the larger represents the recombinative decay of the 

t Approximate steady-state solutions which exhibit small-signal nonlinearity 
have been given by Tolpygo and Rashba. 33 

J A multilevel model for superlinearity has been given by Rose, 29 - 31 ' 35 (and 
Ref. 1, Ch. 1A). See also Bube. 35 

§ See also Ref. 1, Ch. 3A. This chapter also includes some nonlinear cases. 
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carriers in traps and of the carriers of the other kind as these are cap- 
tured. In all these cases, the lifetime decreases monotonically as con- 
centration of centers increases. Under a condition that is usually met 
in the minority-carrier trapping range, this decrease occurs primarily in 
two ranges of concentration of centers, with approximate constancy of 
lifetime in an intermediate range, f 

The photoconductive decay is governed in the general case by non- 
linear differential equations. These are considered for centers of a single 
type in Section 3.2.2. The general single-level problem is rather intractable 
analytically. J Solutions of the nonlinear equations are given for two 
special cases, namely, nonrecombinative trapping and sufficiently small 
concentration of centers or large concentrations of mobile excess carriers 
such that the steady-state lifetimes are substantially equal. The latter 
solution § has the rather restricted general application of the common 
steady-state lifetime function, 31 since it is the integrated form correspond- 
ing to this function. || By solving suitably linearized equations, the decay 
times associated with a small-amplitude pulse of added carriers above a 
steady generation level are evaluated. If, as is often permissible, direct 
recombination may be neglected, 4849 - 50 then the decay in the general 
large-signal limit is exponential with lifetime equal to the steady-state 
large-signal lifetime. During this decay, the concentrations of carriers in 
traps remain constant. This lifetime and the corresponding concentra- 
tions in traps are evaluated for centers of a single type and for the two- 
level cases. A differential equation that is invariant under interchange of 
quantities pertaining to electrons and to holes is derived for centers of a 
single type. It provides a first integral under a condition that holds for 
sufficiently large concentration of centers or concentrations of mobile 
carriers. With this first integral, the decay problem may be formulated 
as a first-order (rather than second-order) nonlinear differential equa- 
tion. The large-signal condition, obtained in this connection, differs from 
the familiar one 10 in that, as a condition for equal electron and hole 
lifetimes, it entails not only relatively large change in conductivity but 

t The approach to constancy with increasing concentration of centers is dis- 
cussed by Wertheim. 41 

| Certain analytical approximations have been considered by Isay. 42 A treat- 
ment which includes numerically computed solutions has been given by Nomura 
and Blakemore. 43 

§ It is equivalent to ones given by Rittner, in Ref. 1, Ch. 3A, and by Guro. 44 

II The decay lifetime has been evaluated as this function by Okada. 46 That the 
nonlinearity according to this function does not account for (small-signal) decay 
in silicon has been observed by Blakemore. 46 This author has fitted dependences 
of lifetime on injection level and temperature assuming two-level recombination 
from one type of center or from two types. The common lifetime function has been 
employed for centers in germanium by Iglitsyn, Kontsevoi and Sidorov. 47 It 
appears that these centers were in the recombination range. 
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also requires saturation of centers, which may be present in relatively 
large concentration. 

Conditions are obtained for centers that give recombination with 
substantially constant lifetime for minority carriers and inappreciable 
trapping. Constant lifetime that applies in the small-signal range also 
applies in the large-signal range, provided the energy level of the centers 
is not too far from the Fermi level towards the majority-carrier band. 
It requires, however, sufficient strongly extrinsic material. In material 
of mixed conductivity type, the recombination rate cannot, in general, 
be specified in terms of a minority-carrier lifetime. But "linear re- 
combination" may apply, characterized by a two-lifetime recombination 
rate that is the sum of contributions respectively proportional to the 
added minority- and majority-carrier concentrations. The assumption of 
general linear recombination is also a convenient notational device: In 
the analysis of models involving nonrecombinative traps in conjunction 
with the recombination centers, it permits deriving results in forms that 
apply for any conductivity of either type. 

The phenomenon of negative photoconductivity, or the decrease in 
conductivity below the equilibrium value upon optical injection, f re- 
sults essentially from excitation of minority carriers from traps with 
recombination in other centers. Theory for this effect is given in Section 
3.2.3, a general expression for mobile-carrier concentrations being derived 
for the linear small-signal case. This result is of comparative formal 
simplicity and shows that the effect tends to be offset by recombination 
in the traps and to be enhanced with deep traps of small capture cross 
section. J 

A general procedure is outlined in Section 3.2.4 for analysis of trapping 
models with a number of discrete energy levels, which relates the various 
decay times to capture cross sections and these energy levels. This 
procedure is applied to observations of Hornbeck and Haynes 54 on elec- 
tron trapping in p-type silicon. § For the sample on which the most 
extensive measurements were made, the decay times ranged from 20 
microseconds to 260 seconds. Their model, that of two kinds of non- 
recombinative traps with recombination in other centers, is found to 
imply a hole-capture cross section of the deep traps and of the shallower 



f This has been analyzed by Stockmann. 51 It has recently been observed in 
silicon by Collins. 52 Infrared quenching of photoconductivity or luminescence 
from short-wavelength excitation, discussed by Rose 29 - 3435 and others, is a closely 

related effect. 

t Excitations involving trapping levels may increase normal photoconduc- 
tivity. 53 

§ See also Ref. 1, Ch. 3F. 
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(but still "deep") traps that is small compared with about 10~ 24 and 
10~ 20 cm 2 , respectively. 

These cross sections are calculated from expressions for the decay 
times for nearly empty traps. Recombination in the deep traps cannot 
account for the observed decay: The hole-capture cross section that 
gives the decay time of 260 seconds for the nearly empty traps would 
give a considerably larger decay time, rather than the observed value 
of 1 second, for the traps nearly full. Recombination in the shallower 
traps, however, can account for the observed decay: The hole-capture 
cross section calculated from the decay time for these nearly empty traps 
is in close agreement with that which fits the entire decay in the deep 
traps. The lifetime with traps filled of 20 microseconds may then be 
ascribed to recombination in the higher level of two-level shallower 
traps. A recalculation of cross sections and energy levels on the basis of 
this model gives hole-capture cross sections large compared with 1.2 X 
10 -17 cm 2 , equal to 2.4 X 10 -20 cm 2 and small compared with 10~ 24 cm 2 , 
respectively, for the recombination level and the shallower and deep 
trapping levels. The corresponding electron-capture cross sections f 
are 2.3 X 10~ 15 , 1.1 X lO" 13 and 2.9 X 10" 14 cm 2 , the last two being 
half an order of magnitude smaller than the ones calculated by Horn- 
beck and Ilaynes. The shallower and deep trapping levels are found to 
lie 0.007 ev above and 0.23 ev below the Fermi level for intrinsic ma- 
terial; the latter trapping level is 0.78 ev below the conduction band. J 
Use is made of the observed straggle effect from the shallower traps, 
comparison being made with the theoretical expression derived in 
Section 3.4.3 for the limiting decay time at fixed location for the tail of 
the distribution from a pulse injected under applied field after the 
maximum has drifted past. It is shown that a model for which the 
trapping levels are levels of centers of a single type cannot account for 
the observations. While the levels found are close to two levels of 
gold, 56 - 57 it is thus unlikely that they result from a single metallic im- 
purity. This conclusion bears on the indications that the deep traps are 
associated with the presence of oxygen as an impurity. 

The steady-state photomagnetoelectric (PME) effect with trapping 
is analyzed in Section 3.3 on the basis of the general formulation with 
applied magnetic field. § Equations formally similar to those for no trap- 
ping apply in terms of redefined quantities that involve the trapping 

t Theory to account for such large cross sections has been given by Lax. 56 
j An energy gap of 1.10 ev at 300°K is used rather than 1.00 ev as in Ref. 54. 
§ A treatment of photoconductance and PME voltage with trapping under ac 
illumination is included in: Lashkarev, Rashba, Romanov and Demidenko. 53 
Mironov 59 deals with the transient decays after removal of steady illumination. 
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ratios. The effect may exhibit small-signal nonlinearity with nonuniform 
lifetime if recombinative deep traps in the minority-carrier trapping 
range are involved. The influence of trapping as such is investigated. 
It is found that nonrecombinative minority-carrier traps (in conjunction 
with recombination centers) increase diffusion-length lifetime by an 
amount proportional to the capture concentration, f Thus, minority- 
carrier trapping decreases PME current. A comparatively slight de- 
crease in r and increase in PME current results from majority-carrier 
trapping. 

Detailed illustrative procedures and related theory are given for the 
determination of capture cross sections, concentrations and energy 
levels from suitable PME and photoconductivity measurements at given 
temperature. With trapping in recombinative traps of a single type, the 
PME current-conductance ratio involves light intensity implicitly 
through its dependence on the lifetime t c that is defined in terms of the 
change in conductivity for a given steady, uniform volume-generation 
rate. The ratio, however, determines a relationship between r and t c , 
a transcendental relationship with the preferred method 11 of the high- 
recombination-vclocity dark surface. This relationship, in conjunction 
with suitable additional conductance measurements also independent of 
light intensity, suffices to determine both t and t c , then the capture 
(and release) frequencies and capture concentration, and finally the 
quantities sought. The linear small-signal theory! is given for recombina- 
tive traps of a single type and also for nonrecombinative traps with 
recombination centers, for which the results are essentially similar 
though somewhat simpler. 

Preliminary to analysis of transport problems, the general ambipolar 
continuity equation is specialized to the linear small-signal case in 
Section 3.4.1. Then, for trapping (and recombination) in centers of a 

t Jonscher 6 " gives an increase of diffusion length with trap concentration which 
is bounded and always essentially negligible, a result at variance with that given 
here. In Jonscher's nonambipolar treatment, the continuity equation does not 
include a term in the second .space derivative of trapped-carrier concentration. 
Though this term is relatively small for sufficiently strongly extrinsic material, 
its neglect significantly affect's the higher-order differential equation for concen- 
tration of mobile minority carriers, obtained by eliminating trapped-carrier con- 
centration, in that it gives a coefficient of the term in the second space derivative 
that is too small by just the factor by which diffusion-length lifetime is increased. 

t Zittcr" 1 discusses the phenomenological dependence for any model of electron 
and hole lifetimes on r,- and lifetime derived from the PME effect (in the thick 
slab). The latter is the same as r , and Zitter relates it to a diffusion length. 
Amith, 82 - 63 has presented the effect of nonrecombinative traps on the PME cur- 
rent-conductance ratio, and has pointed out that the predominant effect is usually 
on conductance. That on PME current is generally negligible in comparison if the 
traps are minority-carrier traps and are present in not too large concentration in 
sufficiently strongly extrinsic material. 
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single type and linear recombination in other centers, the respective 
concentrations are shown to satisfy certain third-order partial differential 
equations, which are of the second order in time. These reduce to the 
same equation if there is no volume generation. They otherwise each 
contain a term proportional to the volume generation function, the 
equations for the mobile-carrier concentrations containing the time 
derivative of this function as well. 

The case analyzed in Section 3.4.2, that of injection into a filament 
in the steady state with applied field, yields qualitative information of 
interest. If a certain frequency v v , the "straggle constant," is positive, 
then the field-opposing and field-aiding solutions in the regions separated 
by the point of injection are sharply varying and gradually varying 
exponentials, as in the no-trapping case. 64 But with negative v v , grad- 
ually varying field-opposing and sharply varying field-aiding solutions 
obtain. In the limit of no diffusion, these give added carrier concentra- 
tions only in the direction opposite to the direction of drift normally 
determined by conductivity type. This "reverse drift" is explained by a 
simple and entirely general criterion, obtained from the fundamental 
equations, for the local direction of drift of a concentration disturbance: 
Normal or reverse drift occurs according to whether injection results in 
proportionately more or fewer minority carriers than there are at thermal 
equilibrium. For no trapping, for example, the concentrations are in- 
creased locally by the same increments, so that proportionately more 
minority carriers result if the material is extrinsic; and zero drift 10 ' 64 
obtains if the material is intrinsic. Conditions for the sign of v v are given. 
It is shown from these that reverse drift, which occurs for sufficiently 
large trap concentration in not too strongly extrinsic material, occurs 
with nonrecombinative trapping if minority carriers are trapped so that 
the fraction of the time they are free is smaller than the equilibrium 
minority-carrier to majority-carrier concentration ratio. 

Drift of a pulse of carriers injected into a filament, with trapping by 
centers of a single type,f is analyzed in detail in Section 3.4.3. Bilateral 



t Fan 37 - 05 has given a solution of this drift problem which applies for negligible 
majority-carrier capture frequency. Clarke 40 has, in effect, pointed out this re- 
striction, to which solutions for the decay of photoconductivity given by Fan 37 
and Rittner 1 are also subject. Jonscher 66 has given solutions for drift of minority 
carriers with recombination and nonrecombinative trapping at variance with solu- 
tions given here. The otherwise plausible neglect by Jonscher, in a nonambipolar 
treatment for strongly extrinsic material, of a term in the continuity equation 
involving the gradient of trapped-carrier concentration is apparently not justified. 
In the differential equation for concentration of mobile minority carriers, it re- 
sults in minority-carrier release frequency only as a factor in the concentration- 
gradient term instead of v v , which, for this case, is substantially the sum of the 
capture and release frequencies. This neglect of the capture frequency is tanta- 
mount to neglect of the capture concentration compared to the equilibrium minor- 
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or two-sided Laplace transforms derived in Appendix A are used to 
obtain solutions of the differential equation for negligible diffusion. These 
solutions are of two main types, according to whether a frequency unit, 
v, connected with the introduction of dimensionless variables and 
parameters, is real or imaginary. From theory concerning v and the 
parameters, it is found that real v implies either the minority-carrier 
trapping range or the recombination range. f Illustrative solutions of 
real v for which nonrecombinative trapping is assumed are presented 
graphically. These show that carriers that remain untrapped appear in a 
comparatively rapidly attenuated pulse that drifts at the ambipolar 
velocity. This remnant of the initial pulse leads a continuous distribu- 
tion, which, as a result of multiple trapping, ultimately spreads as a 
time-dependent gaussian distribution and exhibits a maximum that 
drifts at a fraction of the ambipolar velocity. For nonrecombinative 
trapping, this fraction approaches comparatively slowly a limiting value 
that does not exceed the fraction of the time the carriers are free, J and 
(for imaginary v as well) the fraction of carriers trapped, obtained by 
integrating over the drift range, approaches comparatively rapidly the 
fraction of the time carriers are trapped. Recombination in other centers 
reduces the distance for a maximum at given time and thus the ap- 
parent mobility of the distribution. The decay constant for the straggle 
effect is found to be the (positive) straggle constant, accordingly so 
named. 

Imaginary v obtains over the majority-carrier trapping range and, for 
nonrecombinative trapping, over the reverse-drift range. With re- 
combinative trapping, it obtains also for zero drift and over a normal- 
drift range other than that of majority-carrier trapping. Illustrative 
solutions, for which nonrecombinative trapping is assumed, are presented 
graphically for reverse drift and for majority-carrier trapping. It appears 
that in the reverse-drift range an attenuated pulse of untrapped carriers, 
which drifts at the ambipolar velocity, leads a continuous distribution 

ity-carrier concentration. It presumably leads, for example, to the conclusion of 
this reference that a verv short pulse is transmitted without distortion and is only 
attenuated. Also, the solution given for the steady state of continuous injection 
should properly include the fraction of the time minority carriers are free as a 
factor in the exponent. That is, trapping results in a more gradual decay with 
distance; for this case of no diffusion, a lifetime applies that is equal to the sum of 
the lifetime proper and the (generally much larger) lifetime for multiple trapping, 
which is discussed in Section 3.2.4. 

t As shown in Section 2.2.1.2, these ranges in their entirety together constitute 
the "minoritv-carrier capture range," for which the equilibrium minority- to 
majority-carrier capture frequency ratio exceeds unity. It is shown in Section 
3.4.3 that there is a minority-carrier capture range of imaginary v which includes 
the reverse-drift range. m . 

X Fan" has shown from his solution that this limiting value is, for relatively 
small trap concentration, equal to the free-time fraction. 
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of added mobile minority carriers, which crowds towards the injection 
point as its maximum excursions both above and below the axis increase 
with time. There is local carrier depletion, the distribution being nega- 
tive over part of the drift range after a certain time. The distribution 
approaches a pulse at the injection point of strength equal, for non- 
recombinative trapping, to the initial strength times the free-time 
fraction. It does not exhibit essentially unidirectional drift: The drift of 
added carriers, initially in the direction of the ambipolar velocity, is 
largely in the opposite direction after some trapping has taken place. 
A numerical estimate of the effect of diffusion indicates that negative 
added-carrier concentrations can occur over appreciable distances under 
conditions that can be realized in practice, f The illustrative solution 
for majority-carrier trapping shows that negative added-carrier con- 
centrations occur in this case also. Majority-carrier trapping, however, 
results essentially in drift at the ambipolar velocity and, if it is nonre- 
combinative, the fraction of carriers trapped approaches the trapped- 
time fraction. 

The solution is given for "critical trapping," the borderline case be- 
tween cases of real and imaginary v. For nonrecombinative trapping, it 
is the same as that for zero drift and gives exponential continuous dis- 
tributions that are established progressively as the drift range in- 
creases and otherwise do not change with time. For trapping in intrinsic; 
material without diffusion, drift does not start, and the initial pulse re- 
sults simply in pulses for the concentration increments that remain at 
the injection point, where they change as trapping and recombination 
proceed. With diffusion, ambipolar drift occurs, since the condition for 
zero drift no longer holds in the intrinsic material as carriers are trapped 
away from the injection point. Further physical interpretations for the 
various cases are obtained by evaluating the current density of added 
carriers, which represents the equal departures for given total current 
density of the electron and hole flow densities from their values for no 
added carriers. 

1.3 List of Symbols 

The following list includes most of the symbols to be employed, and 
is largely consistent with previous notation. 10 • u - 64 

t Kaiser 67 has suggested that negative added-earrier concentrations that 
were observed with localized optical injection in silicon under applied field may 
be accounted for through these results. A theoretical discussion of carrier deple- 
tion is included in Ref. 10. 
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a = parameter in distribution of (153). 
A n j , A p j = capture cross sections for electrons and holes for the 
jth energy level or type of center. 
b = iin/np , drift mobility ratio. 
d = coefficient for direct electron-hole recombination given 
in (37). 
C„j , C pj = electron and hole capture coefficients for the jth energy 
level or type of center. 
D = kTn„v P (n + p)/c = (n + p)/(n/D p + p/D„), ambi- 

polar diffusivity for no trapping. 
D' = modified ambipolar diffusivity, defined in (31). 
Di = diffusivity for intrinsic material. 
D„ , D,, = diffusion constants for electrons and holes. 
J) n t J) p = diffusivities defined in (133). 

Do — value of D at thermal equilibrium. 
Do = value of D' at thermal equilibrium. 
e = electronic charge. 
E = electrostatic field. 
S = Fermi level for intrinsic material. 
& F = Fermi level. 

8/ = "reflected Fermi level," the reflection of 8 P about 8. 
Sj = electron energy for the jth. energy level or type of 
center. 
F n , F n , F p = fractions of mobile electrons, trapped electrons, and 
mobile holes for drift of a pulse, given by (164). 
g = rate of volume generation of electron-hole pairs. 

Aff = g - go. 

g a = value of g at thermal equilibrium. 
G = quantity defined by (9). 
A(7 = conductance increase of slab per unit width, given by 

(116), (118) and (129). 
Go = dark conductance of slab per unit width. 

I = I p + I» , total current density. 
Al = current density of added carriers, defined by (19). 
I = diffusion current density, defined by (7). 
I„ , I p = electron and hole current densities. 
h, h = modified Bessel functions, in the notation of Watson. 
J , J l = Bessel functions, in the notation of Watson. 
h = Boltzmann's constant. 
k = unit vector in the direction of magnetic field. 
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K = quantity defined by (10). 
K a = factor in (129) for conductance change of illuminated 

slab, evaluated in (130). 
K T = factor defined in (119) by which diffusion-length life- 
time with nonrecombinative trapping exceeds re- 
combination lifetime. 
L = vot, length unit defined in (143). 
£ = surface rate of generation of electron-hole pairs from 

strongly absorbed radiation. 
L = (D to)*, diffusion length. 
$ = operator symbol for two-sided Laplace transform. 
m=p-\-p = n-\-n. 
Am = m — Wo . 

ma = value of m at thermal equilibrium. 
n = electron concentration. 
An = n — no . 

An = two-sided Laplace transform of An. 
AN = An/((P/L). 
AN = An/((P/L). 

n = concentration of fixed negative charges. 
An = n — no . 

An = two-sided Laplace transform of An. 
AN = An/((P/L). 
91, 9ly = total concentrations of centers. 

m = electron concentration for the Fermi level at the j'th 
trappmg level. 
Uj* = equality density, defined in (54). 
3Xy* = "capture concentration," defined (forj = 1) in (63). 
n 8 = n — po- 
no — value of n at thermal equilibrium. 
n = value of n at thermal equilibrium. 
Ni , N 2 = dimensionless decay constants defined by (147). 
p = hole concentration. 
Ap = p — p . 

Ap = two-sided Laplace transform of Ap. 
AP = Ap/(<?/L). 
AP = Ap/((?/L). 

p = concentration of fixed positive charges. 
Ap = p — po . 
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(P = number per unit area of carrier pairs injected over 

cross section of filament. 
pj = hole concentration for the Fermi level at the jth trap- 
ping level. 
p s * = equality density, defined in (54). 
p = value of p at thermal equilibrium. 
fto = value of ft at thermal equilibrium. 
r„ = dn/dm, steady-state trapping ratio. 
r nj , r„j = trapping ratios for transient photoconductive decay 
modes, evaluated in Section 3.2.1; also similar quan- 
tities that are given in (149). 
r p = dft/dm, steady-state trapping ratio. 
fll„, , (R„ , (ftp = functions specifying rates of decrease of m, n and p 
through trapping and recombination, 
s = Laplace transform variable. 
s„i , s n , s p = surface recombination velocities for m, n and p, re- 
lated by (115). 
T = temperature in degrees absolute. 
U = t/r, dimensionless time variable. 
v = ambipolar drift velocity, defined by (8). 
v' = modified ambipolar drift velocity, defined in (31 ). 
v„ , v p = velocities defined in (133). 

Vu = value of v at thermal equilibrium. 
v ' = value of v' at thermal equilibrium. 
V = electrostatic potential. 
X = x/L, dimensionless distance. 
on , a-> = quantities defined in (78). 
<*io = n /ni = pi/po . 
«20 = Po/pi = n-i/no . 
Ai = quantity defined by (62). 
f = dimensionless parameter defined in (145). 
rj = dimensionless parameter defined by (152). 
6 m d p - B H = d p 4- | 9„ | . 
d n , d p = Hall angles for electrons and holes. 
G = variable defined by (157). 
k = dimensionless parameter defined in (145). 
H„ , n P = drift mobilities for electrons and holes. 

v = frequency unit (real or imaginary) defined in (143). 
v u = quantity defined in (136). 
v uni , v opj — "effective" release frequencies for electrons and holes, 
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defined for acceptor-type (j = 1) and donor-type 
(j = 2) centers in (50), and for two-level centers in 
Section 2.2.2. 
Vij = decay constants defined by (47). 
j/ na , j/ p3 = decay constants for 'linear recombination," defined in 
Section 3.2.2. 
v t — sum of hole and electron capture and release frequencies 
defined by (83). 
v i»j , v tp j = "effective" capture frequencies for electrons and holes, 
defined for acceptor-type (j = 1) and donor-type 
(j = 2) centers in (50), and for two-level centers in 
Section 2.2.2. 
v v = "straggle constant," defined in (136). 
Vl , v> = decay constants for photoconductivity. 
£ = dimensionless parameter defined by (151). 
<r = a,, + c7 p , total conductivity. 
a n , <T P = partial conductivities for electrons and holes. 
t = time unit defined in (143). 
Tc = conductivity lifetime, defined by (117). 
Tgnj , T 0P j = "effective" release times for electrons and holes, the 
reciprocals of v gn j , v gp j . 
r„, = steady-state lifetime or lifetime function for Aw. 
f,„ = lifetime for Am with recombination centers as well as 
traps, evaluated in Section 3.3. 
t„ , Tp = steady-state electron and hole lifetimes or lifetime 
functions. 
Tno , Tp = (Cj3ly) -1 , (Cpfilj)' 1 for a particular j. 

Tr = apparent lifetime from FME current-conductance ratio, 
evaluated in Section 3.3. 
Tt„j , Tipj = "effective" capture times for electrons and holes, the 
reciprocals of v tll j , v lp j . 
7 = "diffusion-length lifetime," evaluated in (35) and (65). 
fo = "diffusion-length lifetime" with nonrecombinative traps 
and recombination centers, given by (119). 
Tl , to = time constants for photoconductive decay, the re- 
ciprocals of vi , v% . 
t 3 = lifetime for decay through recombination centers, the 

reciprocal of v n3 or v p z . 
To, — photoconductive decay time for nearly empty traps. 
* = -e -, 8. 
^ = potential defined by (15). 



CURRENT-CARRIER TRANSPORT AND PHOTOCONDUCTIVITY 533 

II. GENERAL FORMULATION 

2.1 The Transport Equations 

The general neutrality condition may be written as 

m = p + p = n + n, (1) 

which states that the total concentration of positive charges, the sum 
of the concentrations of mobile holes and fixed positive charges, is equal 
to the corresponding total concentration of negative charges. It is, as 
will appear, of advantage to deal with the total concentration m of 
charges of either kind. 

By way of extension of the familiar (nonambipolar) continuity equa- 
tions for holes and for electrons that apply for no trapping, two forms 
of the continuity equation for m may be written: 

dm/dt = dp/dl 4- dp/dt = -e~ l div I p + g - (ft,,, 

- dn/dt + dn/dt = e~ l div I B + g - (ft,,, . (2) 

Here, for simplicity, the same volume-generation-rate function g is as- 
sumed for both holes and electrons; generalization to include excitation 
to or from trapping levels (as well as interband excitation) is given in 
Section 2.2.3. The volume rate (R m is associated with trapping and 
recombination. It depends directly on the various concentrations and 
not explicitly on coordinates and time; dp/dt and dn/dt contribute only 
to (Jim, and, if these are respectively subtracted from (2), then con- 
tinuity equations for holes and electrons, namely, 

dp/dt = -e~' div I p + g - (ft,, , (ft p = (ft,,, + dp/dt, 
dn/dt = e~ l div l a + g - (ft,, , (ft,, = (ft,,, + dn/dt 

result. The same volume rate (R,, is properly used in each of the equa- 
tions (2) since it depends directly only on concentrations; it must 
apply, in particular, in the case of zero current densities. As (2) shows, 
this use of ol,„ is consistent with the neutrality condition and with the 

ion 

div 1 = 0, I = I p + I„ , (4) 

which applies in regions containing no sources or sinks of current. 
Differing volume rates for p + p and n + h are properly introduced 
only if there is appreciable space charge. 
The familiar current -density equations, 
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I„ = o-pE — eD p grad p 

and (5) 

I„ = o-„E + eD n grad n, 

apply under the assumption of Boltzmann statistics, which imply also 
proportionality of the hole and electron mobilities /x p and n„ to the 
corresponding diffusion constants D p and D„ in accordance with Ein- 
stein's relation. By use of these equations and the neutrality condition, 
(1), a continuity equation for m of ambipolar form may be derived: 
The hole and electron current densities I„ and I„ are eliminated from 
(2) ; the electrostatic field E is also eliminated by means of the expression 
for E involving the total current density I that is obtained by adding 
the equations in (5) ; and use is made of (4), the condition of solenoidal 
I. This procedure is similar to that previously employed in the no-trap- 
ping case 10 except that, for the required generality, p and n are treated 
formally as unrelated variables. The single continuity equation for m 
that results from (2) may be written in various forms as follows :t 

dm/dt - g + (R m = -e _1 div I D - v-G 
= -e _1 (divI D + K-I) 
= -e" 1 div I*, (6) 

where 
I D = —e<T~\(TpD n grad n + a n D„ grad p) 

= —e[D grad m — a~ (a p D n grad n + <r n D p grad p)] 
= —eD(n + p)~ grad np = — ekTn n ^ P a~ grad np, (7) 

v = efi„n P (n - p)a~% (8) 

G = (n — p)~\?i grad p — p grad n) 

= grad m — (n — p) -1 (n grad p — p grad n), (9) 

K = (Z>„ — D p y l grad D = e'n„n P (T 2 {n grad p — p grad n) 
- cUp" 1 - iO" 1 grad [a~\n + p)] 

= -ein- 1 + n.-V grad [^(n - p)}, (10) 

and 

I* ■ &' + 01", a + = 1, (11) 

t This equation specialized to the case of Ap = can be shown to be consistent 
with a continuity equation for Ap derived by Rittner 1 under the assumption of 
a common lifetime function for electrons and holes. 
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with 
I' = l D - e ( Mp _1 + M«" , )~ 1 o-~ 1 (w - p)I 

= -ein- 1 + M-'^^Kn - p)E + (fc!T/e) grad (n + p)] 

= (6+ D-'CHp -I.) (12) 

I" = l D + e ( Mj ,- 1 - ^-^"'^(n + p)I 

= e( M p _1 - M»~T'[(« + P)E + (*r/«) grad (n - p)] 

= (6 - lr^Mp + In). 

The diffusivity D, a definition of which is contained in the last of the 
equations in (7), is the general ambipolar concentration-dependent 
diffusivity, which occurs in the theory for no trapping. It is used here 
simply for notational convenience. The velocity v of (8) is properly 
interpreted in the continuity equation as two velocities, nv/(w — p) 
for drift of Ap and v v /(V ~ n ) for drift of An. It is otherwise formally 
similar to the ambipolar velocity of the theory for no trapping, except 
that now n — p is not a constant concentration. The first two right- 
hand forms of the continuity equation, (6), exhibit terms associated 
with diffusion and drift, respectively, as comparison with the continuity 
equation for the no-trapping case shows. The current density I*, as given 
by (11), is introduced for generality; with solenoidal I, the divergences 
of I' and I" are equal, t From the expressions for these current densities 
in terms of I p , I„ and the drift mobility ratio 6 given in ( 12), it may be 
verified that I* may be chosen as I p or — 1„ , as in (2). Indeed, as is 
otherwise evident, I* may be written simply as a linear combination of 
I p and — 1„ , normalized as in (11), since a linear combination, so nor- 
malized, of any two I* is also an I*. The current densities I' and I" are 
introduced because their use is frequently convenient. 

The mobile-carrier concentrations n and p are, in accordance with 
(1), properly written as m — n and m - p where they occur explicitly 
and in the diffusivity D, in the electron and hole conductivities a n = 
en„n and a p = e-itpp and in the total conductivity a = a n + o> . For ac- 
ceptor and donor centers of single types, dn/dt and dp/dt in terms of 
the various concentrations provide, with the continuity equation, three 
simultaneous differential equations in the dependent variable m, n and 
p. For more than single types of acceptor and donor centers, n and p 
are sums of fixed-charge concentrations. Equations are then written for 
the rates of increase of each of these concentrations, and the number of 



t Note that these divergences equal that of I„ and (-I«); also, I" - I' equals 
2b(b 3 - 1)-'I. 
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simultaneous differential equations exceeds by one the total number of 
types of centers present. These are the differential equations for the 
general transport problem with trapping and recombination provided I 
is a known function of the space coordinates and time. 

In some cases, I must be determined from boundary conditions. Use 
is then made of the fundamental differential equation, (4), which may 
be written to involve the electrostatic potential V as additional de- 
pendent variable. With 

o-E = -a grad V = I - e grad (D n n - D p p), (13) 

it follows that (4) may be written in the form 

div [a grad V - (kT/e) grad (<r„ - a p )] = 0, (14) 

in which <r„ , a„ and a are to be expressed in terms of m, n and p. In this 
formulation, V is introduced into the continuity equation through the 
elimination of I by means of (13). Another procedure, of advantage in 
some connections, involves use of the potential 

¥ = V - (kT/e)(b - 1)(6 + l)" 1 In (<x/<x ) (15) 

instead of V as dependent variable. Then I is given by 

I = —o- grad ^ + eDi grad(p — n), (16) 

where /),- = 2(D P ~ 1 + D„ -1 ) -1 is the diffusivity in intrinsic material. 
Aside from the effect of trapping on I, as given by the second term of 
(16), SP is the potential that "drives" the total current density. This 
may be described as the electrostatic potential modified by the Dember 
potential. The latter gives the field associated with diffusion of carriers 
of differing diffusion constants. 
Electrostatic field given by 

E = e-\n + V r\L,/nn + I P /fh, ~ kT grad(p - n)] f (17) 

an equation somewhat analogous to (16), is a result obtained by solving 
for E in the equations for I" in (12). As (17) shows, E in the absence of 
trapping (and of appreciable space charge) may be written in a form 
that does not involve concentration gradients explicitly. 

In the ambipolar form of the present treatment, the equations of (5) 
are 

Ip = (ffpA)I + h> = (<r„o/<ro)I + AI 

and (18) 

I» = (<r»/<01 -Id = (<TnoAo)I - AI, 
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in which zero subscripts denote values at thermal equilibrium; AI is 
defined by 

AI = e 2 n„n P <To~ 1 (T~\nop — pm)l + I/> 

= Wo^U + Mp)(wo + Po)I' + (n„ - /t P )(no - ]h)l"} 



-i, 



= (To [Clip — (Tpfll,,). (19) 

Application of one or more of (13) or (17) and (18) is frequently re- 
quired in connection with boundary conditions. The ambipolar dif- 
fusion current density I,, includes the effect of the Dember field and 
contributes the same particle flow density to both I„ and I p . Use of the 
expressions for I„ and I,, that involve AI is of particular advantage for 
physical interpretations and in small-signal cases, since AI is the current 
density of excess mobile carriers. 10 For given total current density I, 
it represents the equal electron and hole flow densities, that are the 
departures from the thermal-equilibrium flow densities and that do 
not contribute to I. Note that I* may also be chosen as AI. 

2.1.1 Extension for Applied Magnetic Field 

The current densities for Hall angles 6 n and d p small are given in 
general by Equations (10) and (13) of a previous paper. 11 These result 

in 

e(dm/dt - g + <R«) = -div I„ = div I„ 

= — div(apE) + eD p div grad p 

- 0„lgrad a p , E, k] 

= div(o-„E) + eD n div grad n 

+ Mgrad <r» . E, k}, (20) 

in which k is a unit vector in the direction of the magnetic field and the 
heavy brackets denote scalar triple products. With n and p treated 
formally as unrelated variables, multiplying respectively by a„ and <r p , 
adding and simplifying gives 

dm/dl - g + (R*. = -c _I (div l„ + K-I) 

- e 2 M„M„o-~ :t [0!(MpP 2 S racl n + P* 71 * grad p),I,k]| 

+ {9 p a n + d n a p )(D n n - D p p)fgrad n, grad p,k}], (21) 

where I n and K are defined in (7) and (10) and 6 is the sum of the 
magnitudes of the Hall angles, d„ — 9„ . In deriving this continuity 
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equation, use is made of curl E = 0, which holds for steady applied 
magnetic field; time dependence of I generally has a quite negligible 
effect. 11 Use is also made of the relationships 

o-E = I — e grad (D„n — D p p) 

- <T~\d p a p + d n a„)l X k - 6l D X k (22) 

and 

l(8 p n grad p - 6 n p grad n),E,kl = 

a~ l {(0 P n grad p - d n p grad n),I,k] (23) 

+ e(T\d p D n n - d n D p p)lgrad n, grad p,k}, 

which hold with the neglect of terms quadratic in Hall angles. Equation 
(23) is obtained in a straightforward manner from (22), which is ob- 
tained by writing total current density as 

I = exE + e grad(D„n - D p p) + (0 p a p + 0„<r n )E X k 

(24) 
- e grnd(d p D p p - n D„n) X k 

and then solving for E. The terms on the right-hand side of (24) repre- 
sent, respectively, drift, Dember, Hall and PME contributions. 

A differential equation that expresses the solenoidal property of I is, 
from (24), 

div[«rE + (JcT/e) grad((r„ - * r )] 

(25) 
+ fgrad(V P + e n a n ),E,k} = 0. 

If direct use must be made of this fundamental equation, then it is well 
to eliminate I from (21) by means of (24), and to employ the electro- 
static potential V as one of the dependent variables. 
The current densities are given in ambipolar form by 

l P = (*,/*)! + I- 

and (26) 

In = (o-„/<r)I — I-, 

where, if terms quadratic in Hall angles are neglected, 

1= a I D + e(a n a p /<T 2 )l X k - <r _1 (V» + n <r p )I D X k. (27) 

Components of total current density perpendicular to the applied 
magnetic field are 
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I, = <jE x + e i- (D„n - Z) p p) + o- _1 ( Vp + ^h)7„ + 01*,, (28) 

and /„ , which is given by a similar expression obtained by inter- 
changing x's and ?/s and (to retain a right-handed coordinate system by 
effectively reversing the direction of the z axis) changing the signs of 
the Hall angles. One way of deriving (28) is to substitute the expression 
obtained by solving for E y in the equation for l y , for E y in the equation 
for I x obtained from (24), and to neglect terms quadratic in Hall angles. 

2.1.2 Formulation for the Steady State in Terms of Trapping Ratios 

A number of results for the steady state can be established from the 
general differential equations without specifying in detail the trapping 
and recombination statistics. Differential "trapping ratios" 

r„ = dn/dm, r p = dp /dm (29) 

are introduced. These apply since, in the steady state, n and p each 
depend directly only on total concentration m of negative or positive 
charges. In the immediate context, r n and r p will be considered simply 
as factors that depend in general on m, which, multiplying grad m, give 
grad n and grad p, respectively. They apply, of course, for any number 
of types of centers present. Their evaluation for particular models is 
given in Section 3.1 in connection with the more detailed analysis of 
the steady state. 

With (29), it follows from (6) through (9) that the continuity equa- 
tion for the steady state may be written as 

div (D' grad Am) — v'-grad Am + Ag — Am/r m = 0, (30) 

in which D' and v' are modified ambipolar diffusivity and drift velocity 
that are given by 

D' = kTn n y, p <j-\(\ - r p )n + (1 - r n )p] 

= [1 - (r p n + r n p)/(n + p)}D, 

v' = en nl x p a~ 2 [(l - r p )n - (1 - r n )p]I 

= [I - (r p n - r n p)/(n - p)]v, 

and in which the net generation rate g — (R m has been written as the 
increment in this rate over thermal equilibrium, Ag — Aw/t„, , with 
Ag and Am being the corresponding increments in g and m and t„, a 
lifetime function for Am. The modified diffusivity and velocity do not 



540 THE BELL SYSTEM TECHNICAL JOURNAL, MAY 1960 

apply to time-dependent cases; v' would, for example, give the effect of 
applied field on apparent diffusion length, but is not, as will appear in 
Section 3.4.3, drift velocity for an injected pulse. 

Expressing the concentration gradients for the steady state in terms 
of r„ or r p and grad Am also formally simplifies (14), the differential 
equation that must also be used if the current-flow geometry is not 
known. In connection with (18), the current-density equations, the 
procedure results in ambipolar diffusion current density given by 

I„ = -eD' grad Am (32) 

for the steady state. 

The trapping ratios defined by (29) can assume negative as well as 
positive values: If centers of a given type trap mostly carriers of the 
opposite charge, then a negative trapping ratio obtains. Consider, for 
example, trapping in centers of the acceptor type, which are neutral or 
negatively charged. For these, positive r„ cannot exceed unity; it nearly 
equals unity if electron trapping is the predominant process, so that the 
excess trapped electron and mobile hole concentrations are substantially 
equal. If, however, hole trapping is the predominant process, then r„ 
is a large negative number, the increment in concentration of fixed 
negative charges being negative and balanced by the excess mobile 
electron concentration, so that m retains substantially its thermal-equi- 
librium value. Similar considerations apply to r p for centers of the donor 
type. Thus, the trapping ratio is close to unity or a large negative num- 
ber according to whether the centers predominantly trap carriers of the 
same charge or of the opposite charge. t 

For a large negative trapping ratio, the comparatively small incre- 
ments in m are associated with large magnitudes of D' and v', as (31) 
shows. A concentration variable other than Am may then be more suit- 
able. The equation in the linear combination A An + BAp (with constant 
A and B) of the excess mobile-carrier concentrations that results from 
(30) has diffusivity and velocity equal to D' and v' each divided by 
A (I - r») + B(l - r p ), since, from (29), grad Am is [A(l - r n ) + 
B( 1 — r p )] -1 grad (A An + BAp) . In general, they are bounded in magni- 
tude for all values of the trapping ratios that can occur. J In this equa- 

t For acceptor centers, say, of total concentration 3I ( , the trapping ratio 
r p = d(3li — ii)/d(p + 9li — n) for holes may be defined. The two ratios are 
symmetrically related: They may be interchanged in r,,' = r „/(/•„ — 1) and, as 
one increases to unity, the other becomes negatively infinite. 

t Note that A(l - r„) + B(\ - r P ) equals r p - r„ for A = 1 and B = -1 and, 
if r„ and r p are constants (as obtains under suitable small-signal restriction), also 
for A = r„ and B = —r„. It follows that n — p and r p n — r n p are (under this re- 
striction) both subject to diffusivity D '/(r p — /„) and velocity Vo'/Op — ''")• 
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tion, the recombination term is properly written as —(A An + BAp)/ 
(Ar„ + Bt p ), where t„ and t p are lifetime functions for Aw and Ap. 
These functions are respectively equal to An/(R m and Ap/(R m since, from 
(3), (R„ = (R p = (R m holds for the steady state. 

The equilibrium lifetimes for electrons and holes differ in general, but 
are nevertheless always associated with the same diffusion length. This 
result follows readily from (30), whose linear small-signal form is 

D Q ' div grad Am — v '-grad Am + Ag — Am/r m = 0, (33) 

the zero subscripts denoting thermal -equilibrium values, f The lifetime 
function t,„ is here constant; and, since An and Ap equal (1 — r,,)A?n 
and (1 — r p )A?n, with r„ and r p the thermal -equilibrium trapping ratios, 
(33) implies 

(1 - r„) -1 Do' div grad An - (1 - r„)"V-grad An 

+ A0 - An/(1 - ?'„)r m = 

for electrons and a similar equation for holes. Thus, for An the lifetime 
is r,„ multiplied by (1 - r„), while — as may be established in greater 
generality from (29) and (30) —the diffusivity and velocity are those 
for Am multiplied by the reciprocal of this factor, and similarly for Ap. 
It follows, in particular, that the product of equilibrium diffusivity and 
lifetime, which is the square of L , the diffusion length, is the same for 
An, Ap and Am, independently of the particular trapping and recombina- 
tion statistics. 9 A "diffusion-length lifetime" r , based on the unmodi- 
fied ambipolar diffusivity D , may accordingly be defined by" 1 

To = La /D = (Do'/Do)T m = [1 - (r„no + r„p,,)/(rtu + Po)]r„, 

(35) 

= ("oTy, + PoT, t )/(n + p ), 

in wliich t p and r„ are the equilibrium lifetimes for Ap and An. The more 
detailed analysis of Section 3.1.1 includes evaluation of the single diffu- 
sion length and lifetime t that correspond to the (equilibrium) Shock- 
ley-Read electron and hole lifetimes. Diffusion-length lifetime for re- 
combination in the presence of nonrecombinative traps is evaluated in 
Section 3.3. 

For the steady-state formulation that includes applied magnetic field, 
it is readily shown that (21), the continuity equation, assumes the form 
of (30) if v' is redefined in accordance with 



t As shown in Section 3.1.1. the required small-signal restriction may be more 
severe than that given in Ref. 10 for the no-trapping case. 
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v' = en„n P a~ 2 {[(l - r p )n - (1 - r„)p]l 



(36) 



+ 0(e/a)[(l - r p )ix n n z + (1 - r n )n P p]I X k). 

Note that the second scalar product in (21) vanishes, since steady-state 
concentration gradients are collinear vectors. For this case, use of the 
trapping ratios formally simplifies (25), as well as (27) and (28), which 
involve the form for I D of (32). 

2.2 Mass- Action Theory 

2.2.1 Single-Level 'Centers of Two Types 

In this section, centers of both the acceptor and donor types are as- 
sumed to be present, namely centers that can have respectively single 
negative or positive charges or be neutral. By use of a suitable conven- 
tion, the equations apply, in effect, to the more general model of two 
types of centers each of which has two states of charge (which differ 
by one electronic charge). On the basis of equations of this section, 
theory for centers of a single type but with two energy levels or three 
states of charge is given in Section 2.2.2. 

Under the assumption of mass-action interactions, the equations 

g — Sim = g — Cinp — C pi \pn — ?>i(9Ii - to)] 

- C„2[np - to 2 (91 2 - p)], 

dfi/dt = (R„ — (R,„ = CitnOli — to) — TOito] 

(37) 

— C P i[pn - pi(9li - to)], 

dp/dt = a p - (R m = -C„2[np - to 2 (91 2 - p)\ 

+ C p2 [p(9l 2 - p) — p 2 p] 

hold. The first equation gives (R m , and it (as well as the other two) is 
obtained by considering the photoconductive case of uniform concen- 
tration and no transport, g — (R m being the contribution to dm/dt that 
does not involve transport. Four processes are taken into account for 
each type of center. In the second equation, for example, the term C p \pn 
is the volume rate of neutralization of fixed negative charges by holes; 
C p i is a phenomenological capture coefficient, t which depends in general 
on temperature and not on concentration. The second term in the same 
brackets gives the rate for the inverse process, C p ipi being the emission 



t In the terminology of Sah and Shockley 22 this quantity is called a capture 
probability. 
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coefficient for hole emission from a neutral acceptor center. Here 9li is 
the total concentration of the acceptor centers, and the concentration 
pi , constant at given temperature, is defined by the condition that the 
quantity in brackets vanishes at thermal equilibrium, in accordance with 
detailed balance. The preceding brackets relate to the interactions of 
the same centers with electrons, the term C nl n (9li — n) being the volume 
rate of capture of electrons by the neutral acceptor centers and C n \n\ 
being the coefficient for electron emission from the charged ones. The 
concentrations Wi and pi are those of the Hall-Shockley-Read theory, 
and are here introduced without explicit reference to Boltzmann statis- 
tics. 2 ,13,14 The third equation expresses the dependence of dp/dt on the 
analogous processes for the donor centers. In the first equation, which 
includes the rate drip of direct electron-hole recombination, only inter- 
actions that change the total concentration m are involved. 

The sign or magnitude of the charge that a center can assume is not 
of material significance in the analysis of this section ; although written 
symmetrically for fixed charges of both signs, (37) may formally be 
transformed so as to apply to two types of donor or acceptor centers. 
This possibility is related to the circumstance that the fixed charges are 
not properly considered as trapped carriers, since the trapping processes 
are manifest through changes in fixed-charge concentrations rather than 
in these concentrations themselves. For example, centers of the acceptor 
type function as electron or hole traps according to whether the concen- 
tration of the charged centers increases or decreases with carrier injec- 
tion. Consistent with the discussion in Section 2.1.2 of the steady-state 
trapping ratios, either type of center may be considered alternatively as 
an electron trap or a hole trap, under the convention that a change in 
fixed-charge concentration resulting from trapping may be negative as 
well as positive. To establish this result from (37), write the two equa- 
tions that apply for, say, acceptor centers only. Then transform these so 
that the concentration 9li — n of neutral centers becomes concentration 
of fixed positive charges p, and the concentration n of charged centers 
becomes concentration 3L> — P of neutral centers ; note that a given in- 
crease in the original n is equivalent to the corresponding decrease in 
the new p. New equations then result that (with the replacement of 
C„\ and C p i respectively by C n2 and C p2 ) are the ones that follow directly 
from (37) for donor centers only. 

2.2.1.1 Thermal-Equilibrium Relationships. The definitions 

rii = n (9li - n ) /n , Pi = Pon /(Vli - n ), 

(38) 
n 2 = Wopo/(9l2 - po), p-2 = Po(9I 2 - po)/po 
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are required by detailed balance. It is evident from these equations that 
nipi = n^po = ?io7>o — nf (39) 

holds, where n, is the thermal-equilibrium electron or hole concentration 
in intrinsic material. Note that (39) states, in effect, that the product 
— (C,,ini)(Cpipi) or {C .aTio) {C p ipo) — of the electron and hole emission 
coefficients equals n,- 2 times the product of the corresponding capture 

rr- ■ . 2, 22.2S 

coefficients. 

If the concentrations in the right-hand members of (38) are known, 
then ni , pi , w 2 and p* are, of course, determined. Certain relationships 
hold between the concentrations. Since n — po — n s equals po — ?io 
from (1), the neutrality condition, this condition and the last equation 
of (39) give 

Wo = §{[(po - no)' + W]* + (po - wo)} 

and (40) 

Po - £{[(po - no) 2 + 4/i/p - (po - no)}. 

It is readily found from (38) that fractions of charged acceptor and donor 
centers are given respectively by 

rio/9li = (1 + aio" 1 )" 1 , Po/9l 2 = (1 + o*," 1 ) -1 , (41) 

with 



aio = no/ni = pi/po = |wi [(n, + 4n 4 - ) 5 + n s ], 
o- 20 = p /p 2 = n 2 /n = 5P2 _1 [(w s 2 + 4n?)* — n t ]; 



(42) 



the final expressions on the right follow by use of (40). For given semi- 
conductor material at given temperature, n, is known and n g is deter- 
mined by conductivity type and conductivity <r , and ni and n 2 (and 
hence pi and p* also) are accordingly determined by the fractions of 
charged centers. Expressing the thermal-equilibrium concentrations of 
mobile carriers and fixed charges in terms of each other (with other 
concentrations as parameters) thus involves roots of quadratic equations. 
The relationships given apply regardless of the number of kinds of centers 
present, since (40) contains no quantities pertaining to particular cen- 
ters, and each equation of (41) and (42) contains quantities pertaining 
only to a single kind of center. 

On the other hand, the fixed-charge and mobile-carrier concentrations 
for centers of two kinds are obtainable in general in terms of 3li , 312 , n-i , 
Wo and w,- 2 . It will suffice to indicate that the concentrations are roots of 
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(43) 



biquadratics that follow readily from the equations 

m = n Q + WiWo/CKi — «o) = po + Pi(9li — W«)/«o , 
mo = no + w 2 (9l2 — po)/po = po + P2po/(9l2 — po) 

for no and p and 

m = a w [th + 9l]/(l + an,)] - aso[p2 + 912/(1 + a M )], , x 

(44) 

ai0«20 = »2/«l = Pi/ Pi 

for aio and a 50 , and hence n and p . Equations (43) and (44) are ob- 
tained from (38) by eliminating, respectively, the concentrations of 
mobile carriers and fixed charges by use of the neutrality condition, t 
They are equivalent to combining (40) with (41) and (42), which are 
accordingly subject to a requirement of mutual consistency. For exam- 
ple, temperature determines n? for a given semiconductor; specifying 
conductivity also then determines n and p ; specifying further w x and 
n-i determines n /9li and po/9l2 from (41) and (42), but only one of 91 1 
and 9l 2 can now be independently specified, since p — n a must equal 

n — pu • 

Through familiar considerations involving equilibrium Boltzmann 
statistics, the concentration m or pi (and n 2 or p 2 ) has been shown to 
equal electron concentration in the conduction band or hole concentra- 
tion in the valence band for the Fermi level coincident with the energy 
level of the centers.'" The relationship 

n, = nr/ Pi = n ie ^- &,kT = n* M +* lVkrl (45) 

for acceptor centers is here employed, and a similar one is used for donor 
centers. Here 4q = — e _1 £i and * = — e _I S are the equivalent electro- 
static potentials of the energy level £i of the centers and the Fermi energy 
8 for intrinsic material. This relationship is more phenomenological than 
those involving the energies of the conduction-band and valence-band 
edges and which give >h and pi in units of the effective densities of states 
in the bands. Note that the temperature dependence of the energy gap 
is involved through n, , while the difference between the effective densi- 
ties of states or the effective masses with nonspherical energy surfaces 
in momentum space is reflected simply in a difference between <J> and the 
midgap potential. If statistical weights associated with spin degeneracy 
arc taken into account, then the definitions of (38) are of course retained, 
but (45) is modified. The right-hand members (for Wi) are multiplied 

f II is easily soon that rubies result for centers of one kind only, or if complete 
ionization obtains for one of two kinds of centers. 



(46) 
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by two ; the exponentials for pi are multiplied by one-half. In the similar 
result for donor centers, the exponentials for n^ and p 2 are multiplied 
by one-half and two, respectively. For given ni and n 2 , these modifica- 
tions ' produce comparatively minor changes in 81 and 82 or \Fi 
and ^2- 

2.2.1.2 Equations in Concentration Increments; Trapping and Release 
Frequencies and Times. For detailed analysis, it is advantageous to re- 
place (37) by equations in the increments Aw, An and Ap in m, n and p 
over their thermal-equilibrium values. By subtracting from (37) the cor- 
responding thermal-equilibrium equations, in which the time derivatives 
and the quantities in the various square brackets are zero, the result 

Ag — A(R m = Ag + v n Am + v\. 2 An + vnAp — dAnAp 

— CpiApAn — CoAnAp, 
dAn/dt = A(R„ — A(R,„ = v 2i Am + v 22 An -f v 23 Ap 

— (C, a An 4- C p iAp)An, 
dAp/dt = A(R P — A(R m = v-nAm + v n An + v 3i Ap 

— (C n? An 4- C p2 Ap)Ap 

follows, in which the decay constants of what will be referred to as the 
"vij notation" are given by 

Mi = —Ci(n + po) — CpiWo — C„ 2 po , 

"12 = Cip — C P i(po+ pi) + Cn2p0 , 

V13 = C.Wo — C n2 (n + w 2 ) + C pl n , 

"21 = C„i(9li — wo) — C pl n , 

f22 = -C„i(9li - no + rc + n x ) — C r i(p + Pi), 

viz = C P ifiQ , (47) 

"31 = C r2 (Vl2 — po) — C n2 p y 

V 32 — C n2 po , 

v n = -C n2 (n Q + n 2 ) - C p2 (% 2 - po + Po + P2). 

Zero subscripts denote thermal equilibrium values. Note that A(R„, , Afll„ 
and A(R P are respectively (R,„ , (R„ and 6i p minus g = dni . In (46), in 
which Am, An and Ap are to be considered as dependent variables, the 
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quadratic terms have been written compactly with An and Ap, which 
may be replaced by Am — An and Aw — Ap; and n and po in (47) may 
be replaced by m — n and m — po . 

It is desirable to supplement the v t j notation with another notation, 
which, although it often results in less compact expressions, facilitates 
physical interpretations. If volume generation and direct recombination 
are neglected for the present, the respective contributions to dAn/dl and 
dAp/dt other than the terms involving transport processes as such may 
be written as 

— A(R„ = — A(R m — dAn/dt = — itmAn + v gn iAn 

- v t niAn + jv, 2 A(91 2 — P). 
-A(R P = -A(R m - dAp/dt = - v lpi Ap + j> op iA(91i - n) 

- V tp2 Ap + VgpiAJ). 



(48) 



The top and bottom rows of the forms on the right give the respective 
contributions of the acceptor and donor centers. The decay constants 
may be identified as certain capture and release frequencies by com- 
parison with the equations 



— A(R„ = — C„i[(9li — n ) An — n An — AnAn] + C n i>iiAn 

- C„?[poAn + n Ap + AnAp] -f- C n2 n 2 A(9l 2 - p), 
-A(R P = - C P i[n Ap + p An + ApAn] + Cpip J A(9l 1 - n) 

- C p2 [(Vl 2 - po)Ap - p Ap - ApAp] + C p2 p 2 Ap, 



(49) 



which follow from (37). In (49), the magnitudes of the contributions 
involving brackets are capture rates, while the other terms on the right 
are release rates. 

Expression of the capture rates in terms of capture frequencies would 
require writing them with An or Ap as a factor, and would thus necessi- 
tate solution of the particular problem. These physical capture frequen- 
cies would depend in general on coordinates and time. The contributions 
to the capture rates that contain An and Ap as factors are associated, 
however, with trap saturation: These contributions, for carriers of given 
charge, represent the decreases and increases in capture rate with the 
filling of centers that assume, respectively, the same and the opposite 
charges. They may, in a phenomenological sense, be deleted from the 
capture rates and assigned to the release rates. The "effective" capture 
and release rates that result from this procedure are clearly rates in 
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terms of which the capture and release frequencies of (48) and related 
times may be defined as follows: 

Electron capture by neutral acceptors: 

Vmi = r7»i = C„i(3Ii - n ) = C„i9li/(1 + a w ), 

Electron release from charged acceptors: 

Vgnl = Tjnl = C„l(n + 111), 

Hole capture by charged acceptors: 

Vtpi = rjpi = CpiWo = C P i9liaio/(l + «io), 
Hole release from neutral acceptors: 

Vgpl = ^1 = C pl (p + pt), (50) 

Electron capture by charged donors: 

Vtn2 = Un2 = C n 2p0 = C n2 9l 2 0! 2 o/( 1 + « 2 o) , 

Electron release from neutral donors: 

V tl n1 = Tg„-> = C n 2(n + 7l 2 ), 

Hole capture by neutral donors: 

Vt P 2 = rjii = C' p2 (9l2 - p ) = C p2 9l2/(1 + a 20 ), 
Hole release from charged donors: 

V IIP 2 = Tgio = C P 2(P + Pi). 

The second forms given for capture frequencies follow by use of (41). 
Note, for example, that vt„i is an average frequency per electron of 
electron capture by a neutral acceptor center and hence the reciprocal 
of the corresponding electron capture or trapping time, r ( „i ; and that 
Vg H i is an average frequency per charged center of electron release from 
a charged acceptor center and hence the reciprocal of the corresponding 
electron release time, t„„i . The saturation terms that originate from the 
true capture rates appear as the contributions from n and p in the 
"effective" release frequencies, while the "effective" capture frequencies 
do not depend on the injection level. 

It is readily seen that, if direct recombination is neglected, then an 
alternative procedure! for including the quadratic terms in (46) is to 

t Another alternative procedure is to replace n u by n and f>o by p or, more 
generally, by increasing /In and p by a fraction 7 of An and of Ap and n and pn 
by the fraction 1 — 7 of An and of Ap. The definitions of (50), which correspond 
to 7 = 0, are then modified, and depend on 7, the fraction of the quadratic terms 
assigned to capture. 
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generalize the v,j by replacing no and p in (47) by n and p. The formal 
equivalence of (46) and (48) that results then permits expressing these 
generalized va in terms of the effective capture and release frequencies 
of (50), as follows: 

"U = — V tv i — V t n2 , 

"12 = — Vgpl + Vln2 , 

Via = Vtpl — V n2 , 

"21 = I'M — Vtpl , 

"22 = —Vt„i - Vgm - v ap i, (51) 

vn = vtpi , 

"31 = —Vt„2 + Vtpi 

"32 = Vtn'l , 

^33 = — V gn o — V tp 2 — Vgpl . 



Note that 



VU + "23 + "32 = (52) 



holds for this case of no direct recombination. 

The four effective trapping and release times or frequencies for each 
type of center satisfy a fundamental restriction, namely: 

rgnJT, pJ _ VtnjVgpj _ p 1 + Ap/(p + Pj) ■ __ , „ /„>, 

TtnjTgpj Vg n jVtpj 7l 1 + An/ {U + tlj) 

Thus, only three are independent. As will appear, this restriction is 
widely useful for calculations and physical interpretations. It is essen- 
tially a consequence of detailed balance: For thermal equilibrium, it 
follows readily from relationships tantamount to this principle, such as 
(41) and (42) or the definitions of (38). The factor on the right that 
depends on An and Ap results simply from the concentration dependence 
of the effective release frequencies. 

A property easily established from (42) and (50) is the following: 
For centers of given capture coefficients and energy level, if the electron 
and hole capture frequencies are equal for a given conductivity, then 
the equilibrium release frequencies are equal for material of the opposite 
conductivity type and the same value of | no — p \, that is, for material 
such that the values of no and po are, in effect, interchanged. 
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2.2.1.3 Trapping and Recombination Ranges; Shallow and Deep Traps. 
Three linear small-signal ranges, characterized respectively primarily by 
minority-carrier trapping, recombination and majority-carrier trapping, 
may be defined for each type of center by use of (53). The minority- 
carrier trapping range" is defined by the condition that the equilibrium 
minority-carrier to majority-carrier release frequency ratio exceeds unity. 
In p-type material, this ratio, v gH j/v gp j , is C„yW;/C PJ po = C n jno/C P jPj , 
from (39), (41), (42) and (50); and, from (53), v tn j/p t pj is larger by 
the factor p /n . The "majority-carrier trapping range" is defined by the 
condition that the majority- to minority-carrier capture frequency ratio 
exceeds unity, for which the equilibrium majority- to minority-carrier 
release frequency ratio is larger by the factor po/wo for p-type material, 
or by Wo/po for n-type. The "recombination range" is defined as that not 
included in either trapping range. Thus, the recombination range is 
given by n /wy = Pj/po ^ C n j/C P j ^ py/no = po/nj for p-type material, 
the electron-trapping range by C„j/C P j > pj/no = po/nj , and the hole- 
trapping range by C n j/C p j < Pj/po = n /nj . A "minority-carrier capture 
range", which includes the trapping and recombination ranges, may be 
defined by v t nj/v t pi > 1. Similar results, obtainable by interchanging n 
and p, hold for n-type material. Ranges of minority-carrier-dominated 
and majority-carrier-dominated transitions ' are those parts of the 
trapping ranges here considered for which strong inequalities hold. 
Equal capture frequencies, which occur at the boundary between the 
recombination and majority-carrier trapping ranges, result in what will 
be termed "critical recombination", with which, as will be seen, An or 
Ap is zero. 

The three ranges may be specified in terms of the equality densities. 
These are the equilibrium carrier concentrations for the Fermi level 
coincident with the equality level. They are defined in the present 
context by 

nf = CpjPj/Cnj = PQVtpj/vtnj = n V gpj /Vgnj 

and (54) 

Pi* = C nj nj/C p j = n V t nj/vtpj = pOVgnj/pgp} , 

in which the release frequencies are equilibrium values. Thus, the re- 
combination range is given by n ^ pj* ^ pa or p ^ nf ^ no for p- 
type material, the electron-trapping range by n* < n or p* > p and 
the hole-trapping range by n } * > p or pf < n , and similarly for n- 
type material. The ranges may evidently also be specified in terms of 
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the equality level, the Fermi level 8 for intrinsic material, the actual 
Fermi level S F and the "reflected Fermi level," 8/ = 2& - & F , the re- 
flection of 8p about 8: For the recombination range, the equality level 
is between & F and S F ' ; for the minority-carrier trapping range, it is be- 
tween Sf and the edge of the majority-carrier band; and for the major- 
ity-carrier trapping range, it is between 8/ and the edge of the minority- 
carrier band. Note that, if the capture coefficients are equal, then n,* = 
Pi (or p* = n,) holds and the respective trapping ranges are given by 
conditions on the trapping level &, obtained by interchanging those on 
the equality level. 

The volume rates of electron and of hole transitions at equilibrium 
are respectively C„in (9li — wo) = C„inin = n v t ni and CpiPono = 
C P iPi(9li - no) = poftpi for acceptor-type centers. From (53), these 
rates are proportional tot v g „i and v gp i . Hence each definition given for 
a trapping range insures that the transition rate at equilibrium for the 
particular carriers is the larger, and also that the transition rate v tni or 
v tP i per mobile carrier is the larger too. The asymmetrical relationship 
between the definitions for minority- and majority -carrier trapping re- 
flects the circumstance that a transition rate will be the larger if either 
the cross section or the concentration of the particular carriers is suffi- 
ciently large. The recombination range is that for which a larger transi- 
tion rate per mobile minority carrier is associated with a total transition 
rate for majority carriers which is the larger. 

For shallow minority-carrier traps, since relatively few are occupied 
by minority carriers at equilibrium so that they can capture majority 
carriers, the condition for the minority-carrier trapping range may be 
met even though the capture coefficients are comparable in magnitude. 
For deep traps, since relatively few can capture minority carriers, the 
minority-carrier trapping generally requires a minority-carrier capture 
coefficient considerably the larger. Suitable condition for "shallow" and 
"deep" traps are, in view of the condition on C nj /C P j for the electron- 
trapping range, respectively p, « n (or n, » p ) and n, « p (or 
Pj » n ) in p-type material. That is, "shallow" and "deep" traps for 
minority carriers are appreciably removed from the reflected Fermi 
level 8/, towards the edges of the minority- or majority-carrier band. 
Similarly, for majority-carrier trapping, "shallow" and "deep" traps are 
appreciably removed from the Fermi level & F , towards the edges of the 
majority- or minority-carrier band, respectively. 

t They equal v onl and v apl times the capture concentration, as shown by (63) 
in Section 3.1.1. 
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2.2.2 Centers with Two Energy Levels 

The formalism for centers of two types is readily modified to yield 
equations for one type of center with two energy levels. With the assump- 
tion that the centers can each assume single negative or positive charge 
or be neutral, n and p denote concentrations of centers in the respective 
charged states. It is thus clear that the fundamental mass-action equa- 
tions for this case are formally the same as (37), with the modification 
that both 3li — n and 9l 2 — p are replaced by 91 — n — p, where 31 is the 
total concentration of the centers. 

For thermal equilibrium, definitions of n t , pi , n 2 and p 2 apply that 
are equations of (38) with bothOli — n and 9l 2 — po replaced by 91 — 
rip — i)n . Tt follows that the restriction 

?i 2 / n t n 2 = PiPi/p* = no/po (55) 

holds for this two-level case. As is easily verified, (39) and (40) still 
apply, while the fractions of charged centers are 

«o/9l = (1 + ni/n -f ttifto/fto 2 ) -1 

= (1 ■ + Po/pi + po 2 /piPi)~ 1 

= aio/(l + ajo + «2o), 

_ x (56) 

po/9l = (1 + p 2 /po + P1P2/P0 2 ) 

= (1 + n /n 2 + no 2 /^iW 2 ) -1 

= a 20 /(l + «io + a 2 o), 

with aio and aw given by (42). The modifications of (43) and (44) for 
the biquadratics are the replacement of 9lj — no and 9l 2 — po by 91 — 
no — pa andof 9li/(l + am) and9l 2 (l + a 20 ) by 91/(1 + a w + a 20 ). 
Note that (55) is not an independent equation, in that it is implicit in 
the modified (43). f Relationships formally identical with (45) give rii 
and n 2 in terms of the two energy levels. 

The fundamental two-level mass-action equations for no direct re- 
combination yield the equations in concentration increments 

Ag — A(R m = Ag — (v tP i + v^Am + {v tn2 — v gpl — C„ 2 ?i 2 )An 

+ ("*pi - Vgm - C r ipi)Ap, 

dAn/dt = (vmi — v tp i)Am — (v lnl + v „i + v opl )An 

(57) 
+ (vtpi - Can — C pl pi)Ap, 

dAp/dt = (v lp2 — v ln2 )Am + (vmi — C p2 p — C„ 2 n 2 )An 

- (ftp-2 + Vgpl + V an2 )Ap, 

t The corresponding restriction for two types of centers has (fto/po) (91.2 — po)/ 
(iJli —no) as right-hand member. 
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from which appropriate v,j can immediately be identified. Effective cap- 
ture and release frequencies are here employed whose definitions are 
provided by (50) if 91 ! — n a and 3l 2 — p are replaced by 91 — /) — p , 
and9Ii/(l + a l0 ) and9l 2 /(l + a 20 ) replaced by9l 2 /(l + mo + a 2 o). 
Aside from these modified definitions, the equations of (57) are formally 
identical with equations for two kinds of centers except for the additional 
terms in which the capture coefficients appear. These are "constraint" 
terms. The ones in BAn/dt represent the decrease in this rate that results 
from the decrease in the concentration of neutral centers associated 
with an increase in p; neutral centers capturing electrons and emitting 
holes are the two processes that increase n. The rate decrease Ci/iAp 
is that associated with the electron capture, while C,,ipAp is that asso- 
ciated with the hole emission. The condition that the rate decreases for 
these two processes be the same is clearly n t * = n. Similarly, the con- 
straint terms in dAp/dt represent the respective decreases (7 p2 pAn and 
C„ 2 n 2 A?l in the neutral -center hole capture and electron emission rates 
associated with an increase in n; these decreases are equal if p>* = p 
holds. The third forms of (5-i) show that a pair of equal constraint terms 
implies an equilibrium hole-to-electron release-frequency or transition- 
rate ratio for the acceptor or donor levels equal, respectively, to n/no 
or po/p, which are substantially unity near thermal equilibrium. 

For this two-level case, the four effective trapping and release times 
or frequencies associated with each energy level satisfy the fundamental 
restriction that is formally identical with (53). It is also easily verified 
that the various conditions given for the recombination and trapping 
ranges and for shallow and deep traps apply without formal modifica- 
tion. 

By suitable notational generalization of the fundamental mass-action 
equations, the results of this section can be shown to apply to two-level 
centers in general, whose states (differing successively by one electronic 
charge) may include ones that are multiply charged, either positively 
or negatively. Through use of the phenomenological capture coefficients, 
statistical weights associated with multiply charged states do not enter 
explicitly. For example, the concentration 91 — h — p of neutral centers 
may be replaced by concentration p of centers with single positive charge, 
and p used to denote concentration of centers with double positive 
charge. Then h is replaced by the new concentration 91 — p — p of neu- 
tral centers.! Thus, with obvious modifications in the physical descrip- 

t Note that these transformations applied to (1) give Aw = An — Ap — Ap = 
Ap + Ap. While the correct neutrality condition holds, Aw is no longer the incre- 
ment in concentration of total negative or positive charges. 
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tion of capture and release frequencies and other quantities, the theory 
is essentially unchanged. 

2.2.3 Volume Generation with Excitations Involving Trapping Levels 

The excitations associated with the absorption of radiation of wave- 
lengths beyond the limit of intrinsic absorption may be taken into ac- 
count phenomenologically through suitable generation terms in the 
differential equations. To the volume rate g of interband excitations in 
the differential equation for n is added g ci + g C 2 , where g c i is the volume 
rate of hole excitations from the conduction band to centers of type 1 
— that is, electron excitations from these centers to the conduction 
band — and g C 2 is the similar quantity for centers of type 2. Similarly, 
in the differential equation for p, to g is added g vX + g* , each term of 
which is the volume rate of electron excitations from the valence band 
to the centers or hole excitations from the centers to the valence band. 
To g in the differential equation for m is added g vi + ga , and not g ci 
or g v2 , since g ci increases n as it decreases n, while g v2 increases p as it 
decreases p. The generation terms Ag m , Ag n and Ag p in the differential 
equations for Am, An and Ap are thust 

Ag m = Ag + Ag vi + Ag c2 , 

Ag n = Ag + Ag cl + Ag c2 , (58) 

Ag p = Ag + Ag vi + Ag v2 , 

and the generation terms that the equations for An and Ap now contain 
are respectively Ag m - Ag, t = Ag vi - Ag ci and Ag m — Ag p = Ag c2 - Ag v2 . 
The additional generation terms clearly represent the same processes 
as do the emission terms of (37). The distinction implicit in the notation 
is valid, however, consistent with zero values of these additional genera- 
tion terms at thermal equilibrium. Each generation rate of (37) is de- 
termined at equilibrium by the phonons and radiation associated with 
the equal corresponding capture rate. Since detailed balance applies also 
to the radiative part separately, there is no net radiation at equilibrium 
from any given process of capture and the corresponding generation. 

III. DETAILED THEORY AND APPLICATIONS 

3.1 Diffusion Length and Steady -State Lifetime Functions 

3.1.1 JAnear Theory 

The equations of (46) for two types of centers, when written for the 
steady state and linearized by neglect of the quadratic terms, give con- 

t The excitations involving trapping levels only, which may occur for large 
concentrations of centers (presumably with concomitant impurity-band conduc- 
tion), are here neglected. 
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centration increments that are proportional, and solving for An/ Am and 
Ap/Am provides the thermal-equilibrium trapping ratios. These and the 
corresponding lifetime t„, are thus given by 

_ P23P31 ~ P21 ^33 

* / 22' , 33 — ^23 ^32 
^21 ^32 — ^22 ^31 



r„ = 



^22^33 — V-2XV32 

k-1 



r m = — (vn + v l2 r n + v a r p ) (59) 

_ — (f22f33 — ^23^32) 

in terms of which, with the thermal-equilibrium diffusivity D ' from 
(31), the diffusion length can be expressed and the diffusion-length life- 
time evaluated. These results apply also for two-level centers if equilib- 
riums Vij are denned in accordance with (57). 

The case of single-level centers of one type lends itself to more detailed 
analysis. Results from the linearized equations for, say, acceptor centers 
only, for which r p is zero, are: 

Tn = —Vix/Vlt = j ; 

Vinl "P Vg n i + Vg p i 

_ T p0 (9li — n ) — T n0 n 

rpo(9li - n + n + ?ii) + T„ (po + Vi) 

= 9ll*(T <? ,i ~ T M ) 

(Oil* + tto)r,pi + PoTtnl 

= 1 - t„/t p , 
r» = (1 — r*)r m = — (j» 21 + j»22)/Ai 
= (vtpi + v pi + v a „i)/Ai 

T n o(n + Po + Pi) + Tpo(n + ni) 

9li* + n + po (60) 

_ (9ll* + Po)Ttnl +n T tp l 

9li* + n + po ' 

Tp = T„, = — ^22/Ai 

= (vm + f nl + V 0P l)/Ai 

_ t ,,0(9X1 - n + n Q + Mi) + T n0 (po + pi) 
9li* + no + po 

(9li* + n )r <p i + PoTtm 

9li* + n + po 
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The vn and release frequencies are equilibrium values; r n0 and t p0 , given 
by 

t„ = (C„i3li) -1 = (1 - no/9li)Ti»i = (1 + Pi/Po) _1 T(«i 

and (61) 

Tp0 = (C^ftd)' 1 = («. ) /9l 1 )r, pl = (1 + n,/no)"V to i, 

are the respective limiting lifetimes t in strongly extrinsic p- and n-type 
materials (in which they are also T m and t,,,,); A, , given by 

Ai = V n V2". — VnVft. = VtnlVgpl + VtnlVlpl T ''/yjl "*/»l 

= C nl C p i9li(3li* + no + po), 

is always positive if neither C nl nor C pl is zero; and 9li*, which will be 
referred to as the "capture concentration," is given variously by 

SJlj* = v tM r t pi/CniC p &li = niW(no + fh)(po + Pi) 

= noPtnl/Vgnl = PdVtpj/Vopl = ni(.TgnlTgpl/TtnlTtpl} 

= ou(9li - wo)/(l + «io) = Wo/(l + <*io) (63) 

= 9li(no/9li)(l - Tk/Vli) 

= Sliaio/Cl + «io) • 

The different forms for these results are obtained by use of (51) for the 
Vii , definitions of (50) for the capture and release frequencies, and 
equilibrium relationships of (38), (41), (42) and (53). The middle term 
of the second form for Ai is the one that gives rise to 9Ii*, and it follows 
from (53) that the first term is large or small compared with the third 
according to whether p is large or small compared with n . Capture 
concentration 9li* large compared with n + p is, as will be shown in 
Section 3.2.1, the condition that capture frequencies predominate over 
release frequencies. The volume rates of electron and hole transitions 
at equilibrium (see Section 2.2.1.2), n i> t ,a and jhPtja , are equal to 9Ii* 
times the corresponding release frequencies. The equations of (63) show 
that n and 9Zi - no may be written in terms of 01!* and a l0 . The concen- 
tration 9li* is small if the centers are nearly all ionized orun-ionized; the 
last form shows that its largest value is J3lx , which it assumes for a 10 = 1 
or w /9li = h that is, for the Fermi level coincident with the energy level 
of the centers. Entirely similar results, for which obvious notational 
changes are required in some forms, hold for donor centers only. 

t Conditions for these lifetimes are p » 9li* + p. + pi*, »o » 9li* + »i + »i*- 
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The diffusion length L rt and lifetime r corresponding to the electron 
and hole lifetimes of (GO), which are the Shoekley-Read lifetimes, may 
be evaluated from (31) or (35) and (60). These equations give 

I a) = Do't,,, = DuT = Dr,[l — r„po/(«o + Po)]t p 

= kTHnHplTo^inoTtpl + PoTtnl) (64) 

where a,,** and tr p0 are en„n Q and tyipPo . Other forms may be written by 
expressing T t „\ and t, p1 in terms of t„o and r p0 by use of (61). The diffu- 
sion-length lifetime for this case, 

TO = (Vgnl + V a p])/(v tn lVgpl + VtplVg„l) 

= [Tpo(no + »i) + r„o(po + Pi)]/(wo + 7>o) (65) 

= (n T t pi + poTi„i)/(n + Po), 

is formally similar to the familiar common lifetime 31,32 for both electrons 
and holes for the limiting case of 9li small, as inspection of (60) serves 
to verify. f Thus, L and t are, for given t„ and r p0 or T tl ,i and r tp \ , 
independent of 9li . For given energy level and capture coefficients, to 
is proportional to Dlj - . The true L and t apply, of course, in the linear 
part of the small-signal range, in which no appreciable trap saturation 
occurs. With small-signal trap saturation, diffusion length and lifetime 
that are usually considerably larger apply in the saturation range. 
These are evaluated in Section 3.1.2. 

It can be shown that the electron and hole lifetimes of (60) are sub- 
stantially equal to t» if 

I %tnl — V tp i | = | »Vl/"o — "ujil/po| 9*1* « (1 + «)(»«■] + »>„,») (66) 

holds, in which e is the smaller of (9Ii* + n )/po and (9l x * + p )/n , as 
given by the respective conditions \ t p — to \/t <5C 1 and | t„ — t |/t <C 
1. For extrinsic material, e may usually be neglected; including it pro- 
vides an appreciably weaker condition only if 9Zi* is larger than the 
equilibrium majority-carrier concentration. The condition of (66) may 
be severe: It is essentially 3!i* small compared with the equilibrium 
minority-carrier concentration for the minority- to majority -carrier re- 
lease frequency ratio of order unity or larger in extrinsic material, that 
is, for the minority-carrier trapping range defined in Section 2.2.1.2. 

t This formal similarity holds for any number M of kinds of centers, m being 
given by IS/'i (.nortpj + /-'<i7-/, u r l ] -l /("o + po), as may he shown from the first 
form for (ii„, of (71) and the observation that, if two or more different kinds of 
centers arc present, then (R,„ is the sum of similar terms for each kind. 
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General conditions for the validity of this linear analysis may be for- 
mulated as conditions for the neglect of the quadratic terms. For this 
purpose, assume uniform concentrations and volume-generation rate. 
Then (46) and (51) yield 

dAn/dt = Ag - v tal An + v gnl An + CmAnAn = 

and (67) 

dAp/dt = Ag - vt P iAp - v upl An - C pl ApAn = 

for acceptor centers only in the steady state and no direct recombination. 
The conditions may be derived in a self-consistent manner by first ob- 
taining, with the neutrality condition, the concentrations from the 
linearized forms of (67). These concentrations, namely 

An = (v t pi + vgpi + v g „i)Ag/Ai , 

An = (v M - Ptpi)Ag/Ai, (68) 

Ap = (Vtnl + Vgnl + v ap i)Ag/A x , 

are then substituted in (67), so that conditions for negligible quadratic 
terms may be obtained as restrictions on (positive) Ag and, by use of 
(68), as corresponding restrictions on the concentrations. It will suffice 
to give the former restrictions, which for the neglect of C ia AnAn and 
CpiApAn are, respectively, 

Ag/Ai « V t nl/Cnl(Ptnl ~ Ptpl) = (9ll ~ fo)/(Vtnl ~ V tp i) 

and (69) 

Ag/Ai « Ai/C P i(p t ni — Ptpi)(Ptni + Pom + "»pi) 

for v t „i > p tp i . If p tp i > v t ni , the restriction for neglect of one of the 
quadratic terms turns out to be that of (69) for the other one, but with 
subscripts n and p interchanged — an interchange that does not affect 
Ai . This distinction arises because the signs of the approximate linear 
terms in An depend on the sign of p t „i — Pt P i ■ It is easily shown thatf 
for this quantity zero, or the case of "critical recombination," An is 
identically zero and (67) are linear for all Ag. For trapping without 
recombination, (68) does not apply because A x is zero, but the conditions 
may properly be written as the restrictions on the concentrations ob- 
tained by use of these equations. For example, for electron trapping 
with Cpi zero, the condition An « n + n { results, which may be a se- 
vere condition for p-type material. 

t A solution is excluded that does not admit thermal equilibrium, for which 
Ag or Mi is zero for certain negative values of An and Ap or Ap and Ag, respectively. 
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3.1.2 Nonlinear Theory 

For added carrier concentrations resulting from arbitrary injection 
levels, steady-state lifetime functions r„ and t p may be evaluated from 
(R„, . For acceptor centers only and no direct recombination, (R m is given 
by 



(70) 



(R m = Cni[n(9li - n) - mil) 
= C pl [pn — pi(9li - n)], 
which results in 

(R. = VP ~ "* 

T P a{n + Jli) + T„o(p + Pi) 

w(9li — n) + (wo + n t )n + (go + pi)(9l x - n) An 

T P o(9li — n) — T ,Ji " m,i (71) 

= w(l — n/9li) + (n + n x — p — pi)An/9Zi + n + p „ 
T p o(9li — n) — r»oW 

= An/v,, = Ap/r p . 

Eliminating n by means of the second equation of (70) and the use of (61 ) 
results in the first formf for 6l m of (71). This familiar form 31 furnishes 
r„ or t p in terms of An or Ap alone if one of these concentrations is elimi- 
nated by solving the second equation of (70) written with n replaced by 
p — n; and An or Ap may at the same time be related to, say, the 
generation rate Ag = (R„, for steady-state photoconductivity. The algebra 
involves radicals. A better procedure for such analysis employs the 
second or the third form for fll„, ; these result from (70) by elimination 
of n and p with the neutrality condition. Then An and Ap are, with (R m , 
written in term of n or An as independent parameter in accordance with % 

• _ Tpoin-o + nQ + T„ (n + Po + Pi) . - 
Tpo(9li - n) - T„ n 

and (72) 

Ap = Am = T »o(Po + Pi) + T p0 (9Ii - n + n + wi) A? . 
r„o(3li — n) — T„ n 

so that the lifetime functions are given by 

t If two or more different kinds of centers are present, then (R m is the sum of 
similar terms for each kind. For the corresponding lifetimes in terms of An and 
Ap, see Okada. 15 

J Note that An has the sign of the denominator, which is proportional to v tn i — 
v lp \ for Ail small. 



-1 
r„ = 
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TpoUo + tti) + T„o(w + Po + Pi) 

T " " n(l - nfSlx) + (n + n x - Po - Pi)An/9li + n + Po 
and ( 73 ) 

_ T n0 (p + pi) + Tpo(9li - n + n + n x ) 

Tp ~ n(l - n/fJli) + (n + n x - Po - Pi)Art/'Sl l + n o + Po' 

The functions reduce to the Shockley-Read lifetimes of (60) for the 
equilibrium value n of n, as may be verified by use of (41) and (42). 
The range of n is from n to the limiting large-signal value given by 
n/ffli = t p o/(t„ + r p o), for which the denominator in (72) vanishes, 
and for which t„ and t p both equal r„ + t p0 . 

The trapping ratio r„ corresponding to the lifetime functions of (73) 
is given by 

Tpo(3tx - n + n + nQ + T n0 (p + p0 
r P o(9li — n) - Tnon 

' , TpoAn + r n0 Ap "1 Ap 

. + r p o(3li - n + n + n,) + r n0 (po + Pi) J An (74) 

T p o 2 (n + ni) + r, l0 2 (po + Pi) 

_ Ap + rpor„ (9I 1 + n + tti + p + Pi) A? * 

"An [r p0 (9ti — n) — r n on]- 

which is obtained from (72) by differentiating with respect to Am. The 
equilibrium value of r„ , which is that of An/Ap, is the r„ given in (60), 
while the limiting large-signal value is zero, as may be expected. By 
means of (72) and (74), the steady-state continuity equation, (30), may 
be written (for acceptor centers only), with An as independent variable 
and the components of grad An as dependent variables. The second or 
third form for r~ l of (74), with Ap/An given by (72), lends itself to 
this purpose; note that dAn/dAn and dAp/dAn equal r n i - 1 and r n , 
respectively. 

The lifetime function 31 for | An J « An ~ Ap - Am, 

T P o(n + m + Ap) + T n0 (po + Pi + Ap) , 7 s 

T " Tp no+po+Ap 

may be derived most directly from the first form for (R m of (71). By solv- 
ing the second equation of (70) for An, the condition | An | « An 
~ Ap may be written as 

. I Vgnl + Vgy\ ^ | Vtnl ~ Vtpl \ = \ Pgnl/WQ ~ Pgpl/Po \ ^ * ( 7 q\ 

Cnl ~\~ C p i C n l "T C p i Cnl T ^pl 
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Equilibrium release probabilities are here employed. Equation (76) with 
Ap set equal to zero is the condition that An be relatively small for all 
Ap. Since this condition is (66) with e set equal to zero, it subsumes the 
condition for equilibrium lifetimes substantially t , to which it is usually 
equivalent. 

It is readily shown from (75) that, withf Ap « n + Po (and equilib- 
rium lifetimes to), the lifetimes are substantially t if Ap is small com- 
pared with (vy„i + Vg„i)/{C,a + C p i). This condition and the one of 
(66) may be severe conditions under essentially the same circumstances. 
That is, in the minority -carrier trapping range defined in Section 2.2.1.2, 
lifetimes are t for 9li* small compared with minority-carrier concentra- 
tion n or po ; and then, consistent also with the condition of Section 
3.1.1 for the neglect of C^AnAn or C pl ApAn suitably specialized, for Ap 
small compared with n + ih\ or p + pi . If the condition on Ap is not 
met, then, with the condition on "SI*, (75) gives a lifetime that increases 
rapidly with injection level at low injection levels. J But such observed 
behavior with extrinsic material, as these considerations indicate, can- 
not usually be properly analyzed by use of (75). The steady-state life- 
times in the small-signal range generally either result primarily from 
recombination or majority-carrier trapping and are both to and sub- 
stantially constant, or else have distinct equilibrium values given by (60) 
with dependences on (small-signal) injection level obtainable by 
the general procedure described. It will be shown that, in the latter case, 
substantially constant apparent diffusion -length lifetimes given by T gni 
for n-type material or t„„i for p-type generally apply in the small-signal 
range above a certain injection level. Thus, unless trap concentration 
is quite small, (75) has significant application in the former case only 
to the transition from t to the lifetime t„ + t,, for the large-signal 
range. 

Nonconstant small-signal lifetime functions are associated with deep 
traps in the minority-carrier trapping range. Such traps will be saturated 
(in the steady state) even in the presence of a concentration of mobile 
minority carriers that is relatively quite small. From (67), by equating 
dAn/dt and dAp/dt (which are zero also in the immediate context), An 
may be written as 

An = (^ nl An - v tpl Ap)/{v unl + v upl + C nl An + C ;)1 Ap), (77) 



t The more general condition without this restriction includes to ■*«• r»o -f- r p0 
for huge values of Ap. 

X As a result of saturation of centers available for minority-carrier capture, 
this lifetime increases essentially linearly in the small-signal range from the 
equilibrium value t p o(«o + «i)/po or T, l0 (po + P\)/n Q and asymptotically to the 
large-signal value t,,o or r n o . 
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in which the concentration-dependent contributions to the release fre- 
quencies are exhibited separately. Suppose that the transitions in p-type 
material are electron-dominated, so that v lini » v upl . For deep traps, as 
defined in Section 2.2.1.2, wi « p (or p, » no) holds, which gives 
n + fti « po + Pi and therefore implies C«i >» C pl or r p0 »> r„ for the 
present case. Then, for An » n -f- »i , the denominator in (77) is 
C n iAn. Furthermore, v tn i = C nl (9li - no) >» v tP i = C pl ri holds from 
(53). Hence Aw ~ 3li — n follows f for An large compared with n + Wi 
and not too small compared with Ap. If trap concentration is not too 
large, small-signal saturation evidently occurs under the conditions as- 
sumed. If it is large, then a large conductivity increase is associated with 
the majority-carrier concentration corresponding to the saturated traps. 
Nonconstant small-signal lifetime functions apply in either case, whether 
saturation occurs in the small-signal range or not. 

The lifetime functions in the saturation range approach the limiting 
large-signal lifetime, r n0 + r p0 , substantially equal to r p0 . Though r p0 
is otherwise the minority-hole-capture-limited hole lifetime in strongly 
extrinsic n-type material, in this case it is a lifetime limited by majority- 
carrier capture. For small-signal saturation, An changes relatively 
slightly from the small-signal saturation range to its limiting value, 
(v tnl - v m )/{C,a + C p i) = 3li - no - r (1 o9li/(r„o + T p0 ),forthe large- 
signal range. This circumstance might suggest that t„ applies over both 
ranges. In general, it does not: The denominators in (73) are compara- 
tively small (reducing, for example, to n + n x for An = 2Hi — no) and 
are sensitive to very small changes in An; a very small change in con- 
centration of unsaturated traps can affect lifetimes appreciably. As will 
be shown in Section 3.2.2, large-signal lifetime implies relatively large 
increase in conductivity. The equations of (07) for the steady state, 
simplified for relatively small departure of An from 3li - n still (neces- 
sarily) nonlinear, may be solved, in terms of An and for the saturation 
range, for the lifetime functions r„ = An/Ap and t p = Ap/Ag. With 
v g „i » v 0P i + t p o~ 1 , which will still apply in the present case even if Tfii 
is of order p„ + pi , it is found that t„ ~ r p0 /[\ + (p., + 9li - n )/An] 
holds for the saturation range, specified by An » n + ih . Thus, 

t„ ~ r p0 An/(po + 9li — n a ), 

proportional to An, holds for the saturation range of relatively small An. 
An apparent diffusion -length lifetime, t ', may be found by evaluating 

t The equivalent condition Ap » 9li - A + n + wi from (72) takes into ac- 
count Ap > An. The small-signal saturation value of TC, - n may be appreciably 
larger than n + »i , but its limiting large-signal value is small compared with 
no + n, if 9Ii « («v,iA B ,>!)(po + Pi) holds. 
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(D'/D)t p = (1 — r n )T p for small-signal saturation.! With r„ _1 = 1 + 
dAn/dAn, the expression 1 + (w + ni)(3li — n )/An is found for 
1 — r n ; and dividing Ap ~ An + 9l x — n by Ag- gives 

t p ~ [1 + (3li - />o)/An]r„ 

~ t p o[1 + (9li - ri )/An]/[l + (p + 3li - no) /An]. 

It is easily seen that, for small 9li — w of order n + Wj or less, the life- 
time function for to' so obtained is r p0 Att/po , as are t„ and t^ . As may 
be expected, this result is consistent with (75). For 9li — n G >» n + «i , 
however, the lifetime function gives t</ ~ t p0 /(1 + t p0 /t (/p1 ) for 
(9li — no) 2 » An 2 » (n + ni)(9li — no). The condition 9li « p + pi 
then gives r u ' ~ r„ P i for small-signal saturation. If this inequality is 
reversed, then t ' ~ r p0 results, and saturation occurs with relatively 
large increase in hole concentration. 

From (60), the equilibrium electron, hole and diffusion-length life- 
times are, for these cases, generally small compared with r pn . They are 
given by t ~ r» ~ [(no + Wi)/jl ]t p0 = [(no + ni)/(9li - n )]r tfpl and 
t p ~ [(3li — no)/p Q ]Tpo = t op i if small-signal saturation occurs, for 
which ro is also small compared with t ' = r gp \ . The minority -carrier 
and apparent diffusion-length lifetime functions increase with injection 
level, most rapidly as An becomes comparable with n a + ni and the 
traps fill. These results clearly provide a simple model, based on a single 
trapping level, for superlinearity, 29 ' 34 ' 35 ' 36 the more-rapid-than-linear in- 
crease of photoconductivity with injection level. J With small-signal satu- 
ration, two superlinear ranges may occur, the first as diffusion-length 
lifetime increases from r to T tfp] , and the second as it increases from 
T fl pi to T p0 in the large-signal range. With large-signal saturation result- 
ing from large concentration of traps, only one superlinear range occurs, 
since a nearly linear intermediate range is absent. Only one range oc- 
curs also under the condition of (66) for small trap concentration. It is 
evident, however, that with superlinearity this condition is generally 
quite severe. 

For the majority-carrier-dominated case of v tp i » v ln \ (or v (/pl >» v gn i) 
in p-type material, there can be no small-signal saturation. With small 
trap concentration, lifetime r ~ (1 + p\/pa)r„o = Tt„\ , which is limited 

t The second form follows since An/ (no + »i) » (9li — n )/(po + 3li — ra<>) 
holds a fortiori. 

X For small concentration of centers, Ap may exhibit a less-rapid-than-linear, 
a linear, or a superlinear dependence on Ar/, as Rittner 1 has shown using a lifetime 
function tantamount to that of (75). From this equation, superlinearity results, 
as may lie expected, if t m0 + r p0 exceeds m , so that the numerator increases more 
rapidly with Ap than the denominator. See also Ridout, 70 Newman, Woodbury 
and Tyler" and Sandiford. 7 
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by minority-carrier capture and obtains over the entire small-signal 
range, then changes in accordance with (75) to t„o + r p o in the large- 
signal range. This change is a decrease to r„ n » r p0 if the trapping level 
is near the Fermi level or higher. 

The steady-state fractions of ionized centers can be represented by 
simple formal generalizations of the equilibrium relationships of (41) 
and (56). In these equations, n and pa are replaced by n and p and «io 
and a 2 o by 



^ Cm n + C p ipi 1 + An/ (no + nj*) 

m ~ C pl p + C nini ~ ai0 1 + Ap/(po + Pl *) ' 

_ C p2 p + C„ 2 n 2 _ 1 + Ap/( p + p-i*) 

C„ 2 n + C p2 p 2 1 + An/(n + n 2 *) 



(78) 



as can readily be shown f by solving for the ionized fractions from (37) 
and also from the corresponding two-level equations of Section 2.2.2. 

3.2 Photoconductivity 

A number of results for steady-state photoconductivity being implicit 
in Section 3.1, the present section will deal principally with the transient 
decay. 

3.2.1 Linear Theory 

For two types of centers in the linear small-signal case the time de- 
rivatives of Am, An and Ap for photoconductive decay are given respec- 
tively by (46) without Ag and the quadratic terms. The general solution 
is accordingly 



(79) 



in which the Aj are constants determined by the initial conditions, and 
the r nj and r p] are trapping ratios for the respective decay modes deter- 
mined by 

t The equation given in the abstract of the paper of Sah and Shockley 22 re- 
written in the present notation yields rt/(9l — -A — $) = ai and (91 — il — p)/fi = 
or 2 -1 , from which the ionized fractions for the two-level case here given follow as 
solutions of simultaneous linear equations. 



Am 


= T,A ie -» l , 

7=1 


An 


= Jlr nj A j e->! 1 , 

;=1 


Ap 


3 
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(v u + fj) + vnr nj + vnTpj = 0, 

vn + ("22 + r/)r Bi + p 2i r pj = 0, (80) 

"31 + V32T „j + ("33 + "i)^pi = 0, 

with the decay constants v } - being the roots of the equation obtained by 
equating to zero the determinant of (80). The Aj are found in terms of 
the trapping ratios and the initial concentrations Ami , Ani and Api by 
setting t equal to zero in (79) and solving. The solution so obtained ap- 
plies as well to the two-level case, for which the va are defined in accord- 
ance with (57). The decay constants Vj , which are roots of a cubic, are 
always real and (since the coefficients alternate in sign) positive, except 
that one of them may be zero. Establishing these properties involves ex- 
pressing the coefficients in terms of the capture and release frequencies 
by means of (51) or (57) and making use of (53). 

The constant term of the cubic and one decay constant are zero if there 
is trapping only and no recombination. This case can occur essentially 
in two ways: The two types of center may trap, respectively, the two 
kinds of carriers, or they may both trap only one kind. If, say, the ac- 
ceptor centers trap only electrons and the donor centers trap only holes, 
then the vj are readily found to be zero, v tltl + v an \ and v tP 2 + v 0P z . 
The last two decay constants characterize the respective exponential in- 
creases with time of An and Ap to new equilibrium concentrations after 
injection that correspond to the zero decay constant. These equilibrium 
concentrations are fractions of Am (which remains constant) equal to the 
fractions of the time the electrons and the holes are trapped. Indeed, 
since the two types of centers trap independently in this case, the solu- 
tion consists of solutions written independently for each. But if, for ex- 
ample, electrons only are trapped by both types of centers, then this in- 
dependence does not obtain ; electrons released from centers of one type 
may be trapped by centers of the other type. With the convention here 
employed, concentration of electrons trapped by donor centers may then 
be written as negative Ap. The decay constants are found to be zero and 

\{v ln \ + Vgnl + V t „2 + V a „2 ± [(»«„i + V gnl — V tn 2 ~ Pgni) 2 + ±V M V lni ] }, 

with an equilibrium value after injection of total trapped electron concen- 
tration equal to Ap/[1 + v unl v gn 2/(v t „iv on -> + v^iw)]. 

The general linear small-signal case for one type of center is readily 
evaluated in detail. For acceptor centers only, the solution is given by 
the first two equations of (79), all terms withy = 3 being omitted. The 
trapping ratios are given by 

r»j = — "2l/("-'2 + Vj) = (V M — Vtpl)/(v t nl + JW + "„pl — Vj) 

= -("n + Pj)/ph = {vj - vtp\)/v ap i, (81) 
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and the vj are the roots of 

vf - v. Vi + Ai = 0, j = 1, 2, (82) 

where v s is denned by 

V s = — (fll + Von) = Vtnl + Vgnl + Ptpl + "opl (83) 

and A x by (62). The decay constants are thus 

ry-MC-l^V+id, J -1,2, (84) 

with 

,/ r = (j/ - 4Al)* = [(^nl + Vgnl - Vtpl - Pgpl) + 4:V g7 ,lV gp l]' . (85) 

The corresponding time constants n = ^i -1 and t 2 = »»2~ are also equal 
respectively to i» 2 /Ai and ^i/Ai . Nonoscillatory decay is easily verified 
for this case: The second form for v T shows that the Vj are real; and, since 
v r < v 8 , the vj are positive. 

A subcase that provides some physical interpretations is that of 9li 
sufficiently small so that capture frequencies are small compared with 
release frequencies. As (62) and (83) show, the condition j// » 4A X then 
holds, and expansion of the radical in (85) gives 

Ti ~ P s ~ = T„nlT,,pl/ (t„„i + r (/p i) 
«To ~ V>/A\ ~ To . (86) 

Thus, for this subcase, t 2 is the steady-state lifetime r of (65). It is 
large compared with n , the time constant for the adjustment of An to a 
fixed fraction of Ap substantially equal to the equilibrium trapping ratio, 
?•„ . This interpretation of n follows from solutions for the concentra- 
tions: The last form for the r nj of (81) and r„ from (60) give 

r„l ~1 + Vgnl/ Pgpl , 

T„2 ~ {Vtnl — Vipl)/{v„ n \ + V api ) ~ V n J (87) 

and the result 

Aw/Api = [r n Vg„i/(v a „i + v 0P i)}e~ n 

+ [1 - r n Vgnl/(Vgnl + foplW'' 2 , 

Aw/A Pl = r n (-e-" Ti + <T ,/T2 ), (88) 

Ap/Api = — [r„iVi»i/(*W + v aP i)]e~ Tl 

+ [1 + r n v apl /(v unl + v opl )]e~ 11 
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holds for initial concentration Awi zero. Since r n is small, the mobile- 
carrier concentrations are mostly in the second or lifetime decay mode 
and differ only slightly. For this subcase, if the initial trapping ratio is 
r n rather than zero, then the first decay modes are not present and An is 
r n Ap = r n Apie~" Ti . The first decay modes do not occur either for "criti- 
cal recombination," with which An remains identically zero as a result 
of equal capture frequencies v ln i and v tp i or, for this subcase, equal cap- 
ture rates for An and Ap. For small 9Zi , the capture rates are in all cases 
substantially in the ratio vtni/vtpi • In linear cases, they also decay in 
the lifetime mode after this mode predominates. The release rates be- 
have similarly, their ratio being equal to v gn \/v gp i , or to (no/po) (v tn i/v tp i) 
in accordance with (53). 

The condition for neglect of the capture frequencies may be severe. 
The approximate form of v a applies if v tni + vt P \ « v gn \ + v gp i holds, 
which implies that 

9li* « (»W + v gpl )/(v an i/nQ + p gp i/p<i) 

= (n v t „i + poVt P i)/(vtni + v tp i), 
a condition which subsumes 

91,* « n + p (90) 

for neglect of vtnivtpi in Ai . The conditions of (89) and that for steady- 
state lifetimes equal to r of (66) are the same for the minority-carrier 
trapping range defined in Section 2.2.1.2, for which they are 9li* small 
compared with the equilibrium minority-carrier concentration. The con- 
dition v„ » 4Ai is 

9Ii* « in, 2 (w + Po)~ l {v unl / v„ pl + v giA /v unl + 2) 

i 2/ , s-i/ , n2, v^ 1 -* 

= jUi (Wo + Po) {Tgnl + T gpi ) /TgnlTgpl 

if (89) holds, and it can be shown to be weaker than (89) in general if 
the minority- to majority-carrier release frequency ratio exceeds a num- 
ber that is about three for extrinsic material and about six for intrinsic 
material. f Equations (89) and (91) are both subsumed a fortiori by 
9li* « n?/(n + po), which is 9Xi « (n + n0(p + Pi)/(n + p ). 

The release frequencies may be neglected under the condition of (90) 
but with the inequality signs reversed. The solution is then simplyt 
An/Ani = e~" T,n} and Ap/Ap t = e~" Tlpl . For 9li* large, (60) shows that 

t Equation (91) gives a stronger or weaker condition according to whether 
nuf ai ,i/pi, n i + pop i ,ni/i' ll pi is smaller or larger than 3(n + po). 

t This result easily follows directly from the differential equations. Or, note 
that the radicand in (85) is (vt„i — vt P i) 2 . 
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Tt„i and T tp i are respectively the steady-state lifetimes t„ and r p . The 
condition v? » 4A X is accordingly \{tJt p 4- t p /t„) 4- \ » 1 — namely, 
that one of r„ or t p be small compared with the other. If r„ or t„ is the 
smaller, then substantially all of An or Ap, respectively, is transformed 
comparatively rapidly into positive or negative An, after which a slower 
recombinative decay of Aw and the concentration of the other mobile 
carriers takes place as these carriers are captured. 

The condition v 2 » 4Ai implies n « t 2 , with n essentially a charac- 
teristic time for trapping and r 2 essentially a lifetime. This interpreta- 
tion does not apply if v 2 and 4Ai are comparable so that ti and r 2 do not 
differ by much. For small 9li and the majority-carrier trapping range, 
for example, n ~ t 2 may hold; (89) may apply, but not (91) (see foot- 
note on previous page). The case of v 2 ~ 4A X for 3l x large, for which 
t\ , t-> , Ttni , T tp \ , r n and t p axe all substantially equal, is a case of re- 
combination with but slight trapping. 

The general trapping time and lifetime obtained from (84) and related 
equations are 

ti = v' 1 = r onl T flP i/[(9li7Po + 1)t„„i + (9li*/«o + lKpl] 

= r lnlTlpl /[(l + po/3li*)r ( „i + (14- n /9li*)r (p] ] 

«r 2 = y,/Ai = (Oli* + n 4- p ) -1 

'[n Tg„i 4- p T op i + n?(T gn i 4" r„i)/9li*] 

= (3li* + n 4- Po) _1 

■[(9li* + Po)r tBl + (9li* 4- no)T„i]. 

Comparison with (60) and (65) shows that this lifetime t 2 is larger than 
the steady-state lifetimes r„ , r p and to ; all are equal in the limit of 9l x 
small. For 9li large in intrinsic material, t 2 equals 2r . Furthermore, 
these lifetimes all decrease monotonically to zero as 9li increases in- 
definitely. 

The decrease of r 2 with increasing 91! may, however, proceed essen- 
tially in two ranges, with approximate constancy of t 2 in an intermediate 
range. 41 From the first form for t 2 of (92), this intermediate range occurs 
provided there are capture concentrations 9li* that are small compared 
with n + po and also large compared with 

(Vg„l + V gpl ) I '(Pgnl/rk + Vgpl/po) ] 

that is, if the strong inequality 

• T gnl 4- T 0P i <« (wo/po)t„»i 4- (po/n )T gp i (93) 
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holds. It can hold for sufficiently strongly extrinsic material if the ma- 
jority-carrier release time is not too small. For small 9li , t 2 varies in- 
versely with 3li , as (65) for t shows. For large 3li such that 3l x * » 
«o + Pa , to varies similarly, equalling the value (n T g „i + Pot op i)/ 
(w + po) of approximate constancy divided by 9li*/( w o + Po)- With the 
third or fourth form for 9l t * of (63), this n reduces to r tn i + T tp i . Since 
ti for large 3li* is the harmonic mean of T tn \ and r (p i , ti is the smaller 
of these capture times and t 2 the larger, as previously discussed for this 
case. It can be shown that, for the minority-carrier trapping range, the 
inequalities that 9li* must satisfy for approximate constancy of n gen- 
erally imply the condition v t » 4Ai on which the calculation is based. \ 
A similar situation has been shown to obtain with the inequality for the 
case of negligible capture frequencies. But since this case involves a con- 
dition for neglect of the capture frequencies that is usually severe in 
the minority-carrier trapping range, it is the present case that would 
usually apply in practice in this range. 

3.2.2 Nonlinear Theory 

Although the general problem of photoconductive decay is intractable 
analytically, some special cases can be solved and certain techniques of 
approximation are effective. From (46) and (51 ), general equations that 
apply for centers of the acceptor type may be written as 

dAn/dt = Ag + v gnl An - [v tlll - (C tll - d)An\An - dAn, 

dAn/dt = [pmi - v tpl - (C nl + C pl )An]An 

— {vtpi + v B ni + v„i)An - C P iAn (94) 
= [vtni ~ vtpi - {Cni + C p i)An]Ap 

- (vmi + pgni + v 0P i)An + C„iAn 2 , 

dAp/dt = Ag — v gpl An - [v tpl + (C pl - d)An]Ap — dAp 2 , 

in which equilibrium values of release frequencies are employed. Since 
| Aw | is bounded (by a concentration that cannot exceed 9li), it is clear 
that, if the initial concentration Ap x is sufficiently large, then the decay 
proceeds with An after a short transient substantially equal to 

(v M - v tp i)/(da + C pl ), 

t One inequality is the reverse of that of (89); hence it is the condition for the 
neglect of the release frequencies in v. . The other inequality is (90). It follows 
that k„ 2 » 4d| is 91,* » 4(/i _1 + pC^VvnWapil '("am/no + v gP i/po) 2 . This condi- 
tion is weaker than the reverse of (89) under the same circumstances that make 
(91) a weaker condition than (89). 
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for which the coefficient of An or Ap in dAn/dt is zero. This value of An 
corresponds to maximum (but not necessarily complete) saturation or 
neutralization in the traps. Although direct recombination, characterized 
by the quasi-hyperbolic decay law 

Ap/Api = (CtApit + l) _i (95) 

predominates in principle for very large Api , it can frequently be neg- 
lected. 48,49,50 The large-signal decay is then exponential with lifetime 
(Cni + C p i)/(Cniv tp i + C pl p lni ) = t„o + r p0 , the limiting large-signal 
steady-state lifetime. With the limiting value of An, this result follows 
from, say, the last equation of (94), in which neglect of the release rate 
j> 07 ,iA(31i — n) is consistent with Ap large. 

Examination in further detail of the large- and small-signal decay is 
facilitated by the equation for negligible direct recombination, 

TnoArC 1 dAn/dt + T p0 Ap~ l dAp/dt 

+ 1 + [(po + pi)(l - An/Ap) + (n + n0(l - Ap/An)]/3li 

= KoAn- 1 + r.oAp- 1 ^, (96) 

which is readily obtained from (67) or (94) as a linear combination of 
dAn/dt = d(Ap — An)/dt and dAp/dt that eliminates the quadratic 
terms. For example, consistent with results for the linear case of small 
9li , assuming in this equation the steady-state trapping ratio r„ of (87) 
and (so that An and Ap are proportional) a single decay time, this time 
is given as the lifetime t 2 ~ t of (86). If either An ~ Ap or 9li is suf- 
ficiently large, then the term with brackets, which arises from the terms 
involving release frequencies, may evidently be neglected. It is other- 
wise plausible that release does not appreciably affect large mobile car- 
rier concentrations, while capture predominates with large trap concen- 
tration. For no volume generation, (96) may then be integrated, Avith the 
result 

(AnAp) i (An/Ap)- i<c » l - Cpl)/(c " l + c ' ,l) = A e ^-"^ + ^. (g?) 

in which A is a constant determined by the initial concentrations. It is 
easily verified that, besides furnishing the large-signal lifetime, (97) is 
consistent with the linear solution for large 9li , for which the release 
frequencies may be neglected and the decay times are T ln \ and r lp \ . 
This equation is a first integral of (94) for a case of large 9li , one which 
can accordingly be formulated as a first-order (rather than second-order) 
nonlinear differential equation. 

The condition under which (97) holds is, from (96), 
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91, » | Aw[(«o + ni)/An - (p + Pi)/Ap] I . (98) 

From (96) and (98), the large-signal lifetime t, i0 + t p o obtains if 
An ~ Ap » A;i holds and Ap » 3li _1 | An(n + »i — po — Pi) I also. 
With the An for Aw ~ Ap from (77), these conditions are respectively 
that of (76) and 

Ap + (V 0ni + V„pl)/(Cnl + C„l) 

- Ap + to(-«o+ Po)/(r„o + t^) 

» | Wo + Po — ("„nl + V up i)/(C„i + C p i) I 

= («o + po) I 1 — to/(t„o + r p o) | • 

Since the left-hand sides of (76) and (99) are the same, comparison of 
the right-hand sides will indicate which condition is the more restrictive 
in any particular case.f The condition that corresponds to (99) obtained 
from the lifetime for An ~ Ap of (75) is similar except that the constant 
term on the left is replaced by Wo + po . Setting Ap equal to zero in either 
gives the condition that t„ + t„ apply for all Ap, which is that it 
equal to . 

The decay times associated with a small-amplitude pulse of added 
carriers above a steady generation level A</ are readily evaluated. The 
equations for dbn/dt and dfrp/dt, linear in the concentration increments 
bn and 5p that result from the pulse, may be obtained from (67) . Written 
with capture and release frequencies that are concentration-dependent, 
they are formally the same as the linear small-signal ones for dAn/dt 
and dAp/dt. For the release frequencies, the definitions of (50) apply; 
for the capture frequencies, n in these definitions is replaced by ii. The 
condition v s 2 » 4A, of Section 3.2.1 generalized in this way is the condi- 
tion for a lifetime f-> for fin and 8p equal to the generalized ratio i»,/Ai 
and large compared with the corresponding time constant for trapping. 
The lifetime f 2 depends on the steady-state values of Aw, An and Ap; 
(72) gives An and Ap in terms of Aw, and (71) relates An to Ag = (R„, . 
It reduces to t 2 of (92) for the linear small-signal case and to r„o + Tpo 
for Ag large. 

The approximation An ~ Ap applied to (94) gives the differential 
equation 

**t* + " + (»° +_p;)"'f A f_ = o, (ioo) 

[1 + (r n0 + r p0 )r Kwo + Po) 'Apjro 

t In t lie minority-carrier trapping range, for example, (70) and (99) may re- 
quire approximately that Ap be large compared with trap concentration and with 
equilibrium majority-carrier concentration, respectively. 
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whose solution 1 ' 44 may be written as 

Ap p + (», + po)-^ ] 1 '""- 1 ^^'-' = - W r > (10] ) 

Api |_1 + (wo + po) _1 ApJ 

As may be expected, this solution is the integrated form that corresponds 
to the lifetime of (75) ; the approximation results in the steady-state An 
of (77), and the condition required is that of (76). Casual inspection of 
( 101) might suggest that t and t b0 + r p0 apply respectively for Ap small 
and large compared with n + po ■ This conclusion is, of course, illusory: 
For the minority-carrier trapping range, the exponent in (101) is large 
and to applies only if Ap is restricted as explained in Section 3.1.2. 

The limiting large-signal An, Ap and lifetime for two kinds of centers 
are easily evaluated from (46) and (51), and the extension to any num- 
ber of kinds of centers is obvious. As may be expected, the values 
(vtni - v tv x)/(Cn\ + C p i) or (v [p2 - v tn2 )/(C n2 + C p2 ) of An or Ap are 
as if the acceptor or donor centers alone were present; and the lifetime 
is the harmonic mean of lifetimes t„o + t p0 for each kind of center, the 
decay constant being the sum of the separate decay constants. This result 
does not apply to the two-level case: From (57), An and Ap are found to 
equal C nl C^ffl/(C n iC n 2 + C n? C pl 4- C pl C p2 ) - n and C pl C p2 Vl/(C nl C ni 
+ C n2 C pl + C pl C p2 ) - po , with (1 + C nl /C pl + C p2 /C n2 )/(C nl + C p2 )9l 
as the large-signal lifetime. 

General solutions for trapping only and no recombination can be ob- 
tained without difficulty. For, say, electron trapping by acceptor centers, 
Ap maintains its initial value Api , and the nonlinear equation for An 
that results from replacing An by Api — A?i in the first equation of (94) 
has the solution 

. Vt 4~ Vtn\ — Vynl ~ ClApi -y,t 

An = An 2 = V ' ± V '» l ~ > " 1 ~ C " l * Vl (102) 

1 __ — Vt + Vtnl + Vgnl + UllAgl -» t l 
Vt + V tn i 4" Vgnl + CnlApi 

for initial value Api , if direct recombination is neglected, f with 

Vt = [(vtnl - Vgnl - C„iApi) 2 + ^PhOWmaft 

An 2 = \Cm~\vt - V M - V un l + C rl l&Pl). 

The concentration An 2 is the new equilibrium concentration which An 

t The general form of this solution is not changed if direct recombination is 
included. 



(104) 
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approaches asymptotically after injection. For Api large, these equations 
give 

An = [Ap, - (l - '^f) (9li - *) ] 

.M _i_ 9^i — no c -c Hl * Pl t\ 

Thus, An = Api — An rapidly increases to the limiting value 
[1 - (no + w,)/Api](3li - no) ~ 3li - no = v tn i/C n i , 

which corresponds to substantially all traps charged. For Api small, the 
equations give 

An/Ap, = h„i/(".»i + »w)] [1 + (*«/***) c- ( """ +y " nl> '] , (105) 

which may be obtained also by suitable specialization of results for the 
general linear small-signal case. According to (105), An/Api decreases 
from unity to T tn \/(T tn i + T l/nl ), the fraction of the time electrons are 
free, while An/Api increases to T gn \/{Ttn\ + r„„i), the fraction of the time 
electrons are trapped. An effect of slight recombination on An would in 
all cases be a comparatively slow decay from a value approximately 
equal to the equilibrium value A/i 2 for trapping only. 

It is sometimes relevant to deal with a model involving centers that 
provide nonrecombinative trapping in conjunction with other centers, 
of a suitably idealized type, that provide only recombination that can be 
specified simply in terms of a constant lifetime. Such centers would in 
general be present in comparatively small concentration, so that the 
amplitude of their trapping transient is negligible. Furthermore, this 
transient would be comparatively brief, so that steady-state lifetime 
applies after negligible time. 

With certain restrictions, the idealized centers may be centers that 
function in the recombination range or in the majority-carrier trapping 
range, f The An or Ap for these centers obtained by setting dAn/dt or 
dAp/dt equal to zero results in a contribution to both dAn/dt and dAp/dt 
that is the negative of a steady-state recombination rate similar to that 
of the first form for (R„, of (7 1 ) . With subscripts ' '3" employed to denote the 
recombination centers, this recombination rate may be written as 
v tna An/(l + Vgnz/vg P z) + i>, p3 Ap/(l + v gp3 /v un3 ), in which the release 
frequencies are concentration-dependent. If now v llp3 » v„ n3 holds in p- 

t For the minority-carrier trapping range, the lifetime function (t„o + r,,o)Ap/ 
(n + po) would apply for Ah ~ Ap « n + po , as can be shown by use of (100). 



574 THE BELL SYSTEM TECHNICAL JOURNAL, MAY 1960 

type material, it may be a consequence of the condition v lp z > v t „z for 
hole trapping; or, with v tv a < v M , it may imply the recombination 
range. In either case, from (53) the small signal recombination rate is 
substantially vjnsAn provided An is not too small because of strong elec- 
tron trapping in the trapping centers. Lifetime t„ is then T ln3 ~ r . 
While the large-signal lifetime differs in principle, the assumed inequality 
implies r ~ t„ » r p0 if the energy level of the recombination centers 
is not too far from the Fermi level towards the valence band. In general, 
v op s » v„ni gives a small-signal recombination rate equal to 

C*„39t3(poA?i + n Ap)/(p Q + p 3 ) 

and thus a lifetime that cannot properly be associated with either An 
or Ap alone. In intrinsic material, for example, it is An + Ap with which 
a lifetime may be associated. 

These considerations suggest the formal representation of "linear re- 
combination" by including in dAm/dt or in both dAn/dt and dAp/dt the 
negative of a recombination rate v n3 An + v^Ap. This procedure is use- 
ful in deriving results in forms that apply symmetrically without ref- 
erence to conductivity type. For the p-type case here discussed, v nZ and 
j/ P 3 equal v ln z and (n /po)vin3 ■ The one that corresponds to the majority 
carrier can usually be set equal to zero for sufficiently strong extrinsic 
material. 

3.2.3 Negative Photoconductivity 

Under certain conditions, optical generation with excitations involving 
trapping levels will cause a decrease in conductivity below the thermal- 
equilibrium value. 51,52 This negative photoconductivity will be con- 
sidered for a simple model — that of two types of centers, of which one 
gives trapping and the other only recombination. For traps of the ac- 
ceptor type, (94) gives dAn/dt and dAp/dt, except that suitable genera- 
tion and recombination terms must be included. From Section 2.2.3, 
generation terms are respectively Ag„ and Ag p ; and, from Section 3.2.2, 
the linear recombination term -{v n3 An + v pA Ap) may be included in 
both equations. For simplicity, direct recombination and the quadratic 
terms will be neglected, and the concentrations evaluated for the steady 
state. The result is 

= (Vg P l ~ Vna)Ag n + (Vtnl + Vpnl + Vns)Ag p ^^ 

Al + (v ll3 + V p 3){v nl + V op i) + V n3 V lp x + V,oVtnl 

with Ai defined by (62); a similar expression for An is obtainable by 
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interchanging subscripts n and p, a transformation that does not change 
the denominator. 

This result verifies the conclusion that Ap may be negative as a result 
of excitations from the traps to the conduction band in conjunction with 
recombination, and similarly for An, with excitations from the valence 
band to the traps. As a simple case, consider p-type material with v v % 
zero. If there is also trapping only of electrons and excitation only of 
electrons from traps to the conduction band, then C v \ and Ag p are zero, 
whence An is zero and Ap is —Ag„/v g „i , with Ag n equal to Ag c i from (58). 
Recombination, in this case, produces negative Ap, which compensates 
the reduction by the excitation of the concentration of (negatively) 
charged traps, and the effect tends to be enhanced with deep traps of 
small capture cross section. 

3.2.4 Further Theory with an Application to Experiment 

Illustrative application will be made to observations of Hornbeck and 
Haynes on electron trapping in p-type silicon. 54 In this work, techniques 
were devised to measure the various time constants in the decay of pho- 
toconductivity, which, for certain samples, covered a range of about 10 
in relative value. Evidence for two trapping levels was found, and elec- 
tron capture cross sections and energy levels were estimated from the 
data, the model employed being that of two types of traps that capture 
only electrons, a lifetime being associated with recombination in cen- 
ters of another type. The sample f for which there is most detailed in- 
formation exhibited a 20-microsecond photoconductive decay, attributed 
to recombination, for sufficiently high injection levels; a decay of time 
constant about 10 milliseconds, attributed to decay in comparatively 
shallow traps that were initially filled in concentration of 2 X 10 cm - ; 
and decay in deep traps that were initially filled in concentration of 
10 13 cm -3 whose time constant varied from 1 second for the traps nearly 
full to 260 seconds for the traps nearly empty. Both types of traps are 
"deep" traps, as defined in Section 2.2.1.2. The present theory will be 
used to calculate the upper limits for the hole-capture cross sections im- 
plied by this model, and it will be shown how the conclusions are modi- 
fied if an alternative model is assumed. 

In outline, the general procedure here employed involves first assign- 
ing trial values to the energy levels of the traps, and then calculating 
expressions for decay constants from the equations, suitably linearized 
for particular ranges. These decay constants are roots of algebraic equa- 

t Data and results for sample 223B are given in the text and various figures of 
the Hornbeck and Haynes paper. 54 
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tions and, assuming them well separated, may be obtained as the 
magnitudes of ratios of successive coefficients. The coefficients are homo- 
geneous expressions in capture and release frequencies (and also "con- 
straint frequencies," if multiple-level traps are involved) — that of the 
highest power being unity, followed by linear, quadratic and higher- 
order forms for the successively lower powers. With assumed trapping 
levels and known equilibrium concentrations of carriers and unoccupied 
traps, the coefficients provide corresponding homogeneous forms in the 
capture coefficients. These coefficients remain to be found from observed 
decay constants. To each product of capture coefficients that occurs, a 
number of products of frequencies generally contribute, but, for the par- 
ticular semiconductor material and trapping model, these usually differ 
by orders of magnitude and a single one predominates. With this con- 
siderable simplification, decay constants can be expressed in terms of the 
frequencies so that physical mechanisms involved can often be readily 
identified. If a sufficient number of distinct decay constants are known 
from experiment, the energies of the trapping levels may also be deter- 
mined. The consistency of the assumed trapping model may then be 
checked; the energies found should clearly not differ from those assumed 
by so much that the particular simplified expressions employed for the 
decay constants do not apply. 

Consistent with the notation here employed, the deep traps may be 
assumed to be of the donor type and the shallower traps of the acceptor 
type. Trial values of the energy levels will be taken as 0.23 ev below the 
Fermi level 8 hi intrinsic material and at 8. These levels are approxi- 
mately 0.78 ev below the conduction band for an energy gap of 1.10 ev 
and at midgap, substantially in accord with the locations determined 
by Hornbeck and Haynesf The values 72,73 at 300°K of 1500 and 570 cm 2 
volt -1 sec -1 for the electron and hole mobilities and 1.73 X 10 20 cm -6 for 
n/ give n„ = 4.3 X 10 B cm -3 and po = 4.1 X 10 14 cm -3 for the 27-ohm-cm 
p-type sample, with m = pi = 1.32 X 10 10 cm -3 , n 2 = 2.8 X 10 6 cm" 3 
and p 2 = 6.1 X 10 13 cm -3 . 

For the two kinds of traps with recombination at the rate v n3 Aw = 
An/r 3 only in other centers of the idealized type discussed in Section 
3.2.2, the outlined procedure applied to the equations written for the 
linear small-signal range gives the longest decay time t m as 

T m = T 3 + Tgnl + T gn2 + (t u „i/ T tn l) T 3 + ( T„ n2 / Tlnl)^ . (107) 

The fourth and fifth terms represent recombination with multiple trap- 

t An energy gap of 1.10 ev at 300°K is employed rather than 1.00 ev as in Ref. 
54. The trial values employed originated in a two-level analysis (which later 
appeared inapplicable), the n 2 being that for 91 = 1.15 X 10 13 , or p /9l = 0.87. 
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ping in the shallower and deep traps, respectively, and the latter pre- 
dominates for the case under examination.! In (107), T uni and T„ n2 are, 
of course, the "effective" equilibrium release times, and are not cor- 
rectly interpreted as the physically proper ones, J namely (C„ini) _1 and 
(Crf*)"*. 

An upper limit for the coefficient C p2 for the capture of holes by oc- 
cupied deep traps may be obtained by assuming recombination in deep 
traps only and then calculating t x from the linear small-signal equations. 
The result is 

Cfnl + Vg P 2 fio + V op2 /Cn2 (1(\Q\ 

Too = = I ~ Tgpl \i-yfo) 

VlnlVgpl POfgpi 

in the simplified form obtained by the outlined procedure. The final ap- 
proximation on the right § applies with p = 10 cm - and r M = 260 
seconds, provided merely that C„ 2 »4X 10 - cm seconds -1 holds. Then 
C P 2 = 8 X 10" cm seconds - follows, since po + y-i is 4.7 X 10 cm - ; 
and the cross section for hole capture A p2 , obtained by dividing by 
mean thermal velocity, is 8 X 10 -25 cm 2 . Recombination in the deep 
traps that gives t k cannot account for the observed decay. It can be 
shown 1 1 that the decay time for the traps nearly full would then be large 
compared with r (/i)2 = 2G0 seconds rather than 1 second. The actual 
A„2 may thus be considered small compared with about 10 -2 cm". 

For an upper limit to C pi , recombination only in the shallower traps 
is assumed. For this case, the rather lengthy general expression for t w 
simplifies to give 

r M = r„ n2 + (r nn o/r tn ^[T tnl (l + Wtw)]. (109) 

The contribution T fln2 is the time constant for the initial decay in the 
deep traps, obtainable as the longest decay time from the equations 
linearized for nearly full deep traps and nearly empty shallower ones. 
This release time represents recombination of electrons in the shallower 

t See footnote 20 of Ref. 54. 

j The t of Ref. 54 should he identified as t b „\ and r gn i ■ For the deep traps, 
St is accordingly (/in + res) -1 , which increases with the p-type conductivity and 
is not a property of the traps alone. The formula employed for locating trapping 
levels relative to a band edge holds if r in it is the physically proper release time. 
With t b the "effective" release time, it holds only if no is negligible compared with 
«i or «o • In Equations (1) of Ref. 54, dn/dt lacks the term dinoA,,. 

§ The equilibrium t,, of (60) and the lifetime t«_ of (92) written for the deep 
traps, for which 9l 2 * ~ O.lDl™ ~ O.lpo , also reduce to t BP 2 if C n i is not too small. 

|| The lifetime fj of Section 3.2.2 evaluated for the deep traps nearly saturated 
is, by use of approximations for near-saturation of Section 3.1.2, found to be 
given by f •> ~ ((' P 23ln)~ 1 An 2 /(n + «2)po2>T tfp2 for (ho + n->)p y> A« 2 » (n + tin) p„ . 
The inequality on the right is equivalent to An ~2> p; the one on the left is largely 
consistent numerically with An « No + «i , for which the shallower traps are 
substantially empty. 
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traps without their recapture in the deep ones. The contribution (r on2 / 
Tt„2)T t ,a represents recombination in the shallower traps with multiple 
trapping in the deep traps. The contribution with T gp i/T gn i as a factor 
represents recombination in the shallower traps with multiple trapping 
involving both levels, f This contribution predominates since, as shown 
in Section 3.1.2, the small-signal saturation of the shallower traps implies 
t„,>\ » Tom • The only capture coefficient it contains is C pi , which is found 
to equal 2.5 X 10 -13 cm 3 seconds -1 for t„ = 260 seconds. The correspond- 
ing hole-capture cross section A pl is 2.4 X 10 -20 cm 2 . The actual cross 
section is small compared with this value if recombination occurs pri- 
marily in centers of a third kind. 

Recombination in the shallower traps can account for the observed 
deep-trap decay. Indeed, as may be expected, if the equations are lin- 
earized for small departures from a concentration p x of unoccupied deep 
traps, then a (longest) decay time 

Td = T tf „2 + (T on2 /T ta 2)[T*„l(l + T„pi/T„„i)](pi/po) 

= r„ n2 +(t.- r on2 )(pi/p ) 

results, which increases from r ffn2 to t m as pi increases from zero to p 
and is of the same form as that employed by Hornbeck and Haynes to 
fit their data. J The observed decay in the shallower traps can also be 
accounted for through C pl . The equations for nearly full deep traps 
and nearly empty shallower ones give r ffP i(l + r»„i/T BB i) ~ T opl for 
decay in the shallower traps as intermediate time constant, T lnX for elec- 
tron capture being the shortest and r„ n2 for the initial deep-trap decay 
being the longest. § The C pX obtained by setting t (;p1 equal to 10~ 2 second 
is 3 per cent smaller than the value obtained from t x and is thus in rather 
fortuitously close agreement. 

If a model with this C p i is to account for experiment, then the assump- 
tion that the shallower traps are two-level traps, which gave the ob- 
served lifetime II of 20 microseconds through recombination in the higher 



t The quantity in brackets in (109) can be shown, from (53), to be the t of (65) 
for the shallower traps in the p-type material. Thus, to itself may be said to entail 
multiple trapping through (t op i/t ,,i)t,„i , the major contribution to r in the 
minority-carrier trapping range. 

J The interpretation differs, since t„„ 2 is an "effective" release time. In the nota- 
tion of Hornbeck and Haynes, p\/j3o is 1 — y. 

§ Note that t op i results also from assuming filled deep traps and negligible 
An, or Ap ~ po + Art. 

|| For 33 microseconds, as given in Table I of Hornbeck and Haynes, 61 cal- 
culated capture cross section for recombination would be smaller in proportion. 
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level, seems necessary. Otherwise, a near-saturation lifetime much larger 
even than r, /p i would obtain. f From (60), (65) and (66) written for the 
higher level, the recombination lifetime will be t ~ r„ throughout the 
small-signal range and with no small-signal saturation if p » Cans/ 
C p3 = p 3 * and p » C n3 9li/C p3 = (C^Tno)' 1 hold. The latter condition} 
is A P 3 » 1.2 X 10 -17 cm 2 . The coefficient C„s and cross section A n3 for 
electron capture in the higher level are found to equal 2.5 X 10 -8 cm 3 
second -1 and 2.3 X 10~ 15 cm 2 . 

For the further analysis on the hypothesis that this model applies, 
the energy levels are properly treated as unknowns. The contribution 
Typiirtni/Tyni) to the time constant for decay in the shallower traps equals 
T^Ttno/Tgni) — r M (no + n 2 )/po from the expression f or t„, and will ac- 
cordingly be small compared with 10~ 2 second for deep traps sufficiently 
deep so that w + n 2 « 4 X 10 8 cm -3 holds. § Then r„a is 10 -2 second and, 
with pi « po , C P i has the value already found. An additional datum is 
available from experiment, namely the decay constant for the straggle 
effect: With the shallower traps nearly filled, multiple trapping results 
in an extended tail in the distribution of carriers from an injected pulse 
that are caused to drift past a fixed detector, at which the decay with 
time is measured. As shown in Section 3.4.3, the decay constant is the 
"straggle constant" v x , which is substantially v onX + v lp i + v up \ for 
p » n . Since v opX is 10 2 second -1 , the observed value, 2 X 10 second - , 
is to be equated to i» ff »i <^> C„\n\ . With this result, the value for t op i 
and 9li — no ~ 9ti = 2 X 10 12 cm" 3 , the equation for t m contains only 
C„i or ni and n 2 as unknowns. It fixes, say, ni/(n + n 2 ) and thus approxi- 
mately the separation between the energy levels, but there are not suf- 
ficient data with the model assumed to determine each level separately. 
It appears, however, from measurements relating to deep traps in sam- 
ples of various conductivities, || that the location considered for these 
traps is substantially correct. With the trial value 2.8 X 10 6 cm -3 of n 2 , 
the value obtained for C nX is 1.2 X 10 -6 cm 3 second -1 , and the value for 
ni is 1.7 X 10 10 cm -3 , corresponding to an energy level for the shallower 
traps 0.007 ev above the trial location at the Fermi level in intrinsic ma- 
terial. With Vg n2 = 1 second -1 , the value obtained for C„ 2 is 3.1 X 10 - 
cm 3 second -1 . The cross sections A, a = 1.1 X 10 -13 cm 2 and A n2 = 



t With appropriate notational changes, the result in the footnote on page 577 
for (recombinative) deep traps applies to the shallower traps. 

t This controls if the higher energy level is further than about 0.42 ev from the 
conduction baud. 

§ This condition holds by a factor of about 10 2 for the trial value of n« . 

|| See Fig. 13 of Ref. 54. 
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2.9 X 10 -14 cm 2 that result are half an order of magnitude smaller than 
the ones calculated by Hornbeck and Haynes.t 

It seems likely that the photoconductivity under illumination intense 
enough to give the shorter decay times was quite appreciably offset as a 
result of heating of the sample. t Significant error from this source seems 
unlikely, however: The time constant for cooling was very probably com- 
parable with the longer decay times, and measurements concerning these 
were made with considerably less intense initial illumination. § There 
were presumably no pronounced effects of nonuniform generation in the 
thickness of the sample, the generation rate at the dark surface being at 
least 40 per cent of that at the illuminated surface, as calculated from 
the diffusion length for the shortest decay time.|| 

Work has been done towards the identification of the impurities in sili- 
con that occasion these trapping effects. ' It might be noted that the 
energy levels suggest gold. 56,57 But there is evidence that gold gives a 
single center with two (or possibly more) levels, and such a center can- 
not account for the saturation of the shallower traps at a concentration 
less than that of the deep traps. Consider the assumption that a two-level 
model does apply, with shallow traps only partly filled in the experi- 
ments. Then t op i applies for the decay at the shallower level and rmx = 
1/C„i(9l — Wo) for the time constant of 20 microseconds observed with 
the spark source. If It follows that C„i is 4 X 10 -9 cm 3 second -1 for the 91 
of 1.15 po and the negligible w that the trial levels give. But, with this 
C,a , the initial Aw immediately after the steady illumination that is shut 



t Compared with the value from ten samples that they calculated in connection 
with Table I of their paper, .4 „ 2 is one order of magnitude smaller. 

t Perhaps this heating accounts for apparent concentrations of normally 
empty traps determined from Fig. 4 of Ref . 54, which are about 0.7 of the values 
quoted. 

§ Buck 71 has found a positive temperature coefficient of resistance in 38 ohm-cm 
p-type and 350 ohm-cm n-type silicon of 0.8 per cent per °C at room temperature, 
and has observed time constants for the cooling of the samples, similarly supported 
by wire leads and of comparable size and geometry, of the order of 100 seconds. 
The thermal time constant equals the heat capacity divided by the thermal dis- 
sipation constant, or power input per unit temperature elevation. For the sample 
here considered, 0.2 cm square and 2 cm long, power input is 8.7 X 10~ 4 watt for 
10 lfi photons per cm 2 second absorbed, since 1 microwatt corresponds to (5.1 X 10 18 ) X 
photons per second of wavelength X, and effective X for the tungsten illumination 
is about 9 X 10 -5 cm. The dissipation constant for a temperature elevation of 1°C 
with this power input in conjunction with the heat capacity of the sample of 
0.14 joule per °C gives a thermal time constant of about 160 seconds. Haynes 78 
has estimated a temperature elevation of no more than a few degrees for the 
more intense illuminations employed; heating of 3°C would decrease conductivity 
by an amount comparable with the total photoconductivity of Fig. 4 of Ref. 54. 

|| This diffusion length is 0.17 cm. In measurements on n-type silicon, a silicon 
filter and a constant-temperature enclosure were used. 

If Note that t,,,/ and t ( „i ■ 1/C„i(9l — Ho — po) are the times for electron 
capture at the shallower level respectively for filled and empty deep traps. 
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off by the shutter is two orders of magnitude smaller than the apparent 
.saturation value 2 X 10 12 cm -3 of An, and thus appreciably smaller than 
the initial value 6 X 10 n cm -3 of An estimated from the initial con- 
ductivity change, f Also, the apparent saturation value of An is propor- 
tional to Ag, which is not the fairly well-defined saturation observed.! 
Moreover, the value of C nX gives a decay constant for the straggle effect, 
evaluated as v gn \ + v ov \ , equal to 157 second - and the same two orders 
of magnitude smaller than that observed. Inconsistencies largely similar 
result from the assumption of two-level trapping with recombination 
entirely in other centers. 

No evidence has, indeed, been found that the observed trapping 
effects result through metallic impurities or through lattice defects pro- 
duced either by mechanical deformation or by bombardment with high- 
energy electrons. 4,76,77 ' 78 Present indications are that the deep traps in 
silicon are associated with the presence of oxygen as an impurity; ' 
but these traps, as well as the specific reactions instrumental in their 
formation, have not yet been physically identified. Concentrations of the 
traps and of certain donor centers due to oxygen 9 ' have been found to 
be correlated.! Both traps 76 and donors 85 ' 86,87 are much more numerous 
in crystals grown (from quartz crucibles) with rotation of the seed than 
in those grown without, may be considerably increased in concentration 
by comparatively prolonged heating at 450°C, and may be largely re- 
moved quite rapidly by heating at temperatures above 500°C.|| Concen- 
trations of the shallower traps do not exhibit this dependence. The cor- 
relation is qualitative in that donor concentration is the more dependent 
on heating at 450°C; appreciable trap concentration may occur in an 
untreated crystal grown with rotation, and may assume a value con- 
siderably smaller than the donor concentration after heating at the lower 
temperature. 79 ,80,88 H It should further be noted that, while these observa- 
tions have been mostly confined to n-type silicon (because the donors 
tend to convert p-type to n-type), observations concerning the deep 
and shallow traps which occur in p-type silicon indicate that a common 
mechanism is operative. 

t See Figs. 4 and 5 of Ref. 54. 

j The steady-state equations give initial concentrations An = r,„i'(\ + T ap \/ 
7-„„i)A<7 and Art = T„ P \Ag. The apparent saturation value of Ail is the sum of 
these concentrations, since An is trapped rapidly, with time constant t,,,/. 

§ Determinations of oxygen content from infrared absorption at 9 microns in 
combination with resistivity measurements on crystals heat treated at 450°C and 
1000°C have shown that formation of these donors is associated with oxygen. 8081 

|| There seems to be an indication that the trap concentration is increased by 
water vapor but not by oxygen in the gaseous ambient. 

1f Deep traps originally present have been largely removed by heating only 5 
seconds at 700°C and subsequently have been introduced in a concentration larger 
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3.3 The Photomagneloelectric Effect 

The steady-state effect with trapping will be analyzed on the basis of 
the general theory of Sections 2.1.1 and 2.1.2. For the PME field along 
an infinite slab to the faces of which the applied magnetic field is parallel 
and the y axis perpendicular, (28) and (32) give 

E x = tT\l m - BId*) = o-~\l x + deD'dAm/dy). (Ill) 

The total short-circuit current per unit width of slab along the magnetic 
field is accordingly given by 

Iso = / I x dy = -de D —j— dy, (112) 

and the field along the slab under the open-circuit condition is related to 
Tsc as previously derived. f To evaluate the integral in (112), Am is first 
found from the continuity equation 

d(D f dAm/dy)/dy - A?n/r m = 0, (113) 

which follows from (30) and (36) ; the drift term is either zero or of order 
2 for the short-circuit or open-circuit condition. Since, for the slab, 
l PV = —In,, = lay , boundary conditions arc 

£ — D'dAm/dy = s„iAn = s pX Ap = s m iAm, y = y , 

D'dAm/dy = s, l2 An = s p2 Ap = s m2 Am, y = —yo, 

in which £ is the surface rate of generation of electron-hole pairs by 
strongly absorbed radiation and the right-hand members give surface 
recombination rates. For the linear small-signal case, the velocity func- 
tions s n and s p (with second subscripts "1" and "2" for the respective sur- 
faces) are constants, with 

s,„ = (1 — r n )s n = (1 - r p )s p (115) 

the surface recombination velocity for Am. 

The increase in conductance of the slab is given by 

AG = e I (unAn + n P Ap) dy = e(fx n + n P ) I (r c /r„,)Am dy. (116) 

The second form follows from An/r n = Ap/r p = Am/r m = (R m , with 

than the original one by heating 16 hours at 470°C. 76 It is not yet known whether 
prolonged heating at 1000°C, which prevents appreciable subsequent introduction 
of donors at the lower temperature 82 ' 85 ' 86 ' 87 would similarly prevent the introduc- 
tion of deep traps. 

t See Ref. 11, Equation (39). 



CURRENT-CARRIER TRANSPORT AND PHOTOCONDUCTIVITY 583 

t c = Aa/einn + p p )Ag = (n n T n + vw 9 )/(nn + lh>) ( 117 ) 

a lifetime function that determines the conductivity increase Ac for the 
uniform volume generation rate Ag = (R„, . For the linear small -signal 
case, the lifetime functions are constants and 

AG = e(fi„ + n^TcDo'idAm/dy) \ U V Z'\ 

(Ho) 
= e(ju„ + p p )tc(£ - s„ a Am | y=Wo - s„ l2 Aro |„=_ Wo ) 

follows by use of (113) and (114). 

These results show how the theory previously given for the PME ef- 
fect without trapping" is readily generalized to include trapping by writ- 
ing equations in terms of Aw and the diffusivity D', lifetime function 
r m and surface recombination function s m , and employing suitably gen- 
eralized AG. Experiment may determine Z> ', t„, and s„, . In accordance 
with (115) and the results of Section 2.1.2, each of these gives rise, as 
determined by the trapping ratios, to corresponding quantities for elec- 
trons and for holes. 

In its dependence on trapping and recombination in centers of a single 
type, the PME effect is generally nonlinear if deep traps in the minority- 
carrier trapping range are involved. Then trap saturation occurs in the 
small-signal range, as described in Section 3.1.2, and the lifetime may 
be nonuniform: From the illuminated surface into the slab, it may de- 
crease from a saturation value to a much smaller linear small-signal 
value, a transition value at a given depth being sharply dependent on 
light intensity. 

The influence on the PME effect of trapping as such may be investi- 
gated by assuming traps that may be nonrecombinative in conjunction 
with recombination on the dark surface, or with recombination in the 
volume of the idealized type discussed in Section 3.2.2. With the latter 
procedure, the linear recombination term ( — v„aAn — v p3 Ap) is included 
in the continuity equation. For the linear small-signal case, t„, for the 
traps is thus replaced by f m = [r m ~ l + (1 — r n )v„ 3 + (1 — r p )v p3 ] . 
For p-type material, v p3 is set equal to zero, and r 3 = v„-.r l is introduced. 
Then, for nonrecombinative electron traps of the acceptor type, r n = 
vtni/ivM + >W), >'„ = 0, A/ = [1 - Po9li7(«o + Po)(9li* + nu)]D 
and t„, = (1 + Vmi/vg„i)Ta are obtained by use of the first of (31) and 
(60). If the traps are of the donor type, then r„ = 0, r„ = —v ln2 /v on 2 , 
]) {) ' = [I + 91 2 */(ho + p () )].Do and f m = t 3 are obtained. Essentially the 
same diffusion-length lifetime associated with Z) , namely 



T„ 



- [1 + 9l//(n + Po)]r 3 = K,t 3 , j = 1, 2, (119) 
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results for both types of traps. Thus, minority-carrier trapping increases 
the diffusion -length lifetime and hence decreases PME current. The ef- 
fect is appreciable in cases for which the capture concentration 91/* is at 
least comparable with the equilibrium concentration of majority car- 
riers.! Similar analysis for nonrecombinative majority-carrier traps gives 
a K T which is that of (110) modified by division by 1 + dlj*/n for n-type 
material or by 1 + 9l,-*/po for p-type. Thus, majority-carrier trapping de- 
creases the diffusion-length lifetime and increases PME current, but 
this increase is only that for a K T no smaller than (1 + po/n )~ or 
(1 -f- no/po) -1 , respectively. 

Capture cross sections, concentrations and energy levels of traps may 
be found from suitable PME and photoconductivity measurements at a 
single temperature. Theory for trapping and recombination in traps of a 
single type, which holds whatever the method be for determining diffu- 
sion length Lo and lifetime t () will be considered first; while the PME 
method has certain advantages, any one of a number of other methods 
may also be employed.! In view of the fundamental restriction of (53) 
to which the four capture and release frequencies are subject, it will be 
convenient to deal with the capture frequencies v t „i and v lp x and the 
capture concentration 9li* of (63) as independent parameters. To de- 
termine these parameters, three quantities must be measured. Suppose, 
for example, that from suitable linear small -signal measurements, r , t c 
and the lifetime t-> of (92) for decay of photoconductivity are known. 
Solving (60) and (92) for t,„i , t, p i and 9li* gives§ 

Ttnl = [(t 2 — Tp)/(T» — T u )]r , 

r w - [(r 2 - Tn )/(T 2 - r a )]r 0> (120) 

9li* - (no + Po)(r, - T a )/(T n + r„ - r,), 

in which r„ is denned by 

Ta = (WoTfl + Pi>Tp)/(n + p ). (121) 

Then, with 

t„ = [(fi n + f*p)noT e — n P (n + Po)To\/(fJ. n n — m p ?>o), 

(122) 

T„ = [Hn(n + Po)tq — (fl„ + Hp)poT c ]/(t*nn ~ /ipPo), 



t It can he shown that, if different types of traps are present, the 9l>* in (119) 
is replaced by the sum of the respective capture concentrations. 

X See, for example, van Roosbroeck and Buck. 89 

§ Note, from (60) and (92), that n is larger than t„ , t„ and t„ , and smaller 
than r n + t p . 
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subject tot jx„n 7* n P po , which are obtained by solving the equations 
defining t and r , (35) and (117), t m , T tpl and 9li* can be found from 
experiment. It should be noted that (117) does not represent an advan- 
tageous method for determining r c ; an indirect method will be given 
that obviates the necessity of knowing the intensity of absorbed radia- 
tion. 

An additional independent datum is required to determine C„i , C p \ , 
Dli and Wi or pi . As (41), (42), (50) and (63) show, t„ , r p0 , or equi- 
librium concentration of empty traps would serve. % Thus, a measure- 
ment involving the saturation range is required in addition to those in 
the linear small-signal range. It is, in fact, desirable that two such meas- 
urements be made, for reasons that will be discussed. Suppose, for ex- 
ample, that there is small-signal trap saturation in p-type material and 
that the decay frequency C P i9li = t p0 _1 in the saturation range and 3li — 
w are known in addition to t and t 2 . It follows then, from the first 
equation of (63), that 9li* is v tpX times a known constant: 

9li* = (CpiTi p i) , 

C pl = t p o _1 (91i - no)" 1 = (1 + Pi/po)C p i . 

Eliminating w t, p i + poT tn i from t 2 of (92) by use of the third form for 
to of (65) results in an equation linear in r tn i and T tpi , after 9li* has 
been eliminated by use of (123). This linear equation and the one for 
to may be solved for T t »i and t (p i , and, with (50), this solution gives 

no/9li - (1 + Wi/no)" 1 

(po — TOo)/(po + n ) — C p iPq(t 2 — Tp) 
P ° Pot 2 /(Po + n ) — TO 

fti = (3li - no) -1 (124) 

(pa — Wo)/(p + Wo) — C p iPq(t 2 — To) 

To — W T 2 /(po + Wo) — C' p l(po + Wo)(t 2 — To)to' 

C'pi = Cpi — T p0 -1 (9li — w )~'wo/3li • 

With trapping, the PME current-conductance ratio does not deter- 
mine to but depends also on t c (which differs from t because of trapping) , 
and direct determination of t c requires knowledge of light intensity. 

t With the denominator /i„n - n P po equal to zero, to = t c follows, and the 
numerators are also zero. 

J The neutrality condition would serve in cases of trapping by acceptors which 
determine the (p-type) conductivity, for which ii = po — «o holds; but n and 
p — n that differ phenomenologi'cally must generally be considered to ob- 
tain. 
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This circumstance need not, however, really vitiate what in the no- 
trapping case is a primary motivation for dealing with this ratio: What 
is determined independently of light intensity is a relationship between 
ro and r c , and a further relationship between these lifetimes will serve 
to determine both. For C' p i known, the relationship 

C pl = (r^Sli*)" 1 = (r„ + r p - r 2 )/(po + n )(r 2 - r n )r (125) 

obtained from (120) is, with (122), linear in t c and readily solved for 
this lifetime. 

The PME method of the high-recombination-velocity dark surface 
is best employed, since it generally provides better accuracy for the con- 
ductance change than does the thick-slab method which it otherwise 
subsumes as a limiting case.f Optimum slab thickness is about one or 
two diffusion lengths. For large dark-surface recombination velocity, the 
small-signal results for no trapping! give, for the present case, 

/„ = - 6e£L (S 1 + coth 2Fo) _1 

(126) 
= -0(/i» + n v ) (£<A c )(coth Y )AG, 

in which diffusion length L is (Doro)*, Y is yo/Lo and Si is s m iL /D '; 
note that AG now involves r c as a factor. Thus, t c is given by 

2F cothF / 10 -v 

T °-[4/(AG/G )] T0 ' U27) 

where 

d/(AG/G ) = -2y (nn + ix p )I BC /eD AG (128) 

is the dimensionless PME current-conductance ratio. In (127), t enters 
also through Y , and, with (125), both t and r c may be found. Note 
that apparent lifetime t t on the assumption of no trapping, obtained by 
equating 0/(AG/G o ) to [2yo/(DoT T ) *] coth [ya/(D T r ) ], is related to r 
and t c by t t tanh 2 [^o/(^or r ) J ] = (t c 2 /t ) tanh 2 Y , and equals t*/t only 
for the thick slab, for which the hyperbolic tangents are unity. 

If the model that applies is that of nonrecombinative traps with re- 
combination in other centers, then (119) gives the lifetime fo upon which 
the linear small-signal I 8C for minority-carrier trapping depends. For 
AG, (116) holds for the linear small-signal case, for which T c /T m is 
[(1 — r n )n n + (1 — r p )ju p ]/(/x„ + Mp)- The solution for Am is readily 
obtained by comparison with that for the corresponding no-trapping 

f See Ref. 11, Section 3.42. 
t See Ref. 11, Equation (50). 
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case.f With the factor Z/ 2 /D ' replaced by f m , it is found that 

r = K G e(ii n + Mp)r 3 £(cosh 2F - 1) , 19g . 

"'' S x sinh 2F„ + cosh 2F 

is the linear small-signal AG for large dark-surface recombination ve- 
locity, in which Fo and »Si are defined as above, but in terms of L = 
(Dofo)'. The factor K G = (T e /r m )(f m /T 3 ) is t c /t„ in general for p-type 
material, for which f „, is t 3 /(1 — r„) ; for n-type material, K is r c /r p . 
The expression that K G multiplies also depends on trapping, since L 
does. Equations (60) and (117) give, for p-type material,! 

K a = 1 + (6 + l) _1 9ly*/n (electron trapping) 

Ka m i + (i + t-rV/P. (hole trapping) . d30) 

1 + VLj*/p 

For hole and electron trapping, respectively, in n-type material, n and 
p in these equations are interchanged and b = m»/mp replaced by its 
reciprocal. 

With I sc for this model given by the first form of (126) with the re- 
defined Lo , 

tf/(AG/Go) = (Kr/K G ) 2F coth Y (131) 

follows by use of (119), (128) and (129). Apparent lifetime r r is ac- 
cordingly given by t t tanh 2 2/ /(Z>or r ) 4 = (Kq/K t )t 3 tanh 2 F , and 
equals (Kq/K t )tz for the thick slab. As trap concentration increases, 
diffusion length increases and a slab of any given thickness becomes a 
"thin" slab, for which F coth F ~ 1; and 3/(AG/Go) approaches a 
constant value that is independent of the thickness. For example, if 
the half -thickness y is of order (Dots)*, then K T » 1 or 91,-* » n + p 
also gives small F . From the expressions for K r and (130) for K a it is 
found that d(AG/G ) approaches 2(6 + l)n /(n + Po) for electron 
trapping and 2(6 + l)po/b(n + po) for hole trapping, regardless of 
conductivity type. On the other hand, if the slab is so thick that y 2 » 
DoTi^Lj*/(n -f- po) holds, then the condition 91/* » ?? + Po for large 



t In Equation (44) of Ref. 11, Ap is replaced by Am; the Do that appears ex- 
plicitly originates from the boundary conditions and is replaced by Do'; and Si 
and .S'-i are the velocities for Am multiplied by Lo/Do'. 

t Note that K, and K a are equal (for electron or hole trapping) in p-typc ma- 
terial for which n„tlo = n P po holds. 
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trap concentration gives f r 3 /r r equal to K T /K G 2 or (b + l) 2 n 2 / 
(»o + Po)9l>* for electron trapping and (b + l) J po/6 2 (/?o + Po) for hole 
trapping in p-type material, with similar results for n-type obtained as 
in connection with (130). 

With the aid of suitable saturation-range measurements, the cross 
section, concentration and energy level of the nonrecombinative traps 
are readily found. For traps of the acceptor type, C„i , 9li and n\ or pi 
are to be determined, and these can easily be calculated from values from 
experiment of 9li — rk , saturation -range lifetime t 3 , lifetime t „i = 
[C»i(n + Wi)] _1 for the traps nearly full, and any one of to , r M or 9li*. 
These last three quantities are not independent; from (107), lifetime t„ 
for the traps nearly empty is [1 + (9li — n )/(n + ni)]r 3 = 
[1 + 9li*/w ]r 3 . Measurement of 3/ (AG /Go) serves to determine f : By 
means of (119) and (130), K T /K a may be written as (6 + l)n [n + 
Po — (po — 6n )r 3 /fo] -1 for electron trapping in p-type material, or as 
an analogous expression for hole trapping in n-type, so that (131) in- 
volves only to as unknown. 

3.4 Transport of Injected Carriers 

3.4.1 The Linear Differential Equations 

The general differential equations of Sections 2.1 and 2.2 are here 
specialized to the linear small-signal case of trapping (and recombina- 
tion) in centers of a single type, for which certain specific transport 
problems will be considered. From (6) through (9) and (46), the linear 
continuity equation for centers of the acceptor type is 

dAm/dt = Dp div grad Ap + D» div grad An 

— Vp-grad Ap — v„-grad An -\- Ay — (R,„ 

(132) 
= dAp/dt = Do div grad Ap — D n div grad An 

— Vo-grad Ap + v„-grad An + A<7 + ^nAp -f- v v >An, 

the first form being that which applies for the linear case in general. The 
diffusivity and velocity with minority-carrier subscripts of those defined 

»y 

t It has been shown by Amith 62 - 63 - 90 that, for minority-carrier trapping in the 
(hick slab, T 3 /r r is proportional to 91/ -2 for large 31, , if K r is taken as unity. This 
dependence obtains in the intermediate range in which 91,* is large compared 
with minority-carrier concentration n or p„ but small compared with p or no so 
that the change in diffusion length may be neglected. For majority-carrier trap- 
ping in general, tj and r r are substantially equal in this range. 
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t) p = n D /(n + po), D n = p Do/(n + p ), 

(133) 
v p = n Vo/(n - p ), v n s — PoVo/(mo - po) 

are, in sufficiently strongly extrinsic material, substantially the minority- 
carrier diffusivity and velocity, as are D and v ; while those with ma- 
jority-carrier subscripts are comparatively small. The linear equation 

dAn/dt = v n Ap + v-, 2 Aii (134) 

holds for interactions with the traps. Eliminating An from (132) and 
(134), substituting from (133) and making use of (51), (02) and (63) 
results in 

d'Ap/df - Dodiv grad(dAp/3/) + v„-grad(dAp/d/) + v s dAp/dt 

- i> D D div grad Ap + i» r v -grad Ap + A t Ap (135) 

= dAg/dt + (is - vt P i)Ag, 

with i/„ defined by (83), and v D and v v by 

r B » [1 + 9li*/(«o + p )](v g „i 4- ^ pl ) (136) 

Vt = V an \ + Vg P i 4" (v gn i — C 0P l)9ll*/(^O _ Po)- 

The frequency v„ will be referred to as the "straggle constant". It is readily 
shown that the linear differential equations that An and An satisfy are 
entirely similar to (135) except for suitable modifications of the right- 
hand member; all the concentrations satisfy the same equation if there 
is no volume generation. For An it suffices to replace v tp \ where it occurs 
explicitly by v t „i , while for An only the generation term (v, n i — vi p i)Ag 
occurs, dAg/dt being absent. It is also readily shown that linear re- 
combination in other centers can be taken into account by adding 
v„3 4- p,* to the coefficient v s of dAn/dt and dAp/dt and (v, — vt„i)v n t + 
{v, — vipi)v P 3 to the coefficient A x of An and Ap. 

3.4.2 Steady-State Transport; Reverse Drift 

A simple case that yields qualitative information of interest is that 
of injection into a filament in the steady state with applied field. For 
this case, 

voDtffAp/dx 2 - v&dAp/dx - A,Ap = (137) 

is to be solved for, say, Ap zero for distance x along the filament nega- 
tively or positively infinite and continuous at the origin at which there 
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is carrier injection with zero injected total-current density. Equation 
(137) is easily shown to be equivalent to (30) and (31) specialized for 
no volume generation and acceptor centers only; v D and v v are ( — Pn) = 
v, — v,„i times D '/D and vo'/v , respectively, and t p , from (60), is 
(v s — pi p i)/Ai . 

The solutions in the semi-infinite regions separated by the origin are 
e r ' r and e 2X , where ?'i and r 2 are given by 



^ {v v v /v D Df\ 
\— Ai/v v v /, 



(138) 



as obtained from (137). The case of recombination without appreciable 
trapping 64 presents no unfamiliar features; the approximation given, 
which is that for Ai small, as may result from one of the capture coeffi- 
cients small, will accordingly be considered. The magnitude of n is thus 
large compared with that of r 2 . With the condition v > 0, which may 
be assumed without loss of generality, e ix gives the familiar sharply 
varying field-opposing solution to the left of the origin and e r2X gives 
the corresponding gradually varying field-aiding solution to the right, 
provided v v is positive; then, rj and r 2 are respectively positive and nega- 
tive. But negative v v can occur, for which an anomalous behavior ob- 
tains, the field-opposing and field-aiding solutions then being respectively 
the gradually and sharply varying exponentials e 2X and e nx . For this 
case, in the limit of no diffusion, added carrier concentration appears 
only in the direction opposite to that of the ambipolar drift velocity, 
that is, opposite to the direction of drift normally determined by con- 
ductivity type. 

This "reverse drift" associated with trapping may be understood in 
terms of properties of the current density AI of added carriers. From 
(19), added carriers drift in the direction of the total current density, 
or the contribution to AI from drift has the sign of I if noAp — p An 
or Ap/An — po/no is positive; that is, if injection results in proportion- 
ately more holes than electrons than is the case at thermal equilibrium. 
This behavior is, of course, that which normally occurs in n-type mate- 
rial; with no trapping, Ap/An equals unity and added carriers drift 
with or opposite to I according to whether the semiconductor is n-type 
or p-type, with no drift in intrinsic material. ' Thus, the normal be- 
havior requires the conditions that Ap/An — po/no be positive in n-type 
and negative in p-type. It is easily shown, by writing these conditions 
by means of (60) for the steady-state value (1 — r„) -1 of Ap/An, that 
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both are tantamount in the steady state to the single condition, v v > 0. 
This condition clearly always holds for the majority-carrier trapping 
range, while reverse drift results for sufficient minority-carrier trapping 
in not too strongly extrinsic material. From (136), v v > gives 

P0 — »l0 > (V gn \ — Vgpl)Sli*/(v„ nl + Vgpl) (139) 

for p-type material, and a similar inequality for n-type is obtainable by 
changing the sign of each side. Equating the two sides gives the condi- 
tion for no drift, which, for no trapping, holds for intrinsic material. 
From (54), the right-hand side of (139) may be written as 

(n - n!*)9li*/(no + «i*) = (pi* - Po)9li*/(Pi* + Po). 

It reduces to 9li* for electron trapping without recombination, since 
then v„pi and n* are zero. For this case, since 31* equals nov M /v gnl from 
(63), reverse drift obtains if rto/po in p-type material exceeds rt»i/ 
(r t „i + Tg,a), the fraction of the time electrons are free. A similar result 
holds for hole trapping in n-type material. 

3.4.3 Drift of an Injected Pulse 

The differential equation for drift with negligible diffusion and no 
volume generation in one cartesian dimension with trapping by centers 
of a single type is 

d~Ap/dt~ + Vod 2 Ap/dxdt + v s dAp/dt + v v v dAp/dx + AiAp = 0, (140) 

from (135). For a pulse of carriers injected into a doubly infinite fila- 
ment, a suitable technique of solution is that of the bilateral or two- 
sided Laplace transform 91,92 with respect to the distance variable, for 
which the notations 

F(s, U) = [ K e~*f(\ U) dX = 8{/(X, U)} m f(X, U) (141) 

J— co 

are here employed. Dimensionless independent variables 

X s x/L, U a t/r (142) 

are introduced, and with distance and time units given byf 

L = V T, 

r m ( | v 2 | )"\ 

v- = 4[v v (v 3 - v v ) - AJ (143) 

■ 4nt{wba — vt P i)~(po — n )~ 2 
•[(",„i + i»i P i)(po - n )/(v hl i — vt P i)3ii* — 1], 
t The second form for v- follows by use of (52), (62) and (63). 
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subject to the restrictions v 7^ and n 5^ p , the reduced equation 
d 2 Ap/dU 2 + d 2 Ap/dXdU + {dAp/dU 

+ |(f 4- K )dAp/dX + i(r - k" =F l)Ap = 
results, where k and f are the parameters 

K = (2v, - v s )r 

= K»i + '•* + ("0 + poX'tai - v w )/(n - po)]r, (145) 

f = V,T = (»|»i + P„,a + Pjpl + ^pl)?" > 0. 

Coefficient unity for the second term of (144) results from the definition 
of L. The double sign in the last term of the equation results from the 
necessity of defining a real (and positive) t, the upper and lower signs 
applying respectively for positive and negative v . 
Laplace transformation of (144) gives 

(fAi/dl? + (s + t)dAJi/dU 

+ [|(r + *)s + Kf 2 - « 2 =F Dfip - 0. 

As has governed the choices of L and r, the roots ( — Ni) and ( — N 2 ) of 
the associated quadratic reduce to 



(-^) = §{±[( s -k) 2 ± u 5 - (« + r)}, 



(147) 



in which the double sign inside the radical here and in what follows 
relates only to the sign of v. Equation (146) holds for each of the trans- 
formed concentrations, as does (140) for each of the original ones. Gen- 
eral solutions are thus 

A~p = JLA i c-^\ 

__ V (148) 

An = lLr nj Aje- N > u . 
1-1 

From these solutions in conjunction with the Laplace transform of 
dAn/dll obtained from (134), it is found that the r„> are given by 

r*--^<» + Wr) (U9) 

= (Vtnl - Vtpl)/(vinl + Vgm + V„ pl - Nj/r), j = 1,2. 

With vj replaced by Nj/t, these are formally the same as the trapping 
ratios for the decay of photoconductivity given in the first line of (81 ).f 

f Other forms for the r nj obtainable through Laplace transformation of dAp/dl 
from (132) written for one dimension and no diffusion are not similarly related to 
the forms of the second line of (81), though, for s = 0, the Nj/t reduce to the vj 
and all /•„> to those for photoconductivity. 
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The Aj are easily found in terms of the r,,, and transforms Api and Ani 
of the initial distributions. For no carriers trapped initially, the case 
that will be considered, (147) to (149) give 

( A ; / ^\ = \ {1 =F (s - S)/[(s - 2 ± If}, (150) 

in which the parameter £ is defined by 

£ = ( Vll - vm)t = (f. - 2F tpl )r. (151) 

Corresponding coefficients for the solution for Aft are (1 — r,,i)Ai and 
(1 — r„i)Ao , and these are similar to the ones for Ap if £ in (150) is 
replaced by 

7J = (v n — v-ii — 2v 2 i)t = (v, — 2F( B l)r. (152) 

Solutions for an injected gaussian delta pulse are advantageously derived 
as limiting forms as a approaches zero of solutions for the gaussian initial 
distribution whose transform is given by 

A^/«P/L) = 8{|r-*a-V- zW ) = e oM (153) 

for 6* carrier pairs injected per unit area of cross section. From (147) 
through (153), Ap for this initial gaussian distribution is given by 

AP = Ap/WL) = e a2 > 2 - iu( '™ (cosh {Wi(s - k) 2 ± l] 5 ) - (s - £) 

•[(« - k) 2 ± 1]~* sinh {if/[(s - k) 2 ± 11*}), (154) 

from which AN h An/((P/L) is obtained by replacing £ by 77. 

Certain inverse Laplace transforms that are needed are derived in 
Appendix A, and Appendix B includes some details of their use in cal- 
culation of the solutions from (154). Solutions for the initial delta pulse 
at the origin are found to be 

AP = Ap/((?/L) = [e «*-*«+■>*] 

■ Uu - x) + \ ((« - k) j° o [VxJu - X)) 

- j[ iVXiU^X)] Vx( *_ X) } 1 [X(U - X]), (155) 

AA? ■ An/((?/L) = J({ - v )\e' x -'^ +K)u } 

■ f T ° \VX(U - X)} 1 [X(U - X)}; 

for AN ■ A«/((P/L), £ in AP is replaced by y. The modified Bessel 
functions I and h apply for the upper sign in (144) and (146), that is, 
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for v real, while the Bessel functions J Q and J Y apply for v imaginary. 
The term in AP and A.Y with the delta function 8(U - X) = 
voT8(vd — x) represents a contribution that drifts at velocity vq . The 
continuous distributions are confined to the interval ^ x ^ v t, 
1[X(U — X)] = l[x(v t — x)] being the step function that is respec- 
tively zero and unity for negative and positive values of its argument. 
Note that k and 

§(? - k) = (1 - no/Po) _I ("«ni - vt p i)r, 

(15b) 
|(ij — jc) = (po/n — 1) (vtni — vtp\)r, 

£(£ — V) = ("to! - "tpl) 1 " 

are not restricted as to sign. For numerical computation of solutions 
and further analytical study, it is well to transform (155) by eliminating 
X in accordance with 

X = U sin 2 ^0 = \U{\ - cos 6), U > 0, (157) 

which gives 

ap = 6 -**ff+«oo. e, |[i[/( T _ e)r 1 5(7r - 6) 

+ H(€ - {° (ft/ sin 0) I 7 / (^ sin 6) tan *e] 
Jo «/i 

•l[G(7T-e)]} (158) 



Aiv- = Kf- r,)[ e -' u(f + ,(C08e) ] 7o (|cTsine) 

•/o 

•i[e(r - e)]. 

For AN, £ in AP is replaced by 17. The use of as a variable implies 
the step function of (155), while the step function of (158) simply re- 
stricts as defined by (157) to the interval ^ ^ v. 

Interpretation hi descriptive terms of cases of imaginary v requires fur- 
ther analysis. Illustrative cases of minority-carrier trapping in strongly 
extrinsic material, for which v is real and whose interpretation is com- 
paratively straightforward, will be presented first. For strongly extrinsic 
material, since the parameter £ or 77 for minority carriers is substantially 
equal to k, the minority-carrier concentration does not include the term 
with the Bessel function h . If, also, the trapping is nonrecombinative, 
then f = (v + v t )r and k = (v„ — v t )r hold with v = 4v t v a , where 
v t and Vg are v tn \ or v tjA and v gn \ or v ap i , respectively, and refer to the 
minority carrier. Figs. 1 and 2 show distributions of mobile minority 
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2.5 5.0 7.5 10.0 12.5 15.0 17.5 20.0 

DIMENSIONLESS DISTANCE, X = X/V T 

Fig. 1 — Continuous concentration distributions at different times of mobile 
minority carriers from an injected delta pulse for drift with trapping. A strongly 
extrinsic semiconductor and f = 1, k = are assumed. For nonrecombi native 
trapping, these assumptions imply trapping time t ( and release time t„ both equal 
to twice the time unit t. The pulse at the limit of the drift range is attenuated by 
the factor <r»< r ~"> D = e~» p . 

carriers for this case. For Fig. 1, |" is unity and k zero, as for trapping 
and release times equal, f and the continuous distributions are shown 
for different times after injection at the origin of the neutral delta pulse. 
These distributions are led by a delta pulse, which drifts at the ambi- 
polar velocity v . This remnant of the initial pulse is composed of un- 
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2.5 5.0 7.5 10.0 12.5 

DIMENSIONLESS DISTANCE, X = X/ V T 

Fig. 2 — Continuous concentration distributions of mobile minority carriers 
from an injected delta pulse, for drift with trapping, all at time 10V and for re- 
lease time respectively 4, 1 and \ times trapping time. A strongly extrinsic semi- 
conductor and nonrecombinative trapping are assumed; f equals \, 1, £; k equals 
— 4i 0, |; r equals n , hn = hr g , t„ , respectively. 

j Fig. 1 applies more generally: The values of the parameters do not rule out 
recombination, but imply merely v = v, = %n n \ = 2[(v gn i/i>tni)(vtni 2 — vtpi*)]* for 
minority electrons and similar relations for minority holes. 
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trapped carriers and is rapidly attenuated by the exponential factor, 
which is e~ hu for the particular values of the parameters. The abrupt 
fronts of the distributions for the shorter times result from most carriers 
having been trapped only slightly — at least once, but not much more. 
A relative maximum f appears in the case of the figure for times greater 
than 2 } r. For the longer times, the abrupt front disappears as a result 
of multiple trapping. Furthermore, there is a reduction of apparent 
mobility : The maximum ultimately drifts at velocity \vo , the fraction 
of vo equal to the fraction of the time the carriers are free; this equality 
will be shown to apply for nonrecombinative minority-carrier trapping 
in strongly extrinsic material. 3 This limiting behavior sets in rather 
slowly, as the distribution for I = 20t shows; its maximum occurs some- 
what beyond the middle of the drift range. 

Fig. 2 shows distributions all at time t = 10t for release time respec- 
tively 4, 1 and £ times trapping time. The increasing areas under these 
distributions are associated with decreasing fractions of carriers trapped ; 
it will be shown that, for nonrecombinative trapping, this trapped frac- 
tion rapidly approaches the fraction of the time carriers are trapped. 
The distributions have maxima appreciably beyond the respective values 
for large U of one-fifth, one-half and four-fifths of the drift range, and 
the distribution for the comparatively small release time still exhibits 
a high abrupt front at the time 10t. 

The parameters on which the solutions depend have certain general 
properties. From the first forms for k and £" of (145) and the definitions 
of r and v of (143), 

r = {k ± 1 + 4A!T 2 )* ^ (k 2 ± l) 3 (159) 

follows. The inequality sign is associated with recombination, A x being 
zero for nonrecombinative trapping. The parameter f is real and never 
negative. For v imaginary, so that the lower sign applies, a similar calcu- 
lation gives 

K 2 = 1 + (p 2 a - 4A0r 2 ^ 1, v < 0; (160) 

the condition v a ~ — 4Ai ^ implies real decay constants and holds from 
(85). For v real, k is not restricted. For example, for nonrecombinative 
trapping in strongly extrinsic material k is \[{v„/vt) — (vt/vg) 1 ] and 
can be zero (as for Fig. 1) or have any positive or negative value. Thus, 
f ^ 1 holds for v real and f ^ holds for v imaginary. With (151) and 



t Expressions for dAN/dX = {%U sin 9) -, </A.V/</0 from the Maclaurin's ex- 
pansions for the modified Bessel functions reduce, for f = 1 and k = 17 = 0, to \e~'* u 
andj(l — 5C 2 ) e -1£/ at the origin and at the end of the drift range. 
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( 152) , which define £ and tj, it is readily seen that £ ~ -q ~ k holds under 
the condition for the validity in general of the lifetime function r n ~ t p 
of (75) ; and that £ or 17 is approximately k for n » p or p » n , respec- 
tively. These properties are consistent with the easily verified relation- 
ship, 

woE — Pol = (no — Po)k. (161) 

Also, for nonrecombinative trapping, the parameters do not depend on 
the capture coefficient and f = (V ± l) 3 equals £ or rj according to 
whether electrons or holes are trapped. 

Qualitative and physical distinctions are related principally to the 
signs of v and k. The general condition for real v is 

(vim + Pt P i)(Vv ~ 'h)/(u M - v tpl ) > 91,*, (162) 

from (143), and v is imaginary if the inequality is reversed. From (54), 
the left-hand side may be written as (p + ni*)(p — n )/(po — rii*) = 
(pi* + no)(p — n )/(pi* — no). The condition k > is, for p-type 
material, 

p — >U) > (v t „i — Pt P i)(po + nu)/(jv,i + v gpl ), (163) 

from the first equation of (145) ; changing signs of both sides or reversing 
the inequality gives k > for n-type or k < for p-type. From (54) 
and (63), the right-hand side may be written as 

(1 +n n /po)(po- ».i*)3Xi* = (l + po/noHpi* -n u )9li* 
»o + ni* pi* + p 

If recombinative trapping is excluded, then v is evidently real in the 
limit of strongly extrinsic material. Xote, for example, that as p% in- 
creases indefinitely, v lpX approaches zero while v tn \ approaches C nl 9lx = 
r„o _1 , so that v approaches 2C „i(»i9li)*. Also, since v gpl increases indefi- 
nitely and v„,a approaches CniKi , k becomes positively infinite. If v tn \ 
and v tP i are equal for given p Q , then v tn i > v lp i holds for all larger p . 
Suppose first that these capture frequencies are equal for some po of the 
n-type material. Then, as p it decreases from a large value, k docs like- 
wise. With v M > v tp i , (162) and (163) show that k decreases to zero 
and becomes negative for v still real. Further decrease of p results in k 
becoming negatively infinite, since v approaches zero as p — no ap- 
proaches (v M — r«pi)9li*/("^i + "i P i)> following which k increases to 
— 1 with v imaginary. While v is not defined for intrinsic material, the 
equations show that, in the approach to the intrinsic limit, v is imaginary 
and k approaches ±1, the sign being that of {v M — v tp i)/(n — po). 
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It is evident from (162) and (163) that, for n-type material and v tn i > 
vtpi — that is, in the majority-carrier trapping range — v is imaginary 
and k positive, hence greater than unity. With increasing no , v t ,a and 
v tp i ultimately approach equality, v approaches zero and k becomes 
positively infinite. For p-type material and v tn i > v lp i , or the minority- 
carrier capture range, it is likewise evident that positive k implies real 
v and that imaginary v implies k < — 1 . This latter case includes reverse 
drift. From a result of Section 2.2.1.2, since v t „i equals v lp i in n-type 
material, v„„i equals v gpl in p-type. Decrease of p from the value for 
Vgni = v up i gives v„„i > v apl , and, from (139), as p a approaches n , neg- 
ative v v occurs. In the limit of negatively infinite k, for which v is zero 
(with Ptnj 5^ vtpi), Vv is in general positive. Thus, the reverse-drift range 
is in general the portion of the minority-carrier capture range of imag- 
inary v that results if a certain infinite range of large negative values of 
k is excluded. For nonrecombinative trapping, (139) and (162) both 
yield | p — n \ < 3li*, so that the two ranges coincide. 

If the capture frequencies are equal for some po of the p-type material, 
the initial decrease of k as pa decreases from a large value still obtains; 
but the po for equal capture frequencies is approached for v still real 
and k positive, and k again becomes positively infinite as v approaches 
zero. Imaginary v results with further decrease of p so that v lpl > v tIl i 
results, k decreasing from large positive values to unity and then from 
— 1 to large negative values, the majority-carrier trapping and minority- 
carrier capture ranges, respectively, being realized (for hole-capture 
frequency the larger) in p- and n-type material. It is easily shown, as 
before, that the reverse-drift range applies, with recombination, for a 
finite range of negative values of k less than — 1 in this minority-carrier 
capture range of imaginary v. Increase of n beyond the value for nega- 
tively infinite k given (as before) by n — po = (v tp i — v ln i)/(v t p\ + v tn i) 
results in real v, with which k ranges from large negative values and be- 
comes positively infinite as the material becomes strongly n-type. 

For nonrecombinative trapping, k for real v and strongly extrinsic 
material is, as has been noted, a positive or negative constant. For elec- 
tron trapping, for example, v is 2C n i(n 1 9Zi)* and k is i[(ni/3li) — 
(9li/ni) ] in the limit of large po ; k is positively infinite in the limit of large 
n . As po decreases in the minority -carrier capture range of real v, k be- 
comes negatively infinite as p — no approaches 9li*. With further de- 
crease of po , k increases to — 1 in the minority-carrier capture range of 
imaginary p, which (for nonrecombinative trapping) is the reverse-drift 
range. In the majority-carrier trapping range in general, v is imaginary 
and k greater than unity. 
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Maxima for large U of each of the continuous distributions for cases 
of real v occur substantially together. From (158), the for a maximum 
is found to be given byf tan 9 = — k~ 1 , which gives J X/U = x/vot = 
^[1 + «/(k 2 + 1)*]. For nonrecombinative traps, real v implies minority- 
carrier trapping with positive v„ and (V + 1) equal to f. From (63), 
(136) and (145), X/U accordingly reduces to v r /{v t + v g ), which is 
[t m - n T an i/(p - »o)]/(r t „i + T g „i) for electron trapping or [r, pl - 
Por BP i/(no - po)]/(r tp , + T gp i) for hole trapping. Hence X/U, the factor 
by which the apparent mobility is smaller than the magnitude of the 
ambipolar pseudomobility, 10 is in general less than Tt/{r t + t„), the 
fraction of the time minority carriers are free; but X/U is substantially 
equal to this fraction 37 under the condition | no — po | J> 9li*, obtained 
by use of (63). As | n — po | approaches 9Zi* in the nonrecombinative 
case, X/U approaches zero. Recombination reduces the distance for a 
maximum at given time, and thus reduces the apparent mobility, since, 
for nonrecombinative traps with recombination of lifetime tz in other 
centers, the distribution of the mobile carriers subject to trapping is 
simply that for no recombination multiplied by the decay factor e~ x V ° T3 . 
This factor applies because the carriers which arrive at x at whatever 
time have drifted in the conduction band for time x/v . 

The decay constant for the straggle effect is that of the limiting de- 
cay of the tail of the distribution at fixed x after the maximum has passed. 
It is given by the exponent in (158), and is accordingly v v . This result 
follows from (156), since, from (157), x « v t implies cos ~ 1. Real 
v and hence positive v v obtain in this connection. By use of (63), it is 
easily shown that v v for strongly extrinsic material is substantially v „\ + 
v g pi plus either v M , for w » po , or v tp \ , for p » n . 

Integrals of the solutions of (158) over the drift range are evaluated 
in Appendix C. These integrals give 



f v t .U 

F v = (P _1 Apdx = AP dX 

Jo Jo 

= e-*"lit(.K- ± D^sinhiU 2 ± 1)*U + coshi(« 2 ± !)*£/], 



.v l 



(164) 

1 An dx = AiV dX 

o •'o 

= <f !ry (£ - v)U* ± D~ 5 sinh Uk 2 ± 1)>U 



t Use is made of the approximations I (z) ~ Ii{z) ~ (2irz)~k* for |2l large. 
The distributions for large U are substantially proportional and gaussian in 
shape. For nonrecombinative minority-carrier trapping, they are as if the excess 
majority carriers were subject only to drift and diffusion with diffusivity v L/-i£ 3 . 

X Note that, for the maximum, 9 —ir/2 has the sign of k, so that cos Q has the 
opposite sign. 
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as mobile and trapped fractions at given time of carriers initially in- 
jected. For 

pV H t 



? n = (P _1 f An dx = F p - F n , 

Jo 



£ in F p is replaced by ij. Equation (159) serves to verify that, with re- 
combination, F p , F„ and F n all approach zero as U increases indefinitely. 
For nonrecombinative electron trapping, £ = f = (k 2 ± 1) gives Fp = 1, 
as may be expected, and F n = |(1 - *?/f)(l - e~ tu ). Thus, from (145) 
and (152), the trapped fraction approaches T anl /(T tn i + r„„i), the frac- 
tion of the time electrons are trapped, with time constant r ff „iT<ni/ 
(ti„i + t„„i); the mobile fraction approaches the fraction of the time 
electrons are free. For hole trapping, entirely similar results apply. All 
of these results evidently apply for v imaginary as well as real. 

In Fig. 3 are shown continuous minority-carrier distributions for 
imaginary v, in particular, for nonrecombinative trapping of minority 
carriers in the reverse-drift range and also of majority carriers. In the 
former case, an attenuated delta pulse of untrapped mobile carriers, 
which drifts at velocity v , leads a continuous distribution of minority 
carriers that crowds towards the origin as its maximum excursions 
above and below the axis both increase with time. The distribution of 
added minority carriers is negative over part of the drift range after a 
certain time.f In accordance with (164), it approaches a net positive 
delta pulse at the origin of strength (P, the initial strength, times the 
free-time fraction. In the case of majority-carrier trapping, the pulse of 
untrapped carriers increases in strength as it drifts at velocity Vo . This 
augmentation is appreciable with appreciable equilibrium minority- 
carrier concentration; it is negligible in strongly extrinsic material. The 
pulse leads a largely or entirely negative continuous distribution of mi- 
nority carriers, which crowds towards the pulse as its excursion below 
the axis increases with time. J This distribution approaches a negative 
delta pulse, which, with the minority carriers of the augmented pulse, 
gives a net pulse of strength (P. The corresponding distributions of mo- 
bile and trapped majority carriers approach net pulses of strengths 
equal to (P times the free- and trapped-time fractions, respectively. 

These results exhibit, with due allowance for the neglect of diffusion, 

t Negative concentrations of the trapped and majority carriers also occur. 
For this case, ?? — k is positive for p-type material, and, if it is not too large, 
negative An first appears (as for the case of the figure) at the end of the drift 
range for times greater than 2{t) — k)t, since An equals H«n<f _, »>] ir (j| - k — \U) 
for 9 = ir. With sufficiently large 77 - k, negative An first appears within the range. 

X The majority carriers are similarly distributed: For the case of the figure, 
(negative) dimensionless majority-carrier concentration AP or A.V is §AA7 or 
|AP for X = U = 5 and \*AN or If AP for X = U = 10. 
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Fig. 3 — Continuous concentration distributions at different times of mobile 
minority curriers from an injected delta pulse, for drift with minority-carrier 
I rapping in the reverse-drift range and for drift with majority-carrier trapping. 
Equilibrium majority-carrier concentration is taken as 2n, for both cases, and 
nonrecombinative trapping is assumed, with t u /ti large compared with the equi- 
librium ratio of majority-carrier concent rat ion excess to concentration of the 
(mobile) carriers that arc subject to trapping. For minority-carrier trapping, a 
strong inequality for reverse drift accordingly holds, with t„ » 3n = 4t, f = j, 
k = — I, and attenuation of the pulse by the factor e~ u . For majority-carrier 
trapping, r g » %Tt = t, f = f, k = f hold, and the pulse is augmented by the 
factor c« f '. For both cases, £ = — 77 = § or »j = — £ = f holds for electrons or 
holes trapped, so that the coefficients for the minority carriers of the terms in 
J» for the two cases arc respectively \ and —J. 



the essential behavior that would be realized in practice. In each case, 
the initial neutral pulse of mobile carriers will appear essentially as a 
diffuse pulse of trapped and mobile carriers in substantially the propor- 
tions that obtain for the steady state. Minority-carrier trapping in the 
reverse-drift range does not give largely unidirectional drift, while ma- 
jority-carrier trapping results essentially in drift at the ambipolar ve- 
locity i\) . Increase of applied field, however, makes for the idealized 
behavior shown in the figure — in particular, for continuous distributions 
that are negative over increased actual distance at given time. Numerical 
estimate of the effect of diffusion in the reverse-drift case shows that 
negative added-carrier concentrations can occur over appreciable dis- 
tances under conditions that can be realized in practice for reasonable 
values of trapping time.f 

t For Fig. 3, \r u ,,i 2> t = |t ( „i holds for p-type material, and if t is 10~ 5 second, 
then L = vm is about 0.4 cm for p = 2h, and electron and hole mobilities of 1500 
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A simple descriptive interpretation of the crowding of the distributions 
towards the origin in the case of reverse drift is that the drift of added 
carriers, initially in the direction of v , is largely in the opposite direction, 
the direction of drift of majority carriers, after some trapping has taken 
place. For more detailed general interpretations, consider the current 
density Al of added carriers of (19), which is given under the present 
assumptions by 

AI = en„n P <ro 2 I((?/L)(n AP - p AN) 

= ev ((P/L)[e- w « +KCOae) jftff(w - e)]" 1 ^ - 6) (wg) 



• + fV QU sin 0) tan § G 1[0(tt - G)]>, 



from (8), (158) and (161). With (158), this result shows that the Bessel 
functions of order zero are associated with carriers that neutralize the 
charge of trapped carriers or with the trapped carriers themselves, while 
those of order one are associated with the drift of, in effect, carrier pairs. 
The direction of drift of a mobile-carrier distribution considered in its 
entirety depends on the sign of the net A/, or A7 integrated over the 
drift range. By use of (161) and (164), 

,u 

/ AldX = e 2 n n n P <To 2 I((P/L)(noF p - poF n ) 
Jo 

= ev ((?/L)(e- i!:u )[K(K 2 ± l)" 4 sinh |(«" ± l) h U (166) 

+ cosh^U 2 ± 1) 5 C7] 

results. For nonrecombinative electron trapping, this integral reduces to 
ev ((P/L), the initial A/, times t i(l + */f) = v v /(v M + v gn \) in the 
limit for U infinite, from (136), (145) and (159). As may be expected, 
the limiting integral has the sign of v Q or the opposite sign according to 
whether v v is positive or negative. That the distributions ultimately 
crowd towards the origin in the case of reverse drift is established on a 
more quantitative basis by this result. 

and 570 cm 2 volt -1 second -1 (as for silicon 72 ' 73 at 300°K) and for an applied field 
of 100 volt cm -1 , which is reasonable for a filament of resistivity at least a few 
hundred ohm-cm and of area of cross section about 10 -2 cm or less. The correspond- 
ing diffusion distance (Do/)* for t = 10r is 0.054 cm, an order of magnitude smaller 
than the approximate distance 2L of 0.8 cm over which negative added-carrier 
concentrations occur. This difference is greater for larger ti„i . 

t For real v, which implies positive v v and p-type material, this factor is the 
limiting value of X/U for the maximum of the mobile-electron distribution. 
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It follows from (165) and is otherwise evident that A7 is zero at the 
origin.j Also, for imaginary v, the Bessel function with the minus sign 
applies and Al within the drift range is negative for positive v if U is 
not too large. Hence, in the case of reverse drift, the mobile-carrier dis- 
tributions become steeper even for small times. J This conclusion is con- 
sistent with the increase with U of the exponential factor for given small 
9, for which J(f + * cos 6) ~ v v r is negative, and also with the decrease 
with increasing 8 of this factor for given U, which results because of 
k < — 1 . With electron trapping, An becomes small for sufficiently large 
U at some 6. But, with £ > r\, since Aw <-^ Aw is still positive at this 6, 
A/ is still negative, from (165), I being negative for positive r . Thus, 
An becomes negative. As An becomes increasingly negative, electrons 
tend to be released from traps and Aw in turn becomes small [as a zero 
of Jq{\U sin 9) is approached]; but Al is still negative. Thus, Aw be- 
comes negative too. To the right of the location at which Aw equals 
(po/wo — l)Aw < 0, which is the condition A/ = (for negative Aw, 
Aw and An), A/ is positive and the concentrations tend to increase alge- 
braically. A progressive increase with time of the maximum excursion 
of Aw below the axis is associated with the presence of a (time-depend- 
ent) location at which Al changes sign. 

In the case of majority-carrier trapping, the exponential factor for 
given U increases with 9 because of k > 1; and f — k being negative, § 
this factor is smaller or greater than unity respectively for 9 small or 
^ v. Thus, as U increases, the continuous minority-carrier distribu- 
tion, which is negative for all 9 and U, crowds towards the pulse, which 
drifts at velocity v , and both increase hi amplitude. For electron trap- 
ping in general, An > holds at (and near) the pulse. If the material 
is sufficiently strongly n-type, this An is largely compensated entirely by 
negative An. The pulse, in effect, removes electrons from the semicon- 
ductor and transfers them to traps, and its strength remains substanti- 
ally constant ; the continuous contribution to A/ is negligible because of 
large t. In less strongly n-type material, however, neutrality is main- 
tained at the pulse in part by negative Aw and correspondingly more 
negative Aw, which makes for less trapping for given pulse strength. 
But holes removed from the semiconductor drift at velocity v and thus 
(with an equal number of electrons) cause a progressive increase in the 
strength of the pulse; a Aw obtains that does not depend explicitly on 
??o . The Al of the pulse likewise increases indefinitely. The net or inte- 

t Thus, at the origin, An/Ap maintains the value n /po ■ 

% Note also that (158) gives, for small U > 0, AP or AN at the origin equal to 
£ — k or n — k times (1 — vj), which increases with I in the case of reverse drift. 

§ Note that (159) may be written as (f -f- «)(f — k) = — 1 for nonrecombinative 
trapping and imaginary v. 
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grated Al, however, increases to a limiting value, as given by the factor 
Vv/{vt,a -\- Vgni), which is greater than unity in the present case. The 
largest excursions of the continuous mobile-carrier distributions below 
the axis occur, of course, at the pulse, at which there is a discontinuity 
in Al with change of sign. 

In cases of real v, the exponential factor is less than unity and de- 
creases as U increases.! With the monotonically increasing modified 
Bessel functions and Al everywhere in the direction of v , the continu- 
ous distributions, and, from (165), Al as well, ultimately possess 
relative maxima that drift in this direction with a common velocity. 
While the maxima occur substantially together, the distribution of the 
trapped minority carriers lags, consistent with the Bessel function of 
order zero dominating that of order one; the mobile minority-carrier 
distribution ultimately results entirely from carriers released from traps. 

The case of "critical trapping," the borderline case between cases of 
real and imaginary v, is one whose analysis furnishes further qualitative 
insight. For it, v as denned in (143) is zero with vui j* v lp i , and the 
condition is the equality corresponding to (162). The previous notation 
can be used for this case by choosing the time unit r arbitrarily and 
noting that (144) through (154) then all apply if the terms with the 
double signs are omitted wherever they occur. A calculation essentially 
similar to that given in Appendix B provides the solution, which, for 
the initial delta pulse of (P carrier pairs per unit area at the origin, is 

A V = ((P/L){exp [-|(r - k)U]-8(U - X) 

+ }({ - t) exp [kX - Kf + k)U]-1[X(U - X)}\ 

= <?\ exp [-Uvtm + v tpl )(l + (n + po)/fJli*)t]-S(vot - x) (167) 

+ \po(v la i + vt P i)/vi&li*\ exp [-(v<ffli*)~ l 

•(pot>t,,i + nt,v tp i)x — t/i(r tn i + T tP i)]-l[x(vdt - x)]\, 

An = J((P/L)(* - v) exp [kX - \{k + S)U\-1[X(U - X)] 

= <P[(po - noXriai + vw)/vtfli*\ exp [- (i/ofti*)" 1 

•(PflVtnl + noVtpl)x - t/\{Ttnl + T tp i)]-l[x(vat - as)]. 

For An, ^ is replaced by rj in the first expression for Ap, or n and p as 

t In these cases, f + k has the sign of v,. and is positive, from (156), and (159) 
may be written for nonrecombinative trapping as (f + «)(f — «) = 1. 
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well as subscripts n and p are interchanged in the second. The concen- 
tration increments in terms of the dimensionless quantities are, of course, 
actually independent of r. In the second expressions, quantities involved 
have been written in forms that apply under the condition for critical 
trapping. It follows readily from these equations that the mobile and 
trapped fractions are similar to those of (164) except that (k 2 ± 1)' is 
replaced by k; they are independent of r. For nonrecombinativc trapping, 
the mobile and trapped fractions of carriers respectively approach, as 
before, the free- and trapped-time fractions, since k equals ( — f ).f 

As the second expression for An shows, critical trapping is a case of 
minority-carrier trapping. The continuous distributions are proportional 
and are equal to products of exponential decays with distance and with 
time. The amplitudes of the distributions at the origin are larger and 
the decay with distance is sharper the smaller is v . The time decay 
results from recombination, with which neither T t „\ nor T tp \ is infinite. 
For nonrecombinativc trapping, the condition of critical trapping is the 
same as that of zero drift, and exponential distributions are established 
progressively over the variable range v t that otherwise do not change 
with time. It is well to note that the general case of zero v v is a case of 
imaginary v, with v 1 equal to (— 4Ai). Furthermore, for trapping in 
intrinsic material, for which v,. is not defined, it can be shown that the 
distributions are all identically zero except at the origin. As may be 
expected for this case of no diffusion, the initial pulse results simply in 
pulses for the concentration increments that change as trapping and 
recombination proceed. 
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t Note thai f = u- -f lA,)i and * < (l hold. 
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APPENDIX A 

Derivation of Two- Sided Laplace Transforms 
The transforms 

2 {lxJo [(ff ~ X * )i]m d( *} = (S2 ± l) ~* F[ ~ ( * 2 ± ^ 

and ( 168 ) 

8 lf(X) ± £W - Z 2 )"^M(/3 2 - X 2 ) 5 ]/(/3)d/3| 

= F[-(.s 2 ±l)*] 

entail restrictions on f(X) and on the transform variable s for conver- 
gence of the integrals. In accordance with the definition of (141), the 
first transform is 

J- m J\ Jo 

J-co J-=oJ0 (1(;9) 

= f [/(/3) + /(-«! f r° [(A' - M ^ «& # 



+ 



Jo J—pio 



as obtained by changing the order of integration and changing from 
variable to — in the contribution from the range of negative 0. The 
first of the integrals with respect to X isf 



= r m (m 2 + rt H f° (m) exp [s(m 2 + /3 2 ) 5 ] d/» (170) 

Jo J o 

- (« 2 ±l)- 5 exp[-(s 2 ±l)^]. 

For convergence with the upper function and sign, s is taken as (real 



t Watson, 93 p. 416, Equation (2). 
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and) negative; for convergence with the lower function and sign, the 
real part of s is to be less than — 1. The second is a Gegcnbauer integral: 
With 

X = — /3 cos 4>, 



(171) 





z sin >A = ( | j ft 




z cos \f/ = —isp, 




z = i(s 2 ± 1) 1 I3, 




cot*= -Q S , 


the result 




J-a h 


- 2 ) 5 ]<r x dX 



= (3 T J n (z sin <f> sin +) e " C08< » c ° s * H [ n </ (172 ) 
Jo 

= 2(s 2 ± l)- 4 sinh [(s 2 ± l) s /3] 

follows. f By substituting from (170) and (172) in (1G9), the first trans- 
form of (168) is established. The second transform is established through 

= r/(X)exp[( S 2 ±l) ! X]r/\ (173) 

J — oo 

= F[-(s 2 ±l) 5 ]. 

The original expression reduces to the form on the left upon integrating 
by parts for/(/3) such that /o(/3)/(/3) or J o (0)f(f3) is zero for = oo. 
The result then follows by application of the first transform and another 
integration by parts, which entails /(X) exp [(s 2 ± 1) 5 X] equal to zero 
for positively and negatively infinite X, a condition that subsumes the 
preceding one. 



t Watson 93 p. 379, Equation (1). 
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APPENDIX B 

Solutions for Drift of an Injected Pulse 

The dimensionless concentration AP = Ap/((P/L) is the inverse 
Laplace transform of AP given by (154). It may be evaluated by use of 

r 1 exp (aV - kUa) = W*<T l exp [- (X - hU) 2 /*a*} (174) 

and the similar formula obtained by replacing V by its negative in con- 
junction with (168) and certain elementary rules. Equations (168) and 
(174) give 

rV ± D" ! exp [oris 2 ±1)4- Wis 2 ± l) J ] 

= Itt-V 1 f h ; [(/3 2 - X 2 )*] exp [ - {fi - W)V±a\ d0, 

J.Y «/0 

r 1 exp [a 2 (s 2 ± 1) + JC7(s 2 ± 1)*] (175) 

- frr-V 1 [exp [- (X - if/) 2 /4a 2 ] ± jf tflS 8 - X 2 ) H 

• £ [(/3 2 - X 2 )'] exp [- OS - JW'/VI #} . 

The exponents that occur in (154) may be transformed in accordance 
with 

oV - |tT(« + r) +H'[(s-.) 2 ± 1] J 

= a 2 (/c 2 =F 1) - J(« + t)V + (2^« " 2 f/ )(s - «) (176) 

+ a 2 [(s - K ) 2 ± 1] + §#[(« " «) 2 ± U 1 

and the similar identity that holds with minus signs before the radicals. 
Consideration of the transformations that change the exponent in the 
transforms of (175) to that given by (176) shows that the inverse trans- 
forms or solutions sought contain exp [o 2 (k =F 1) — J(k + t)U + «X] 
as factor, in which exp (kX) results from replacing 8 by s - k. They con- 
tain also X + 2a k - \U in place of X in (175), because of the term 
(2a k - £E/)(s - k) in (176). The solutions are obtained by straight- 
forward application of these results to (154); corresponding to the s in 
the factor s — £ of the second terms, there are contributions obtained 
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by differentiating with respect to X the expressions that originate from 
the first equation of (175). The solutions are: 

AP = JrV 1 {exp [aV =F 1) - |(f + k)U + kX]\ 
Yexp [-(X + 2a\ - tf) 2 /^ 2 ] - *({ - «) 

f / ;{[/i 2 -(X + 2a 2 K-iC/) 2 ] 1 } 

J:r-l-2/i2*— W i/ n 

t/) 2 /4a 2 ]| rf/3 

(177) 



'X+2a 2 K-}(/ '/o 

{exp [-(/3 + hU)*/4a*\ - exp [-(/3 - hU) 2 /4a 2 ]\ <tf 



=F 5 / [0* - (X + 2a\- - it/) 2 ]- 5 



■\ l {{f3 2 - (X + 2a\ - W)*?} 
J i 

•{(X + 2o 2 k - \V - 0) exp [-(0 + Jt/)'/^ 

- (X + 2a 2 « - |C7 + 0) exp [-(/3 - ^f/) 2 /4a 2 ]| r/A 

Replacing £ by t? in (177) gives the solutions for AN, and those for AN 
are accordingly 

AN = |T~*a -1 (ij - £){exp [aV =F l ) - |(f + *)£/ + *X]} 

• f 7 ; j[/3 2 - (X + 2a 2 K - W)V) (178) 

•[exp [-(/3 + H r ) 2 /4a 2 ] - exp [-(0 - |l7) , /4o t |} d0. 

The corresponding limiting solutions of (155) involve the step-function 
factor as result of the requirement that, for contributions in the limit 
of zero a, the gaussian factors in the integrands of (177) and (178) be 
centered at values within the range of integration. 

APPENDIX C 

Integrals Over the Drift Range 

Consider first evaluation of /*'„ , which may be written as 

K = 1(4- i)r[exp(-l.f(1] 

(179) 

• [ [exp (-U-C cos 6)] ! j (H' sin G) sin B (/6 
Jo Jo 
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from (157) and (158). The transformations 

,WB *-VW (180) 

z = %(k ± \fu, 

cot \1/ = ( . ) 
\ikJ 

result in the single form, 

fn = \ ($ - y)U[exp (-KU)] f Q [exp (iz cos 9 cos *)] 

• J (2 sin 9 sin ^) sin 9 d9; 

the Gegenbauer's integral! reduces to (2ir/z)' l J i (z) = 2z _1 sin z, and, 
with the definitions of z, the respective expressions for F n of (164) follow. 
The contribution to F p from the Bessel function of order zero is clearly 
[(£ _ k)/(£ — ri)]F n , while the contribution from the delta pulse at the 
limit of the drift range is exp[§(« - f) #]• The contribution 

±Mexp[-i(f + «)Ul} 

•f \[x(u - x)rW («a)]JM(x(17 - x))*]<zx 

- ±117 [exp (-UU)] f [exp (-IkU cos 9)] 



/ r 1 (|C/sin9)(l - cos 9) d9 



(182) 



- MW[«p(-*rl7)] 

• / [exp (iz cos 9 cos 4t)]J\(z sin 9 sin ^)(1 - cos 9) de, 

the last form for which follows from (180), remains to be evaluated. A 
method due to Pollak 94 involves first writing the integral in the last 
form, to be denoted by #, in terms of 

11 = z cos \p, V = z sin \j/. (183) 

t Watson 93 p. 379, Equation (1). 
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Then, the operator (1 — z'd/dOl) applied to g replaces, in effect, 1 — cos 
in the integrand by sin 2 0; it gives a Gegenbauer integral t which reduces 
to (2t/z)* sin \p J\(z) = 2Vz~ 2 (z~ l sins — cosz). Since this is 
-i[exp(-i c \i)]d(3 exp i € \l)/d e \L, 

3=3 1^=0 exp (-ill) + i2V [exp (-«Ol)] 

[^ _, _, . (184) 

• / z ~(z sin z — cos z) exp (ill) dot 
•'o 

follows, in which z is (Ol 2 + c 2 ) i . The integral in this equation may be 
evaluated by writing sin z and cos z in terms of exponentials and intro- 
ducing 01 + (Ol 2 + c U 2 ) i and 01 - (Ol 2 + I) 2 )* as variables of integration. 
It is found to equal 

Mexp (;0l)]{[0l 2 + V 2 + 0l(0l 2 + V 2 )*]' 1 exp [?:(0l 2 + v 2 ) 1 ') 

+ [oi 2 + v 2 - oi(or + v 2 ) ] ] 1 exp [-?"(oi 2 + v 2 y\\ - iv~ 2 cos v. 

With # I 01=0 equal to 2 c lT 1 (l - cos V) from an integral of Sonine,t^J is 
thus found in terms of 01 and V and, by use of (180), in terms of k 
and U. The respective expressions for T^p of (164) then follow. 
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